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COVERAGE OF JOURNAL 

THE aim of this journal of abstracts is to give complete coverage of papers in the field of statistical 
theory and new contributions to statistical method. Papers which report only applications or examples 
of existing statistical theory and method will not be included. There are approximately two hundred 
and fifty journals published in various parts of the world which are wholly or partly devoted to 
the field of statistical theory and method and which will be brought within the scope of this journal of 
abstracts. A complete list of journals covered is printed in the annual Index Supplement. In the case of the 
following journals, however, being those which are wholly devoted to statistical theory—all contributions, 
whether a paper, note or miscellanea, will be abstracted : 


Annals of Mathematical Statistics 
Biometrika 

Journal, Royal Statistical Society (Series B) 
Bulletin of Mathematical Statistics 

Annals, Institute of Statistical Mathematics 
Sankhya (Series A) 


Within the larger group of journals, which are not wholly devoted to statistical theory and method, there 
are some journals which have the vast majority of their contributions in this field. These journals, there- 
fore, will be abstracted on a virtually complete basis : 


Biometrics 

Metrika 

Metron 

Review, International Statistical Institute 
Technometrics 

Sankhya (Series B) 


If any reader of this journal discovers a paper which happens to have been overlooked, the General Editor 
. will be pleased to be informed so that the appropriate abstract can be made: always provided that the 
date of publication is after 1st October 1958, when the abstracting for this journal commenced. 

In addition to the ordinary journals, there are two kinds of publication which fall within the scope 
of this journal of abstracts. ‘They are the experiment and other research station reports—which occur 
particularly in the North American region—and the reports of conferences, symposia and seminars. 
Whilst these latter may be included in the book review sections of journals it is unusual for any individual 
contribution to be noted at any length. These publications are, in effect, special collections of papers 
and for this reason the appropriate arrangements will be made for them to be included in this journal. 
By the same token, abstracts of papers given at conferences and reproduced in an appropriate journal will 
be disregarded until the definitive publication is available. 


FORMAT OF JOURNAL 
The abstracts will be up to 400 words long—the recommendation of UNESCO for the “long ” 


abstract service: they will be in the English language although the language of the original paper will 
be indicated on the abstract together with the name of the abstractor. In addition, the address of the 
author(s) will be given in sufficient detail to facilitate contact in order to obtain further detail or request 


an offprint. Suitable indexes will be provided annually in a supplement. 
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The scheme of classification has been developed upon lines that will facilitate the transfer to punched 
cards of the code numbers allocated to each abstract: to allow for future development it is suggested that use 
is made of 4-column fields. Each abstract will have two classification numbers: the primary number in 
heavy type to indicate the basic topic of the paper and the secondary number in brackets to take account of 
the most important cross-reference. The sheets of the journal are colour-coded according to the twelve main 
sections of the classification and it should be noted that it is the main section number of the primary 
classification which determines the colour code for each abstract. It is believed that this method of 
colour-coding the pages will provide a distinctive visual aid to the identification of abstracts both when 
the journal is in bound form or dismantled and filed on cards or in binders. The format and simple 
binding allows of the following alternative treatments by users of the journal : 

(a) Leave intact as a shelf-periodical. 

(b) Split and filed in page form according to the main sections of the classification. 

(c) Split and guillotine (single cut) each page ready for : 
(i) pasting onto standard index cards. 
(ii) filing in loose-leaf or other binders—for which the appropriate holes are punched. 


It will also be possible for those users who need a completely referenced file to insert skeleton cards or 
sheets according to the secondary classification number provided on each abstract. 
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Volume 2. CORRECTIONS 


2/447 1.4 Palasti Author’s christian name is, Ilona. 
Author Index: Vol. 2, No. 1 Filippo, E; delete and add, Emmanuelli, F. 
a Bs ‘3 No. 2 Palasti; delete initial and insert the name, Ilona. 


Note The symposium from which papers are reported at 2/500, 501-2, 503, 510 and 
~--4545 was held at Leyden and published in Amsterdam (North Holland 
Publishing Co.) 
The Congress from which papers are reported at 2/433, 442, 444, 445, 447, 
451, 452, 453, 454, 471, 479, 482, 491, 492, 511, 517, 518, 527, 558, 571, 576, 
579-80, 587, 592, 741 and 767 was held at Budapest and published there by 
Akadémiai Kiado. 


CANDLER, W. (University of New England, Australia) 


0.8 (0.7) 


A “ short-cut ’” method for the complete solution of game theory and 


feed-mix problems—ZIn English 


Econometrica (1960) 28, 618-634 (6 references, 3 tables, 4 figures) 


The author states that this “ short-cut’ method is an 
extension of the graphical short-cut methods presented 
by Waugh & Burrows [Econometrica (1955) 23, 18-29] 
and by Boles [Journal of Farm Economics (1956) 38, 981]. 

Commencing by defining the terms in the title of his 
paper, the author remarks that he is using the term 
“complete solution” to imply that with this method 
it is not necessary to employ the simplex method after 
the completion of the “short-cut” method. Even in 
the problem where one is dealing with a large number 
of columns and rows, a solution can be reached by use 
of the “ short-cut ’’. Two-person zero-sum games are 
used in the discussion of game theory problems: pro- 
gramming problems with all positive resource supplies 
and with either positive or negative prices are used in 
the feed-mix problems. 

Following a description of the notation employed, 
the short-cut method is discussed and illustrated by an 
example dealing with game theory. Concise tables are 
given both for computing the successive steps in the 
“‘ short-cut’ solution to the problem posed and to 
illustrate the algebraic representation of this table. 
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DINKELBACH, W. & STEFFENS, F. 


(K6In University) 


After a discussion of strategy, an outline is given of 
various difficulties which may be encountered and the 
manner in which strategies can be both improved and 
eliminated. The feed-mix problem used to illustrate the 
use of this method is then described and its solution 
discussed. 

In conclusion the author remarks that it is now 
apparent that, at its present level of development, the 
short-cut programme appears to require a higher level 
of decision-making than is usual in the case where the 
simplex method is employed. He considers therefore, 
that the simplex method is still superior as a method 
for use by clerical assistants and high speed computors. 
The short-cut method can usefully be employed when 
computers are not available. 

The author’s main object in writing this paper has, 
he states, been to develop a complete non-simplex 
method for the solutions of a wide range of problems 
in the field of programming: further research is intended 
into short-cut methods both of general and parametric 
programming problems. 


(W. R. Buckland) 


0.8 (-.-) 


The solution of some decision problems by mixed integer linear programming—In German 
Unternehmensforschung (1961) 5, 3-14 (7 references) 


In addition to usual quantitative restrictions in linear 
programming the following qualitative restrictions are 
considered. Certain subsets only of a set of variables 
can be non-zero in the final solution of the linear 
programme. 

The following four models are formulated as mixed 
integer linear programming problems: 


(i) The final solution consists of at most exactly 
one variable. 

(ii) The final solution consists of at most exactly 
k (k <n) variables. 

(iii) Each activity P; is divided in n; sub-activities 
P;,; the corresponding levels of the activities 
are x, and x,; 7 = 1,....m,r = 1, .-., Mj The 
final solution contains at most one x, of P;; 
the total number of variables x, is at most 
k (k<n). 

(iv) The set of variables is divided in two subsets. 
The final solution contains only variables of one 


of these subsets. 
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It is proposed to solve the models by the method of 
Gomory, see “An algorithm of the mixed integer 
problem ” [RAND Corp. Papers (1960) p. 1885]. This 
method is explained in the second part of the paper. 
An example of production planning demonstrates how 
the solution of a linear programming is modified by 
qualitative restrictions as given in the third model. 


(W. Piesch) 
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DANSKIN, J. M. (Rutgers University, New Jersey, USA) 


0.7 (10.5) 


Calculus of variations and operations research—In English 
ARL Technical Note 60-151 (1960), Aero. Res. Lab., Wright-Patterson Air Force Base. 


iv+19 pp. (7 references) 


The author discusses the role of mathematics, especially 
the calculus of variations, in operational research, and 
exhibits some problems of operational research. The 
first problem is that of distributing a given amount of 
available resources among various alternatives so as to 
obtain the maximum return. If the number n of alter- 
natives is finite, the problem may be posed as follows: 
maximise 24(x7,)d..4°-x2) Subject 'to. 0%, =X; 7,20: 
where x,/X denotes the proportion of the resources 
assigned to the ith alternative. If F is strictly convex, 
the solution is trivial; one selects the largest of the 
Comersvalest# (1, 0,,..:; 0), F(O,. 1, 0,-.... 0), ..., FO, 
meres tS Mie, 25 x) SD f.(x;),. wherestheaf;aré 
strictly concave, a unique solution is given. If the f; 
are concave, but not strictly concave, the solution is 
no longer unique. A solution is given for the case in 
which the f; are S-shaped; first convex, then concave. 
These solutions of the finite problem are generalised 
by the use of calculus of variations to solutions of the 
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DANSKIN, J. M. 


(Rutgers University, New Jersey, USA) 
Calculus of variations and operations research—In English 


analogous infinite problem; stated as follows: maximise 
{r [x (2), t] dt subject to [> (t) dt = X, x (t)=0, where 


the x’s are measurable functions of 1, with f(x, t) and 
f. (x, 1) continuous and f, (x, 2)>0. 

Examples of the application of the results are given, 
together with game theory and topology, to the solution 
of various problems. Many military problems amount 
to zero-sum, two-person games in which each side 
disposes his forces according to some distribution and 
then awaits events. The result is embodied in a payoff 
function 2 (x, y). The game is said to have a solution 
if there exist a value vy and optimal strategies x° and 
y® such that (x, y°)2v9 for all x and 2(x°, y)2v, 
for all y. Military problems solved by these methods 
include: firstly, a stockpiling problem in which it is 
desired to stockpile a certain item for a contingency, 
for example war, the latter happening only once; 


continued 


0.7 (10.5) 
continued 


ARL Technical Note 60-151 (1960), Aero. Res. Lab., Wright-Patterson Air Force Base. 


iv-+19 pp. (7 references) 


secondly, a two-machine-gun duel problem in which a 
fighter closes on a bomber at a constant rate, so that the 
rate of fire of his machine gun may be described as a 
function x (r) of the range. Similarly the bomber has a 
machine gun, with rate y(r). The guns have limited 


total amounts of ammunition, hence [xo Oh =X 


> (r) dr = Y; thirdly, a problem in which each side 


places military forces linearly, whose value at point ¢ is 
v(t). At each point the greater force wins. Similar 
games have mixed strategies; since if the weaker player 
discloses his strategy he loses all. Suggestions are given 
for solving nontrivial min-max problems in which 
player x chooses his strategy and player y, knowing 
x’s strategy, chooses his own to minimise the payoff 
LOEX. 


(Hols Harter) 
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FERSCHL, F. (Bundeskammer der Gewerblichen Wirtschaft, Vienna) 0.8 (-.-) 
A calculation scheme for the simplex method—ZJn German 
Unternehmensforschung (1961) 5, 32-41 (2 tables, 1 diagram) 


The author starts by describing the normal form of 
linear programming, generally known as the ordinary 
maximisation problem. Thereafter a detailed calcula- 
tion scheme is given to solve the problems mentioned. 
It differs from the usual ones, as given for instance by 
Charnes, Cooper ‘and Henderson [An Introduction to 
Linear Programming (1953) New York: Wiley] by an 
auxiliary calculation scheme which facilitates the com- 
putation and the storage of provisional results. The 
calculation is carefully explained by a flow diagram. It 
does not make use of matrix calculus and may be 
applied by non-mathematicians. Two examples illus- 
trate the scheme. 

In the appendix the author mentions some methods 
by which the number of computing errors can be reduced 
and by which some more general linear programmes 
can be transformed into those of the normal form. 


(W. Dinkelbach) 
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FINCH, P. D. (London School of Economics) 0.6 (10.6) 


On the covariance determinants of moving average and autoregression 
models—ZIn English 


Biometrika (1960) 47, 194-195 (2 references) 


The weakly stationary process defined by 


p q 
yi A,X,» = » Deoray 


r=0 r=0 
where the {e,} are uncorrelated errors of unit variance, 
has the covariance determinant YX, = | &(x,x,) || for 
Xo, +++) X, The author finds lim X, as noo in terms 
of the roots of La,z” = 0 and 2 b,z" = 0 and shows it 
to have the same value when the a’s and b’s are inter- 
changed. 

The result is obtained by expressing & (x,x,) in terms 
of the spectral density function of the process, noting 
that these may be regarded as Toeplitz forms and hence 
applying a theorem of Grenander & Rosenblatt on the 
limit of the determinant of such forms. 


(D. E. Barton) 
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FORTET, R. (Sorbonne University, France) 


0.0 (0.8) 


Boolean algebra and its applications in operational research—In French 
Cahiers Cent. Etudes Recherche Opérat. (1959) 1, (10 references) 


The author introduces the Boolean algebra U consisting 
of the elements 0 and 1 and the operations of dis- 
junction, conjunction and negation. A variable taking 
its values in U is called a Boolean variable. A function 
f of n independent Boolean variables x,, x2, ..., Xm 
f (U) is called a Boolean function. Generally a Boolean 
problem is to find amongst the 2” sets of values of 
(x1, X2, -.., X,) those that satisfy a given condition C. 
Two special cases of the condition C are treated: 


(a) C is.an algebraic equation: P (x,, x»,...-;.x,) = 0, 
P being an entire algebraic function, that is a 
function of the form 

Pr Yap xixG 
i ik; 
where 2 is taken in the ordinary, non-Boolean 
sense, II denotes Boolean multiplication, the 
a; are constants and the elements of U are inter- 
preted as ordinary numbers. 

(i) -C is-a Boolean equation: ((%,> x; --- 

f being a Boolean function. 


, X,) = 0, 
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HELLMICH, K. (Graz) 


The route of shortest length or duration—Jn German 


Methods of solution are discussed with the aid of 
simple examples and some practical applications are 
given; firstly to the four-colour problem, secondly to a 
machine problem and also to the problem of finding 
elementary Hamiltonian circuits in a given graph. 

Finally, the problem of finding a set of values which 
satisfies C and which minimises a given entire algebraic 
function is discussed and the author illustrates this with 
an a problem from the field of operational research. 
See also abstract No. 2/430, 0.0. : 


(F. W. Steutel) 


0.8 (-.-) 


Math. Tech. Wirtschaft (1960) 4, 166-175 (7 references, 1 table) 


The author gives a solution of the “‘ travelling salesman 
problem ”’: n given cities are considered to be connected 
by a route of shortest length ending at its starting point. 
The method of solving this problem given by the author 
can also be used for asymmetrical matrices of the 
distances. Moreover, the parts of the route given in 
advance can be taken into account. 

The method is illustrated by a numerical example. 
Firstly, a known method of reduction of the matrix is 
demonstrated, leading to groups of neighbouring cities, 
which are each connected by a closed route. Thereby 
the reduced matrix is divided in submatrices, called 
cycle matrices and distance matrices: the cycle matrices 
represent the preliminary subcycles obtained by the 
reduction. The distance matrices contain all possible 
distances between cities of different subcycles. 

A systematic procedure is given to connect the sub- 
cycles to one shortest route by appropriate solution of 
distances. For this purpose all distances are ordered. 
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Beginning with the shortest distance subsets are system- 
atically formed and investigated; whether they can be 
used for a shortest cycle or not. The given numerical 
example is computed by the selection procedure and a 
computing programme is given for the general case. 


(F. Ferschl) 
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HOUTHAKER, H. S. (Stanford Univ., California, USA) 


0.8 (-.-) 


The capacity method of quadratic programming—ZIn English 
Econometrica (1960) 28, 62-96 (6 references, 3 figures) 


“* Quadratic programming, as considered in this paper, 
is the art of finding non-negative values of several 
variables x,, ..., X, which maximise a polynomial of the 
second degree # (x,, .-., X,) While satisfying given con- 
straints. The latter constraints will be assumed to be 
linear.” After discussing the facts that quadratic 
programming is of interest in economics because of its 
flexibility and realism and does not assume, as does 
linear programming, that the unit cost of each activity 
or process and the availabilities of those inputs or 
requirements for those outputs subject to constraint, 
are all fixed, the author states that possibilities exist 
by postulating appropriate linear dependencies, to expand 
linear programming into quadratic programming. In 
this paper a method is presented for assigning non- 
negative values to several variables so as to maximise 
a given quadratic objective function which satisfies 
given linear constraints with non-negative coefficients. 
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KOSAMBI, D. D. (Poona University) 


An application of stochastic convergence—In English 


The method used by the author is parametric and 
iterative; it leads to a solution in a finite number of 
steps. 

Computational aspects are dealt with at some length 
and include the opinion that, “another problem on 
which more experience is needed is that of numerical 
accuracy, particularly the effect of round-off errors ”’. | 

Finally the author says it will be interesting to compare 
the efficiency of the capacity method with that of other 
methods of quadratic programming. 

An appendix gives computing directions for quadratic | 
programming in the case of non-negative constraints. 


(W. R. Buckland) 


| 
4 
0.4 (0.1) | 


J. Indian Soc. Agric. Statist. (1959) 11, 58-72 (7 references) e . 


The author considers the sequence 
a,= yin ’— ) p “logp 
nsx psx 
where the first summation is over all integers and the 
second is over the primes and where a is a fixed constant. 
Probabilistic techniques are used to prove that limit 


x oO 
a, exists for ¢>4. From this the author deduces that 
if a Riemann zeta function is defined by the equation 


f= Sv n-i(s sot it; o>1) 
1 


and that if it is continued analytically, it has all its 


non-trivial zeros on the vertical line o = 4. 


(S. R. S. Varadhan) 
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KUNZI, H. P. (Universitat Ziirich) 
Non-linear programming—In German 
Ablaufs Planungsforschung (1961) 2, 8-14 


This article is a paper presented at a recent meeting 
in Minster of the AKOR. The author starts with a 
short review on linear programming: he explains the 
mathematical problems and the ideas of two methods 
of solution, the simplex-method of Dantzig and the 
gradient-method of Frisch. 

A problem in production planning in which the prices 
depend on the output is given as an example of a non- 
linear programme. The problem of the quadratic 
programme is formulated. Before discussing the methods 
of solution for a quadratic programme the theorem of 
Kuhn-Tucker is given, which concerns the extrema of 
convex functions under convex side conditions for non- 
negative variables. The necessary and sufficient con- 
ditions of the Kuhn-Tucker theorem can be applied in 
the special case of the quadratic programming provided 
the quadratic objective function is semi-definite. The 
methods of solution of Barankin and Dorfman are 
based on this theorem. Barankin’s method is explained 
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LAND, Ailsa H. & DOIG, Alison G. 


(London School of Economics) 


0.8 (-.-) 


in detail and some connections with the simplex-method 
are pointed out. Afterwards Beale’s method [Proc. 
Camb. Phil. Soc. (1954) 50, 513-523] is considered; this 
also follows the simplex-method. Beale starts with a 
basic solution, but uses variables without restrictions on 
their signs. 

Finally the paper gives a short survey on other 
methods of non-linear programming. Houthakker’s 
capacity-method [Econometrica (1960) 28, 62-96; ab- 
stracted in this journal No. 2/607, 0.8] is fully discussed 
and, furthermore, is a method which starts with an 
unrestricted extremum and varies this by adding side- 


conditions successively: the gradient methods are 
quoted. 
Note: The references to this paper are to appear 


in a later issue of Ablaufs Planungsforschung. 


(F. Ferschl) 


0.8 (0.6) 


An automatic method of solving discrete programming problems—In English 
Econometrica (1960) 28, 497-520 (8 references, 9 figures) 


A generalisation of the classical linear programming 
problem would be to relax the continuity condition 
on the non-negative variables which are subject to a 
system of linear inequalities. This paper presents a 
numerical algorithm for the solution of programming 
problems where all or some of the variables are discrete, 
In some problems the discrete variable constraints are 
both significant and costly. The algorithm in the paper 
is simple and lends itself to automatic computations. 
An alternative approach was given by Gomory [Bull. 
Amer. Math. Soc. (1958) 64, 275-278] and extended by 
Beale [SRG Princeton (1958) Tech. Rep. 19]. 

After sections on the formulation of the problem 
and the description of the method a numerical example 
is given. Two computational routines, one based upon 
the solution of many simple linear programmes and the 
other on a parametric linear programming, are given in 
Appendix I. In Appendix II a solution is given to the 
non-negative integer problem originally stated in a 
paper by Markowitz & Manne [Econometrica (1957) 
25, 84-110]. 
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The main method of the paper is directed towards 
overcoming the difficulty that because a set of feasible 
solutions to the discrete variable problems is not convex, 
most methods demonstrate only the existence of a 
local optimum. The method proposed here is to make 
systematic parallel shifts in the functional hyperplane 
until a print within the ordinary linear programming 
set is found which has integral co-ordinates in the 
specified dimensions. Although this is obvious in 
principle, the difficulty of visualising a hyperplane in 
n-dimensions with a point whose x co-ordinates are all 
integers, makes it necessary to devise rules to switch 
attention between possible points so as not to ignore 
any integer point. 


(W. R. Buckland) 
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TINTNER, G. (Iowa State University, Ames) 
A note on stochastic linear programming—In English 


0.8 (0.6) 


Econometrica (1960) 28, 490-495 (15 references, 3 tables) 


An ordinary linear programming problem has the 
system of inequalities, the availabilities and coefficients 
of the objective function as assumed or given quantities. 
If these parameters become random variables in order 
to improve the realisation of the model, then the problem 
becomes one of stochastic linear programming. 

A passive approach to this problem is to approximate 
the distribution of the objective function and to let 
the decisions be based thereon. The active approach 
is to let the decision variables be the amount of resources 
devoted to the various activities within the scope of the 
programme. This note investigates the numerical 
methods for optimising the various distributions of 
resources. Three examples are given. 


(W. R. Buckland) 
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VOTAW, D. F., Jr. (Yale University) 
Statistical programming—ZJn English 


0.8 (-.-) 


Ann. Math. Statist. (1960) 31, 1077-1083 (8 references) 


The author states, “A ‘statistical programming’ 
problem is encountered when the information about 
one or more constants in a programming problem is 
statistical.”” Two programming problems are mentioned, 
firstly, a special case of the classical Hitchcock-Koopmans 
transportation problem; and secondly, the personnel 
assignment problem. The first problem is to find an 
nxn matrix (x;;) of real numbers such that 


n 
¢;;X;; is minimised, subject to the constraints 


x..>0 where the 


yi 
j=l 
matrix, (c;,;), is given beforehand. The assign- 


ment problem is similar, except that (i) is to be 
maximised, and the nonlinear constraint 

(iii) each x;; = 0 or 1 is introduced. 

Of the two, the assignment problem is discussed in 
detail. In this case (iii) dictates that there are n! admis- 
sible matrices (x;;). The problem is to find a permuta- 
tion (jj, .-., j,) such that 
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(iv) Cy +...+¢nj, = max (C,;,+...4+¢,;,) where 


(Jt. «9 Jn) 
(j1, «++» j,) denotes any permutation of (1, ..., 7). 


The sum (c,;,+...+¢,;,) is termed a “ permuta- 
tion sum ”’. 

The situation treated in this paper arises if the values 
of the c;;’s are not all known but information regarding 
them is available. The notion of “ programming by 
estimation ”’ is introduced. Estimates of the n! permuta- 
tion sums based on observed values are considered and 
the permutation corresponding to the largest such 
estimate is selected. 

Programming by estimation is compared with purely 
random programming, in which a uniform probability 
1/n! is assigned to each permutation. A theorem is 
stated and proved to the effect that under rather general 
conditions, programming by estimation is uniformly at 
least as good as purely random programming. Under 
slightly more restrictive conditions, programming by 
estimation is uniformly better. An example where 
random programming is better is given. A more 
general form of the assignment problem is also briefly 
discussed. 


(R. M. Durstine) 
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BARTON, D. E., DAVID, Florence N. & FIX, Evelyn (Univ. Coll., London 


The polykays of the natural numbers—Jn English 


Biometrika (1960) 47, 53-60 (7 references, 1 table) 


The single-suffix polykays of the finite population Fy, 
of numbers 1, 2, ..., N are obtained explicitly in terms 
of Bernoulli numbers and also the polynomials of 
argument (N+1). This was done by using Aty’s 
formulae, for which a general proof has been given 
by Barton & David [Biometrika (1961) 48, 190-191], 
for the moments of the sample mean of n elements 
taken from a finite population in conjunction with the 
moment generating function of that mean when the 
population is also finite (Wilcoxon’s statistic). This is 
achieved by using a limiting process equivalent to 
“ equating coefficients ’ in a polynomial identity. 

The unitary polykays of Fy (with suffix 1”) are shown 
directly to be simple multiples of Nérlund’s generalised 
Bernoulli numbers of order N+1. A series of relations 
giving polykays of suffixes 2°1", 2°, 32°1", 32", 42” in 
terms of the unitary polykays for general finite popula- 
tions are derived as also are recurrence relations for 
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(Univ. Minnesota, Minneapolis and Ohio State Univ.) 


and Univ. Cal., Berkeley) 


polykays: for example between those of suffices 3’, 
53’~', 473""2 and 3'2 and between those of suffices 
rs and rl’. These are used to build up tables of the 
polykays of Fy to weight eight. 

The use of these tables is exemplified by obtaining 
the first four central moments of Mood’s rank-dispersion 
statistic; that is a polynomial in k, and k, in samples 
from Fy. These are compared with the same moments 
derived by the use of Aty’s formulae, noting that Mood’s 
statistic may also be treated as a sample mean from a 
finite population of squares of natural numbers. Some 
polykays are given for this latter population. 


(D. E. Barton) 


1.0 (10.1) 


On bounded infinitely divisible random variables—Jn English 


Sankhya (1960) 22, 253-260 (2 references) 


Let X be an infinitely divisible random variable. Let 
M(u) and N(u) denote the non-decreasing functions 
determining the representation of the logarithm of the 
characteristic function of X. In this paper the author 
proves that for the existence of a constant A such that 
Pr [X¥>A] = 0, that is to say that X is bounded from 
above, it is necessary and sufficient that 


(i) N(uw) = 0 for u>0, 
(ii):a" = 0, 


(iii) lim | M (wu) du =0, when the region of integra- 
e70 
tion for uis —l<u<-~—e. 

A similar result is stated for the case where YX is 
bounded from below: This generalises a well-known 
result that a bounded infinitely divisible random variable 
is necessarily a constant. The results are then applied 
to give information regarding the monotonicity of the 
sample functions of certain stochastic processes. 


(R. Ranga Rao) 
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BIRNBAUM, Z. W., ESARY, J. D. & SAUNDERS, S. C. (Univ. Washington and 
Boeing Scientific Res. Lab., Seattle) 


1.1 (10.8) 


Multi-component systems and structures and their reliability—In English 
Technometrics (1961) 1, 55-77 (4 references, 1 table, 2 figures) 


The first part of this paper presents formal properties 
of general classes of multi-component systems which 
contain two terminal networks and most other kinds of 
systems considered previously as special cases. The 
second part applies the properties to the calculation 
of reliability where reliability is defined as the probability 
that a structure will perform the task for which it was 
designed. This definition implies that a structure can 
only perform or fail. 

Consider a structure of n components: the state of all 
components is described by a vector x = (Xj, ..., X,) 
where the random variable XY; = 1 if the component 
performs and X; = Oifit fails. Assume Pr (Y; = 1) = p 
and Pr (XY; = 0) = 1—p fori = 1, 2,...,n. For certain 
of the 2” possible vectors the structure will perform 
and for the rest it will fail. Let ¢ (x) = 1 if x is a vector 
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for which the structure performs and ¢ (x) = 0 if x is 
a vector for which it fails. (x) is a random variable 
called the structure function or, in short, structure. 
Reliability is then defined as Pr {# (x) = 1} = & [¢ (x)] 
=h(p). Under mild assumptions which are met by 
most practical situations, h(p) is an S-shaped function. 
Thus, there exists a critical value, py, of p such that 
for p>po the reliability of the system is greater than 
the reliability of a single component and for p<pzp it 
is smaller. It is shown that by using replicas of the 
system instead of single components, a system is obtained 
with reliability arbitrarily close to one if p>po; but 
with reliability arbitrarily close to zero if p<po. 


(A. P. Berens) 


1.6 (10.1) 


The strong law of large numbers for a class of Markoff chains 
Ann. Math. Statist. (1960) 31, 801-803 (2 references) 


In this paper the author proves the following result: 
let (Q, BZ) be a measurable space and let Xo, X,, ... be 
a Markoff chain with a compact Hausdorf state space 
Q and transition probabilities Pr (A | x) with the usual 
assumption of measurability on Q for each Ac Z. Let 
@ be the class of continuous real-valued functions on 
Q and define an operator T on @ into @ by T¢ (x) 
= & (¢(X,)| Xo = x). If the Markoff chain satisfies (1) 
and (ii), where (i) states there exists a unique invariant 


measure a such that 2(A) = [es (A | x) 2 (x) for all 


Ae & and (ii) states if fe Y then Tfe @, it follows that 
for any ¢e @ and x EQ we have 


» O(X,)/N> 6,9 (Xi) 
n= 1 
almost surely P,, as Noo where P,, is the probability 


measure for (X,, X2, ...) given Xo = x. 
The proof is accomplished through two lemmas, 
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N 

firstly, that for any ¢ € @, )) T"$/N converges uniformly 
1 

to &,¢(X,), and secondly, that 


(1/N) Y [é (6 (X) |X.) -€ OX) | Xv] 


tends to zero almost surely for each k. The first, well 
known in linear space theory, is proved for completeness: 
the second follows from the Markoff property and a 
text-book stability theorem for a sequence of random 
variables centred at conditional expectations given the 
predecessors. 

This result established that for learning models 


N 
YX /N>6,X, 
1 


almost surely where x is the unique invariant initial 
distribution. 


(S. C. Saunders) 
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EHRENFEUCHT, A. & FISZ, M. (Math. Inst., Polish Acad. Sci., Warsaw) 


1.5 (1.6) 


A necessary and sufficient condition for the validity of the weak law of large 


numbers—ZJn English 
Bull. Acad. Polon. Sci. III (1960) 8, 583-585 


The following theorem is proved: Let {x;}; i = 1, 2, ..., 
be a sequence of independent and equally distributed 
random variables, and let ¢(t) denote their common 
characteristic function. The relation 


lim Pr [ | (x, +...+x,)/n—a | >e] =0 


n> oo 


holds for any e>0 if and only if @(f) is derivable at 
f= .0 and o’ (0) = ai. 


(L. Kubik) 
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(Math. Inst., Hungarian Academy of Sciences, Budapest) 


1.3 (1.5) 


On the strength of connectedness of a random graph—ZJn English 
Acta Math. Acad. Sci. Hung. (1961) 12, 261-267 (5 references) 


The following problem is posed by the authors: let 


there be given n labelled points V,, V2, ..., V,. Let us 
choose at random WN edges among the (") possible 


edges connecting these n points so that each of the 
n 
((3)) possible choices chould be equiprobable. Denote 


by T,y the random graph thus obtained. Let d(V,) 
denote the valency of a point V,, the number of edges 
starting from V,, and, for a graph G, c(G) = min d(V,). 


1sk<n 
Let c, (G) and c,(G) denote the least numbers k and / 
such that deleting k appropriately chosen vertices, 
together with all edges starting from these vertices or 
1 edges from G the resulting graph is not connected. 
Let v,(G) denote the number of vertices of G which 
have the valency r; r = 0, 1, 2, .... 
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The authors prove that if 
N (n) = 4n log n+4rn log log n+an+o(n) 


where « is a real constant and r a non-negative integer 
then (Theorem 2) 


lim Pr [e, (Un, nm) = r]= lim Pr [e. (Un, nim) = 1] 


Nak F 60 n> 
=" lim iPr fic (1, ny) = = exp Caen a1) 
n> +0 

and the distribution of v, (T,, yin) tends to the Poisson 
distribution with expectation e~ ?%/r! (Theorem 3). 

Theorem 2 is a generalisation of a previous result, 
see Erdés & Rényi, ‘“On random graphs. I.” [Publ. 
Math. Debrecen (1959) 6, 290-297: abstracted in this 
journal No. 2/238, 1.3]. 


(K. Sarkadi) 
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GHELLINCK, G. de (Operational Research Centre, Brussels) 


Sequential decision problems—ZJn French 
Cahiers Cent. Etudes Recherche Opérat. (1960) 2 


This paper is concerned with multi-stage decision 
problems. In such a problem a decision concerning a 
system is to be made from a set of alternatives at 
Successive intervals of time. The outcome of any 
decision is to define a distribution of transition prob- 
abilities as well as an expected reward during the period 
running until the next decision. 

The problem is to find a decision policy which in the 
long run yields the maximum expected reward per 
period. The number of states of the system as well 
as the number of aiternatives available in any state 
are both assumed to be finite. There are three different 
methods to solve this kind of problem; the purpose 
of the paper is to analyse their mutual relations. 

The dynamic programming approach as presented by 
Bellman, see ‘“‘ Functional approximations and dynamic 
programming” [Math. Tab., Wash. (1959) 13, 247-251: 
abstracted in this journal No. 1/523, 11.5], is the more 
general case. This consists in assuming that the process 
ends at a given time. Using recurrently the “ principle 
of optimality” one finds for each state of the system 
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GHELLINCK, G. de (Operatiorial Research Centre, Brussels) 


Sequential decision problems—ZIn French 
Cahiers Cent. Etudes Recherche Opérat. (1960) 2 


paper analyses their common features and states several 
propositions. As an illustrative example a replacement 
problem is successively handled by the three methods. 

Whatever was said about the processes that make 
changes at discrete intervals of time, may be adapted to 
continuous-time processes. The critical feature to be 
maintained, is the finiteness of both the number of 
states and the number of available alternatives. 


(G. de Ghellinck) 
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1.8 (0.8) 


and for each period an optimal decision as well as the 
corresponding maximum expected reward: such a 
procedure converges to a permanent optimal policy. 

In contrast with this “ value-iteration” approach, 
Howard recently presented a “ policy-iteration ” method 
that uses the finite features of the model more efficiently 
and is focussed on permanent policies only. Starting 
from an arbitrary decision policy, the algorithm first 
determines relative “values” associated with every 
state of the system. These values then form a criterion 
that leads toward an improved policy, if any. After a 
finite number of such steps, one reaches the optimal 
policy. 

As a third method, one may state the problem in 
terms of linear programming. The non-negative 
variables are probabilities imposed upon all the available 
alternatives as a mixed strategy. These probabilities 
are constrained by linear equalities in order to be con- 
sistent with steady-state decisions. The expected reward 
per period is the linear objective function. 

After a brief summary of these three methods, the 


continued 


1.8 (0.8) 
continued 


ie PW oe Palate Se eg Se Iiubivoat deny tl Sota 
PUD | hy tary) eee eo hs ae ae 
a Leniityy 2 eetory. 6 Atay qian 
: a » Silt rer = es iam ise Gwe sak 


ha Ine he oN ee Peo iy 
fas vie oa tinecees oe “he Sangeet) A Saas 
aniiagty oe Wee ceil they adhe 
SA ge ite ' aah: ) See 
gf Rise endl CO SNL. i ‘a i 
i = 406.0 Sey besaiepy ty Pare Par gid 


ha, 


i vib 
onal 


Witte? Wray, 
I, OK 
ret 


pee Boney 5.0" 
ite oe a) ae ae" 
Pebiesi Weer ea 1) TR 


2 IPTC Ws A aes, 
alli 


ree 2. ! 
an 


awe 


“sy teal (Wighis Tape lt te he sities, Cuseea Siw ybh: 


Preyer’ 


i vires He ei ba Se eee Soa 
; : Ee 
yi veh apt Avaie eit: aah fs nal apaeg AZ 
Fela trimiysihe we “yeas travel 2 ee / 
Pe WET Fb, ta elutes cee ee eg) ete ah gees be 
Prehary culls STW jy el arpt in gle cmend Sie tp 
ur ot Ol VS OF Stee i ae ay nf Cao 


56 loa 


a) 


Bic: ne we 


eee 


i oh ETE, 4 ae hy 
ohh : 
‘sh a “i gui 


i Bg Teal 
7 WINE "nen peas 
hao naan A ahaa 
‘el Pits ee 


0:1 palate yillaliags ae er, 
one ee a” ane 
See oe sir hala a 
ta. Aha le ONE a arg 
ta ee TL aa Widhul'g asikits yas ort 
Satine ti apnea ie eed al 
“7. 2 a ee Denk ia 
£7 ty +g (in oe 
2 eh ne, og Gia relies A ma 
CCDS TD AT ie Si 
Satin by (ae ee x nai eel 
sowed mit rely dea aeenOe Liles eae 
vapeain] "-all Cone me ers ty . 
cette: 24 ee he ay a he 


AiG: sta Wie < 


i AP aR ie 


\ 


15h SA led ‘tet eae alii an ay 


ecw ehgatg 
i ibe dai 


ees 


& GAN LER 


is. Ui a | 
eee o Vie ae 

“nde deel pelt i Pi 
CO) Sere ae eee lees, Ae 
| rae en > haha 
iO) Cisahan aes 
si site 2 ARS 
| - lotta vie a 


GLADYSZ, S. (Math. Inst., Polish Academy of Sciences, Wroctaw) 


An ergodic paradox—In German 
Coll. Mat. (1960) 7, 245-249 


Let X be a space, € a class of its subsets E, and m (E) 
a probability measure defined on €. The author 
considers a measurable measure-preserving transforma- 
tion Tx and a “ distance”? p(x, y) defined on X and 
such that p (x, Tx) is a measurable function with finite 
expectation. In a sequence of transforms x, Tx, T’x, ... 
of a given point x we can distinguish those points T"x 
which are in a given set E. If n, and n, are two con- 
secutive integers for which it is so, we cali the sum 


p (Ex ge hae sa 
+p CTS te T"*7x)+...4 9 (hea Xe, T"x) 


the length of way made by x between these two con- 
secutive sojourns in E. This gives rise to a sequence 
ry (X), rz (x), -.., of lengths-of-ways made by x between 
consecutive sojourns in £. Under usual conditions 
on T it can be shown with the aid of Birkhoff’s ergodic 
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GYIRES, B. (University of Debrecen, Hungary) 


1.5 (10.0) 


theorem that for a given E with m(E)>0 the mean of 
r’s, that is 


lim + [r, (x) +... +], 
k> 0 k 


exists with probability one and does not depend on x. 

However, as shown by the author, it cannot be so 
for the mean of every sth length. Namely; an example 
is given where for s = 2 the mean 


Mee [r2 (x) +14 (X)H... +14 (x) /k 


exists with probability one, but takes on essentially two 
different values. 


(S. Zubrzycki) 


1.5 (0.6) 


A generalization of the local form of the central limit distribution 


theorem—Jn Hungarian 


Magy. Tud. Akad. III. Oszt. K6zl. (1960) 10, 469-479 (3 references) 


The author generalises the local form of the limit 
distribution theorem concerning a sequence of equally 
distributed independent random variables. He starts 
by observing that the local limit distribution theorem is 
equivalent to the asymptotical investigation of the powers 
of a stochastic matrix, and is infinite in four directions. 
The author investigates a stochastic matrix infinite in 
four directions, which is generalised Toeplitzian in the 
sense that the elements of the matrix are replaced by 
quadratic matrices of finite order. It is shown that, if 
the elements of the mth power of such a matrix are 
denoted by A\”; k =0, +1, +2, ..., then, under certain 
conditions and with suitably chosen constants B,, the 
relation 
; TNO ze 
B, A?” — —— e2 P-0 
APs 

holds for noo, where P is a stochastic projector matrix 
and the convergence is uniform in k. 
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A probabilistic interpretation of this result is also 
given. The results of the paper are obtained by the use 
of a generalisation of the method of characteristic 
functions. 


(T. Balogh) 
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KARLIN, S. & PROSCHAN, F. (Stanford University) 1.0 (2.0) 
Pélya type distributions of convolutions—Jn English 
Ann. Math. Statist. (1960) 31, 721-736 (16 references) 


A function f(x, y) of two real variables is said to be 
totally positive of order k (TP,) if for all 1<m<k, all 
= Noe ky ANG Y,<,<...<y,, the determinant 
|f (xi y;) | 0; 

If a TP, function f(x, y) is a probability density in x 
with respect to a o-finite measure (x) for each fixed y, 
then f (x, y) is said to be of Pélya type of order k (P7;,). 
If a PT, frequency f(x, y) can be written in the form 
Jf (x—y) then f(u) is called a Pélya frequency density of 
order k (PF,_). 

The author studies these classes of functions and 
obtains interesting results concerning the n-fold con- 
volutions of PF,’s. He also studies the variation- 
diminishing and smoothing properties of 7P, functions 
as well as the properties of compound distributions 
composed from PT, densities. 


(E. Lukacs) 
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KRAFT, C. A., PRATT, J. W. & SEIDENBERG, A. (Michigan State Univ., 1.1 (1.0) 
Univ. Chicago & Harvard Univ., and Univ. California, Berkeley) 


Intuitive probability on finite sets—Jn English 
Amn. Math. Statist. (1959) 30, 408-419 (4 references) 


This paper presents further results on the axiomatic By artificially enlarging the set of elementary events, 
development of intuitive probability on finite sets based on a construction in which each event is decom- 
introduced by de Finetti [see L. J. Savage, Foundations posed into two equally likely and exclusive events, and 
of Statistics (1954) New York: Wiley]. If one presumes adding intuitively justifiable axioms, they are able to 
that, for every pair of compound events, one can judge prove in this new system that any ordering is compatible 
which is more likely or whether they are equally likely, with a measure. 
then it is possible to assign a numerical probability to 
each elementary event in a way consistent with these (W. J. Hall) 
judgments; that is, does there exist a strictly, or almost 
strictly, agreeing measure? De Finetti conjectured 
Ves: 

The authors show that the answer is “‘ no ” by present- 
ing counter examples, in which there are five elementary 
events. They develop a necessary and sufficient condi- 
tion that an ordering of the compound events which 
satisfies de Finetti’s axioms does arise from a measure. 


2/624 


Mao ee A 8 ee ‘a Borie i 
boiih myirah. Pan iG eee pe: an 


= ee ole ll 


é Hag) Geka CERO 
LOM he 
_ nite (othe al Ree a. 
Bios mera agen “a ¢ a 
2 Li ee: PhS ‘ae cin 
ve 1) Lit? tk aa He 
‘teh. aly Eliana anaes oe 


=v volitional! co inte 
ie ET 


. 8 .) 


S°e dag. =am slave ae fairy 
2 sige f ed 7 
' ) Ae oe HE i€ Beas _ aa 


Ai wae aie allel ni 


” sth eee i” (Brut ms | ie. 
a or 


ee "oe ¢° ares Kilt ‘ay * dvi 
ly : < 7 ae ef nt. Ate ‘oe nai 

ie ating “ 5 * ten ion a7 Yea\ 

seen om 5 eee pee et ee = a 

a 

se ew 

Babee hi: Bh i. ; 


ONICESCU, O. (University “C. I. Parhon ”, Bucharest) 


Sum functions on a Boolean algebra—In French 


1.1 (1.0) 


Bull. Math. Soc. Sci. Math. Phys. R.P.R. (1959) 3, 77-91 


A new foundation of probability theory is given, replac- 
ing random variables by sum functions. Thus, let 
(Q, P) be a probabilised Boolean o-algebra. By defini- 
tion, a sum function on Q is a real valued completely 
additive function F on Q; a P-sum function on Q is a 
sum function F absolutely continuous with respect to 
P. A P-sum function is simple if it is of the form 


n 
F(@)= ¥ a;P(6;A@) where a; 1SjSn, are given 
j=1 
real numbers and 0;, 1<j<n, is a finite partition of the 
unit element of Q. 

Operations with simple P-sum functions as additivity, 
multiplication and convolution are defined ; moreover, 


a theorem concerning unicity of the representation of 


simple P-sum functions is proved. It is also shown that 
every P-sum function F of the form G—H, G(w)=0, 
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PALASTI, Iona. (Math. Inst., Hungarian Academy of Sciences, Budapest) 


On some random space filling problems—ZIn English 


H (@)20, wm €Q, is the uniform limit of a sequence of 
simple P-sum functions. Further, a canonical represent- 
ation of simple P-sum functions is provided and its 
properties are investigated. 

Finally, two theorems concerning the local and 
integral representation of P-sum functions are given 
and operations with P-sum functions in the wide sense 
are defined. 


(R. Theodorescu) 


1.4 (11.7) 


Publ. Math. Inst. Hung. Acad. Sci. (1960) 5, 353-360 (3 references, 1 table, 1 figure) 


The following two-dimensional model of random space 
filling is considered by the author in this paper. Domains 
which are congruent and parallel with a given domain D 
of unit area, are placed at random into a rectangle with 
sides x>1, y>1. Ifa domain intersects any of previously 
placed ones it is discarded. The random space filling is 
finished when there will be no further possibility of placing 
a domain without intersection. Let M (x, y) denote the 
mean value of the number of domains placed by applying 
this model. 
The author proves that the limit 


M (x, y) _ 
xy 


a (D) 


lim 
ae 
exists if the following hypothesis is true: there exists a 
constant A>0O depending on D, such that for any 
x,>0, x,>0 and y>0 we have 
| M (x,+X2, y)—M (x1, y)—M (x2, y) | SAy 
respectively for any x>0, y,;>0, y,>0 


| M (x, y1+32)—M (x, ¥1)—M (x, y2) | SAX. 
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The proof is based on a lemma which is a generalisation 
of a theorem of Hyers, related to functional inequalities. 

The author deals with the special case when D is a 
rectangle. An analogous problem in the one-dimensional 
case was solved by Rényi [Publ. Math. Inst. Hung. Acad. 
Sci. (1958) 3, 109-127: abstracted in this journal 
No. 1/18, 1.4], and the limiting relation 

lim M (x)/x =C~0-748 

was found. The author’s conjecture «(D) = C? is in 
very good agreement with the results of some Monte 
Carlo experiments carried out. 

The corresponding n-dimensional (n= 3) random space 
filling problems are said to be treated in the same way 
as the two-dimensional one. 


(G. Bank6vi) 
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SHAPIRO, J. M. (Ohio State University, Columbus, Ohio) 
Sums of small powers of independent random variables—In English 
Ann. Math. Statist. (1960) 31, 222-224 (3 references) 


Given a double sequence of infinitesimal independent 
random variables, {x,,}, the author studies the distribu- 
tion function FY (x) of 


Si = » | Xnk \"—B, (r); 


where 0<r<1 and B, (r) is a suitably chosen constant. 
Necessary and sufficient conditions are given for Fy (x) 
to converge (n> 00) to a distribution function [F. (x)]. 
If F, (x) converges (r>0+) to H (x), H'(x) is the dis- 
tribution function of the sum of two independent 
random variables, one Poisson and the other Gaussian. 

In a previous paper [Ann. Math. Statist. (1958) 29, 
515-522] he considered the case of r>1. 


(R. L. Anderson) 
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TROTTER, H. F. (Kingston, Canada) 
An elementary proof of the central limit theorem—ZJ/n English 
Arch. Math., Karlsruhe (1959) 10, 226-234 (2 references) 


The proof given by the author is similar to that given 
by Lindeberg [Math. Zeit. (1922) 15, 211-225] and is 
elementary in the sense that it avoids the use of char- 
acteristic functions. The similarity ends when the 
author introduces the new idea of the characterisation 
of the distribution of a random variable X by a con- 
tracting operator 7, in the space of continuous functions 
(Txf)@) = & [f(X+9)]. 

This operator has the property Ty,y = Ty7y for the 
independent random variables, X¥ and Y. For con- 
vergence in distribution it is sufficient to have 


(Tx, —Ty) f (0)>0 


for twice differentiable functions /. 
This criterion is applied by means of Taylor’s expan- 
sion and the result is generalised to several dimensions. 


(D. Morgenstern) 


2/628 


1.5 (2.9) 


1.5 (-.-) 


. 


OT Oe eee eee 
é 


oe an tT bes ae 
ae ha ail ci foe fi 


ae ee 
wie ‘Soi Rs . 


a 
watiees ign Adana on 
oP Set iy rk 
He ee 

a alee VE ae 

scape (ent AVENE Ls na an “e's 

1 ale? satan ot yt ean 
or OY VO) 9 rg ilo Ee 

(i an oa 


* roaryag iyi : pie i " a eT 
~ u co on He iy 


: SD bings 
Wer Palle} of aust jaye aT es pee: 


(ont © ss a ii 


«) ea aay die. anil ah Se 
| > oe, See oe a ee aan 1 ie 
oe as REN 
< | “A Ct lat | yeaa 

ag; lit yen tga Og glee =e 


ae 


Vd 


ULLRICH, M. & URBANIK, K. (Inst. Radio Eng. and Electronics, Czech. Acad. Sci., 


1.5 (1.0) 


Prague and Math. Inst., Polish Acad. Sci., Wroclaw) 
A limit theorem for random variables in compact topological groups—In English 


Coll. Mat. (1960) 7, 191-198 


Let X,, X2, ..., be a sequence of independent random 
variables with values in a compact topological group G, 
which are symmetrical. The authors consider probability 
distributions of random variables Y, = X,-X,-...-X,, 
where the product is taken in the sense of group multi- 
plication. The question considered is under what 
conditions imposed on the probability distributions, the 
H; of the X;’s, the probability distributions y, of the 
Y,/s are weakly convergent to the uniform probability 
distribution over G. The authors prove that the 
existence of a positive probability distribution 2 over 
G, that is such that 4(V)>0 for every non-empty open 
subset V of G, and of a sequence a;, a>, ... of numbers 
with OSa,<1 such that for every Borel subset E of G 


we have 
@) uw, (E)24,A (E) 


(ii) a,+a,+... = 0 


and 
is a sufficient condition for this. Moreover the authors 
show that the condition (ii) is essential in the sense that 
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WOLFF, P. M. de. 


(Philips Labs., Irvington-on-Hudson, N.Y.) 


if only the series in it is convergent, then there is a 
sequence of X,’s so that the probability distributions of 
Y,/s do not converge to the uniform probability dis- 
tribution over G. 


(S. Zubrzycki) 


1.2 (1.4) 


Particle statistics in X-ray diffractometry. I. General Theory—IJn English 
Appl. Sci. Res., Hague (1958) 7, 102-112 (1 reference, 3 figures) 


In X-ray diffraction by powders, a statistical error in 
the measured intensity results from the finite number 
of crystals contributing to it. A recent device for 
reducing this error, the rotating sample holder or speci- 
men spinner, is not covered by existing formulae for 
stationary specimens. The present paper establishes 
an expression for this new situation, which has also been 
investigated experimentally in the second part of this 
paper [J. Appl. Phys, (1959) 30, 63-69]. 

The contribution of each crystallite can be written 
as a product of functions of mutually independent 
variables; its volume, and its position and orientation 
coordinates. The resulting root mean-square error is 
expressed in the relative standard deviations of those 
functions. The latter are calculated using somewhat 
simplified physical assumptions. This is done first for 
the stationary specimen, where it is found that the form 
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of the distribution of intensity contributions has surpris- 
ingly little effect on the statistical error. The same 
conclusion is reached for rotating specimens, where a 
very rough estimate of the number of contributing 
crystals leads to an expression differing only in the 
numerical factor, 4 against 6-5, from the final expression 
which is obtained by detailed analysis of the rather 
intricate geometry of this case. 


(P. M. de Wolff) 
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EPSTEIN, B. (Wayne State University and Stanford University) 
Elements of the theory of extreme values—In English 


Technometrics (1960) 2, 27-41 (8 references) 


The purpose of this paper is to describe briefly the basic 
elements of the theory of extreme values. Frequent 
examples are given to illustrate both the derivations and 
applications of the theory. 

Throughout the article, the author deals with samples 
of size n, drawn independently and at random from a 
parent population with density function f(x) or dis- 
tribution function F(x). Of particular interest here is 
the minimum value, y,, and the maximum value, zo 
from such a sample. 

The density function and distribution function for ve 
are derived in terms of F(x). Examples are given for 
parent uniform and exponential populations. 

Using a technique of Cramér, the author asserts that 
the sequence of random variables, y,, converges in 
distribution to a random variable y, whose distribution 
is more easily evaluated. Examples are given where 
the parent populations are selected to illustrate Gumbel’s 
types 1, 2 and 3 asymptotic distributions of smallest 
values. The conditions for each type of asymptotic 
distribution is given. 
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FREIMER, M. & GOLD, B. (Lincoln Lab., Massachusetts Inst. Tech.) 


2.6 (1.5) 


The equivalent results for the largest value, z,, are 
also derived, with corresponding examples and their 
relationship to Gumbel’s classification of asymptotic 
distributions of largest values. 

For the case when f(x) is differentiable, the author 
shows several methods by which the modes of the 
densities of y, and z, may be found. These methods 
are carried out for some of the preceding examples. 

The article concludes with remarks on the applications 
of this theory to a variety of fields. 


(L. A. Foster) 


2.1 (2.9) 


Note on the distribution of locally maximal elements in a random sample—Jn English 


Ann. Math. Statist. (1960) 31, 518 (3 references) 


This note presents a simpler second factorial moment 
of the distribution referred to in the title than presented 
by Glasgow [Anmn. Math. Statist. (1959) 30, 586-590; 
abstracted in this journal No. 1/190, 2.1]. 


(R. L. Anderson) 


eA vel 
‘ hi 


a 


GOLOVIN, N. E. (Natl. Aeronautics & Space Admin., Washington, D. C.) 


2.5 (1.0) 


An approach to a reliability prediction programme—In English 
Annual Cony. Trans., Amer. Soc. Qual. Contr. (1960) 173-182 (4 references, 5 figures) 


Before there is useful operating experience for a system, 
reliability prediction is valuable in forecasting operational 
feasibility and costs versus reliability. A failure is 
_ defined as occurring at the time any relevant physical 
characteristic of the system attains a value outside of a 
specified range. A full reliability analysis would dis- 
tinguish between catastrophic and degradation failures. 
The procedure recommended here is to start with 
system failures and work back through sub-systems and 
components to the data on parts failures. A given set 
of mathematical relationships can be taken which allow 
calculation of system outputs from sub-system outputs. 
In this way, system failures can be traced to sub-system 
failures. The effects of variation in sub-system outputs 
on the system output can be studied. Grouping sub- 
systems into independent groups allows the use of the 
product rule for determining probabilities. Components 
of a sub-system and parts of a component can be treated 
similarly. Any prediction of reliability must be done 
for a specified set of environmental conditions. 

The reliability, that is the probability that a part 
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JOHN, S. (Indian Statistical Institute, Calcutta) 


survives to time f, can be expressed in general as 


R, (t) = exp & { A, (t) a 


where J, (t) is the failure rate of the pth part per unit 
of time. Frequently, 4,(¢) increases rapidly at first 
during an “infant mortality’ period, decreases to a 
fairly constant level, chance failures, during a ‘‘ normal 
operating’’ period, and finally increases during a 
““wear-out’”’ period. For a component with n in- 
dependent parts 


R, (t) exp] - | ys A; (t) a| 


If there are k identical parts in a redundant arrange- 
ment, the reliability of the group is 


R,(t) =1-[1-R, @]*. 


The validity of a reliability prediction depends on the 
ability to identify the independent groups of a system 
and on the information available for the failure rates 
of the groups, J, (¢). 


(D. W. Gaylor) 


2.9 (0.3) 


On the evaluation of the probability integral of a multivariate normal 


distribution—ZJn English 
Sankhya (1959) 21, 367-370 (7 references) 


It is shown by the author in this paper that the probability 
integral of a multivariate normal distribution of order k 
can be expressed as the sum of k terms each of which 
is the product of a univariate normal density and multi- 
variate normal integral of order k—1. To illustrate 
the application of the formula an example is given. 


(T. V. Hanumantha Rao) 


KOTLARSKI, I. (Warsaw) 


On random variables whose quotient follows the Cauchy law—Jn English 


Coll. Mat. (1960) 7, 277-284 


The author considers a class & of random variables € 
with the property that the quotient of two independent 
random variables having the same distribution as € has 
Cauchy’s distribution. Owing to the connections 
between Mellin transforms h(z) = E| é |*, z complex, 
and characteristic functions ¢(t) = &e'*, where ¢ is 
real of random variables €, he characterises the class X 
in terms of Mellin transforms. 


(S. Zubrzycki) 
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PATIL, G. P. (Indian Statist. Inst. and Univ. Michigan, Ann Arbor) 
On the evaluation of the negative binomial distribution with 


examples—ZIn English 
Technometrics (1960) 2, 501-506 (11 references) 


In this paper the author proposes an evaluation for 


Y(r, P; k) = ye. y (x; DP, k) 
of the negative binomial density function 


k+x-1 


kA~x ile 
7 ) phat, p+4 


y (x, pK) =( 
ek =) Voe2, 2. then 
il Bea Rie 
Wasa ke a gs te dL, 
aay ra 


u+v = 1; namely the incomplete beta function. 

A short proof of this is given by the author and is 
as follows: 

Y(r, p, k) = Pr [at most k+r independent trials are 
required to achieve k successes when p is 
the probability of any single success] 

= Pr [at least k successes occur in k+r 
trials | 
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2.9 (0.4) 


2.5 (0.4) 


kt+r , 
= eee which by a well- 


x=k x 


known theorem connecting the binomial 
and beta distribution 


if Dp 
= (rt i= dy: 
p( ) mel, 


A numerical example is given. 


(E. H. Lehman, Jr.) 


MENDENHALL, W. & LEHMAN, E. H., Jr. (North Carolina State College) 
An approximation to the negative moments of the positive binomial useful in life 


testing—ZIn English 


2.4 (4.3) 


Technometrics (1960) 2, 227-242 (10 references, 4 tables, 1 figure) 


The Weibull density function 

F(t; m, a) = (m/a) t™~* exp (—2"*) 
is useful in life testing. Assuming m known, n items are 
placed on test and the test is stopped when time T has 
elapsed. It is assumed that T is determined in such a 


way that the probability of zero failures is negligible. 
Then the maximum likelihood estimator for « is 


a 


= 1 


y+ L(nfr)—1] 7", 


where ¢; is the time until the ith ordered failure and r 
is the number of failures. 
The expectation of &@ is shown to be 


a+T™ [né (1/r)—(q/p)—1] 
where gq = exp(—7J”/a) is the probability of survival 
until time 7, p = 1—g, and r is distributed as a positive 
binomial: & therefore is considerably biased. 
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testing—ZIn English 


(North Carolina State College) 
An approximation to the negative moments of the positive binomial useful in life 


The variance of @ then is 


V (@) =T2™ ly (+5 st) 


eet 


+eer(2) +2n cov (; yx ‘)| 
AN r r 


The covariance term is shown to be zero and then, after 
some manipulation, 


ra rofon()[3)- A) C) 


Hence, the first two negative moments of r, which is 
distributed as a positive binomial, are required: to 
approximate these the ‘‘ best”? beta function is used. 
It is noted that if y has the beta distribution with para- 
meters a = (n—1) p and b = (n—1)gq, then 


e(*) a n—2 (4) _ (@=2)(n—3) 
y ai Ve (a—1)(a—2) 
continued 


2.4 (4.3) 
continued 


Technometrics (1960) 2, 227-242 (10 references, 4 tables, 1 figure) 


Now r can be approximated by ny, for 


é (9) =np = —— = 6 (ny) 


+b 
and ; 
2 2.2 n°a(a+1) = & (ny? 
= = = y°). 
Ug ae ig (a+b)(a+b+1) 


1 1 ‘ 
Therefore & a and & (3) can be approximated 


1 ie (n—2) 
(5) n(a—1) 
iio )- (n—2)(n—3) 

n?y? n?(a—1)(a—2) 


respectively by 


and 
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Then V (1/r) is approximated by 
(n—2) (n—a-1) 
n ta—)? (a2). 
& (1/r) and & (1/r?) are tabulated for n = 5, 10, 20, 30, 


40, 100 and p = 0-1, 0-3, 0-5, 0-7, 0-9. 


The approximation is good to two significant digits 
for np2=5. 


(E. H. Lehman, Jr.) 
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COHN, R., MOSTELLER, F., PRATT, J. W. & TATSUOKA, M. 
(Rutgers Univ., Harvard Univ., Harvard Univ. and Univ. of Hawaii) 


3.8 (5.6) 


Maximising the probability that Adjacent. order statistics of samples from several 


populations form overlapping intervals—Jn English 


Ann. Math. Statist. (1960) 31, 1095-1104 (6 references) 


The problem treated here arose in a sociological study, 
but is stated abstractly as follows, ‘‘ Provide a test of 
the hypothesis that several samples of the same size are 
randomly drawn from possibly different populations, 
against the alternative that the samples are not in- 
dependently and randomly drawn from distributions ”’. 
The authors consider the interval, in each of the k 
samples of size n, between the rth and (r+1)st order- 
statistics. It is determined whether these k intervals 
have a non-empty intersection; or, equivalently, whether 
the least of all the rth order-statistics is greater than or 
equal to the greatest of all the (r+1)st order statistics. 
In this case, one looks to see with what probability it 
would occur had the samples been randomly and in- 
dependently drawn from the same population. If this 
probability is sufficiently small, one rejects the hypothesis 
that the samples were randomly and independently 
drawn from, possibly different, populations, since it is 
shown that the probability of a non-empty intersection, 
when the samples are drawn from k univariate continuous 
cumulative distribution functions on the same real line, 
2/639 


is greatest when the populations are identical. 

After treating the case r = 1 in Sections three and 
four, the authors prove two general theorems in Section 
five. If x,; denotes the rth observation, in descending 
order of magnitude in the ith sample, the theorems are 
firstly, if samples of size n are drawn from each of k 
continuous distribution functions y;(x); i= 1, 2, ..., k, 
the probability that ate egrets Nica LOEra ao 


a, ae —1) is Hater if and only if y; = y2 
= =e. bee secondly under the same conditions, 
the probability that 

min 2b pens Xp+1,¢ and min x,;>max X544,; 


U t 

fora eter i (= ie 2,....N—1)ands(=r7,r+1,...,n—1), 
is) maximised if and only if y, = y, =... = yp. In 
each theorem a combinatorial formula for the maximum 
probability is given. The results extend to three or more 
corresponding intervals. The authors describe obstacles 
in the path of extending the results from the univariate 
to the multivariate case: and remark on some other 
possible extensions. , 

(R. Hooke) 


3.8 (2.7) 


GUPTA, S. S. (Bell Telephone Laboratories, Allentown, Pennsylvania) 
Order statistics from the gamma distribution—IJn English 
Technometrics (1960) 2, 243-262 (7 references, 3 tables) 


This paper is concerned with order statistics from the 
gamma or Pearson Type III distribution. The first 
four moments about the origin and the mean are 
calculated and tabled for the kth order statistic from a 
sample of size N. An expression is derived which would 
enable one to evaluate the covariance between the mth 
and nth order statistics. 

The distribution function of the kth order statistic 
is obtained in terms of the incomplete beta function 
and also in terms of the incomplete gamma function 
and the cumulative Poisson distribution. An iterative 
procedure is presented for finding the mode of the 
kth order statistic and a table is presented evaluating 
this mode with k = 1 and k = N for several values of 
the gamma-distribution parameter r. If B, (k, N—k+1) 
denotes the «-percentage point of the beta distribution, 
then Y,, the «-percentage point of the kth order statistic, 
is the solution of 

G. (Y,).= 8, (k, N-kK+1). 
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An iterative procedure was used to solve this equation 
and Y, is tabled for certain values of N, k, r and «. 
If the gamma density function is considered in the 
form 
Se Wem et | 
Gr (x, 0) a: Soe 
é'T(r) 
where 0 is to be estimated, then a method of finding the 
best linear unbiased estimate of 0 is presented in the 
case of censored sampling. 
Applications to life tests, extreme values, reliability 
and maintenance are described and illustrated in some 
cases. 


ses), DSO; 


(A. P. Berens) 


HOGG, R. V. (University of Iowa, Iowa City, U.S.A.) 
Certain uncorrelated statistics—Jn English 


3.61675) 


J. Amer. Statist. Ass. (1960) 55, 265-267 (1 reference) 


The author proves that an odd location statistic and 
an even location-free statistic are uncorrelated when 
sampling from a symmetric distribution. An odd 
location statistic, T, by definition is transformed similarly 
to the sample space under (i) translation or (ii) multi- 
plication by —1, while an even location statistic, S, by 
definition is invariant under (i) or (ii). The theorem of 
the author is stated as: . 

Theorem. Let X;, X, ..., X, be a random sample 
from either a continuous-type or a discrete-type distribu- 
tion that is symmetric about the point c. If the correla- 
tion coefficient of an odd location statistic T(X,, X>, 
..., X,) and an even location-free statistic S(X,, X>, 
..., X,) exists, it is equal to zero. 

The proof may be outlined as follows: the expectation, 
6 [T], of T is c. Thus the covariance, V, of T and S 
is &[(T—c) S]. Employing the transformation x;=y;+c 
and the definitions, it follows that V = & [TS]. By 
the transformation y; = —z; and the definitions it 
follows that V = —& [TS]. It is thus implied that 
V = 0 and hence the correlation is zero. 
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(Wright Air Development Centre, Dayton, Ohio) 


The author notes that the theorem is actually valid 
for every random vector whose distribution is sym- 
metric about (c, c, ..., ¢) in n-dimensional space. 
Restriction to stochastically independent and identically 
distributed random variables and to either continuous- 
type or discrete-type distribution is unnecessary. 


(S. Krane) 


3.8 (2.9) 


Variance of the median of samples from a Cauchy distribution—Jn English 
J. Amer. Statist. Ass. (1960) 55, 322-323 (2 references, 1 table) 


The author obtains the integral expressing the exact 
variance of the median for samples of size 2k+1 from 
the Cauchy distribution with location parameter 0. 
The variance is computed by numerical integration for 
k = 2(1) 15. Corresponding values of the approximate 
variance formula for large samples are computed for 
comparison. Values of the intrinsic accuracy, that is 
the variance of logarithmic derivative of density function, 
are also given and the three functions are displayed in 
tabular form. 

The large-sample approximation is found to be 
poor for small samples. For k = 2, for example, the 
exact variance is given as 1:221253, the large-sample 
approximation as 0:493480 and the reciprocal of the 
given intrinsic accuracy is 0-572895. For k = 15 the 
corresponding values are 0-087944, 0-079594 and 
0-085545. 


(S. Krane) 


2/642 


SARKADI, K. (Hungarian Academy of Sciences, Budapest, Hungary) 


38 (1.3) 


On the median of the distribution of exceedances—Jn English 
Ann. Math. Statist. (1960) 31, 225-226 (4 references) 


If we have two random samples of sizes n, and n, from 
the same continuous distribution, the number of 
exceedances (x) is defined by the author as the number 
of elements of the second sample which surpass at least 
ny—m-+1 elements of the first, for a fixed natural number 
mSn,. The cumulative distribution of x, W (n,, m, Ny; X) 
is given; ifm, = n, the median is m—1. 
In this paper it is shown that 


W (ny, M4, Nz, m,—1)+ W (ng, mp, ny,m,—1) = 1. 


(R. L. Anderson) 
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SIDDIQUI, M. M. (National Bureau of Standards, Boulder, Col., U.S.A.) 


3.8 (4.4) 


Distribution of quantiles in samples from a bivariate population—ZJn English 
J. Res., Nat. Bur. Stand. B (1960) 64, 145-150 (3 references) 


This paper extends the results of former investigators— 
the exact distribution of quantiles in the univariate case, 
the distribution of medians in the multivariate case— 
to the distribution of quantiles in samples from a 
bivariate population. The general approach follows 
that of Cramér and Mood. 


statistic. 

The author considers the joint distribution F(x, y) 
of (X, Y), possessing a probability density function 
f(x, y), and having F, (x) and F, (y) as the respective 
marginal distribution functions. Let « be a quantile 
of F, (x), and B be a quantile of F,(y). It is assumed 
that the first and second partial derivatives of F(x, y) 
are continuous in a neighbourhood of («, f) and that 
f(@, B) #0. A random sample (X;,, Y,), k = 1, 2, ..., 1, 
is drawn and the values on each variate are ordered 
so that X}<Xj and Y/<Y; ifi<j. Let i and j be the 
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Correlation coefficient: 
- between quantiles is estimated by the use of an auxiliary 


greatest integers such that i/n<F;, («), j/n<F,(f), and 
let M be the number of elements (XY, Y) such that 
X<Xj;and Y<Y;. The joint distribution of (M, X;, Y}) 
is obtained and is shown to be asymptotically trivariate 
normal. The parameters of the asymptotic joint distribu- 
tion of the variates (X’, Y’) and that of the distribution 
of Q, where Q is a continuous random variable replacing 
the discrete random variable M/n such that 


| Q—M/n|<I1/n, 
are given in a theorem. 
Estimates of a, B, F(a, B) and p(a, fB) are given 
together with examples in setting confidence limits on 
these parameters. 


(H. H. Ku) 
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BIRNBAUM, ne & HEALY, on Coe en York Jniversity and Ethyl Corp.) 
Estimates with prescribed Me hala on two-stage sampling—In Be 
Ann, Math, Statist. (1960) 31, 662-676 (17 references, 3 tables\eiean gas te 


In most problems of estimation, hoa based on number Be Se the: ie ee oe 


‘method as Stein’s two-sample procedure ee ge ie par an i it 


general methods are due to Cox and Anscombe bi 
these are based on asymptotic theory and hence 


Here a new method is. siren 
conditions satisfied by 
rules for sampling in| eS s 
unbiased estimator of a given p 
equal to or not exceeding : 
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BERKSON, J. (Mayo Clinic, Rochester, Minnesota) 
Problems recently discussed regarding estimating the logistic curve—In English 


This 
Amer 
work 


Bull. Int. Statist. Inst. (1960) 37, III 207-211 


paper discusses criticism made by Silverstone [J. 
. Statist. Ass. (1959) 52, 567-577] of the recent 
by Berkson in which the maximum likelihood 


estimator for the parameters («, f) of the logistic function 
was compared with several other estimators, including 
the minimum logit y* estimator [J. Amer. Statist. Ass. 


(1953 
130-1 


) 48, 565-599; J. Amer. Statist. Ass. (1955) 50, 
62; Proc. Third Berkeley Symposium on Math. 


Statist. and Prob. (1956), pp. 1-13]. The chief points of 
Silverstone’s paper were: 


i 


am 


The maximum likelihood estimator is sufficient, 
instead of being not sufficient as had been stated. 
The minimum logit y? estimator is not sufficient 
in all conditions, though it may be in the conditions 
explored. 


. The minimum logit y? estimates are not consistent 


““in Fisher’s sense for finite samples,” (F.C.) 
even though they are consistent in the sense of 
convergence in probability (C.P.). 


. The superior efficiency (mean square error) of the 


minimum logit y* estimate was reported for certain 
values of the parameters. However, for values of 


2/647 


BERKSON, J. 
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(c) 


(Mayo Clinic, Rochester, Minnesota) 


Bull. Int. Statist. Inst. (1960) 37, III 207-211 


logit y? estimate: In the conditions explored, 


i.e. n,; = 10, k = 3, and for an infinite number 
of other controllable conditions, the estimates are 
sufficient. It is true, however, that conditions 
can be defined in which the estimates are not 
sufficient. 

With respect to consistency of the minimum logit 
y? estimate: The estimate is consistent (F.C.) as 
well as consistent (C.P.). Rao has proved the 
equivalence of F.C.=C.P. for conditions which 
are fulfilled in these experiments [Bull. Int. Statist. 
Inst. (1957) 35, 25-32]. It is remarked that the 
maximum likelihood estimate is not always con- 
sistent (F.C.) and an example is given. 


(d) With respect to the relative efficiency of the 


estimates: It appears to be considered that the 
minimum logit 7 is superior in conditions where 
the dosages are not exceedingly far from X59 but 
doubt is expressed as regards situations in which 
all P; are near unity (or zero). Actually, experi- 
ments showed the maximum likelihood estimator 
to be increasingly inferior to minimum y estimators 
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4.1 (6.7) 


the parameters ‘“‘ far outside the range ’’ that was 
examined, “it is doubtful whether such a claim 
could be upheld.” 


With regard to these points, the present paper indicates 


(a) With respect to the sufficiency of the maximum 


likelihood estimate: A necessary and sufficient 
condition for the estimates (@, f) to be sufficient 
is that no two or more samples having different 
values for the minimal sufficient statistics ( r;, 
Xr;x;) should yield the same values for (@, f). 
There are in fact, samples having different values 
for (Zr;, Ur;x;) which yield the same estimates 
(a, B); for these samples the estimate is oo for 
either or both of (@, f). This establishes that the 
maximum likelihood estimates (@, B) are not 
sufficient. The proposal of Silverstone of certain 
degenerate curves as distinguishing ‘ solutions ”’ 
for the samples concerned, even if correct, is 
considered irrelevant, for these do not provide 
distinguishing estimates (@, f). 


(b) With respect to the sufficiency of the minimum 


continued 


4.1 (6.7) 
continued 


as one moved away from the X59. There is a 
point beyond which comparison cannot properly 
be made. The maximum likelihood estimator 
yields infinity as estimates for some samples. 
The frequency of these samples increases as the 
disposition of the dosages is removed from the 
X59, and comparisons were made omitting these 
samples, but only up to the point where samples 
insoluble by maximum likelihood are not more 
than 5 per cent. of the sample population. This 
is already beyond where practical assays can 
validly be made. If comparison is to be made 
for its mathematical statistical interest in regions 
with all P; not far from unity (or zero), it will 
have to be for large sample size where the prob- 
ability of infinite estimates by maximum likeli- 
hood is negligibly small. The author believed 
that the inferiority of the maximum likelihood 
estimate found for more centrally placed dosages 
with small numbers would be confirmed. 


(J. Berkson) 


BERKSON, J., & ELVEBACK, Lillian (Mayo Clinic, Rochester, and Tulane Univ., 


4.1 (2.5) 
New Orleans) 


Competing exponential risks with particular reference to the study of smoking and 


lung cancer—In English 


J. Amer. Statist. Ass. (1960) 55, 415-428 (17 references, 3 tables) 


This paper shows that the use of the “ mortality ratio ”’ 
in comparing the numbers of deaths caused by lung 
cancer in smokers and non-smokers can be very mis- 
leading since other causes of death play an important 
role. The authors utilise the notion of ‘“ competing 
risks” and present a model for the effect of smoking 
in causing death from a specific disease in which the 
effect is measured in terms of the net probability of 
death from this disease. 

The probability of death when one risk operates alone 
is called the net risk or rate. When it operates together 
with another risk it is called the crude risk or rate. 
In order to estimate the net risk from the crude risk a 
simplifying and, in this context, reasonable assumption 
is made that each net instantaneous risk or force of 
mortality is constant over the period of observation and 
that the risks are independent. 

For the case in which the numbers of deaths as well 
as the times when the death occurred are available 
the authors derive the maximum likelihood estimates, 


99 
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called “time estimates’? of the net rates and their 
variances. For the case in which the numbers of deaths 
are available but not the times when deaths occurred, 
the maximum likelihood estimates, called “ frequency 
estimates ” and their variances are also provided. The 
entries in Table 1 show that if the crude death rate QO 
from all causes is below 0-6 the variances of the above 
two estimates are close; but for larger values of QO 
the increase of variance is substantially greater when 
the information on time of death is not used. 

Table 2 gives a list of time estimates and frequency 
estimates of the parameters involved in a model of 
exposure to two competing risks such as may be used 
in the analysis of data giving numbers of deaths among 
smokers and non-smokers due to lung cancer and to 
other causes: the formulae for variances are also given. 
Some data taken from the prospective study sponsored 
by the American Cancer Society and reported by 
Hammond and Horn (1958) provides an example. 
When the net probability of death is used to measure 


continued 


4.1 (2.5) 


New Orleans) continued 


Competing exponential risks with particular reference to the study of smoking and 


lung cancer—In English 


J. Amer. Statist. Ass. (1960) 55, 415-428 (17 references, 3 tables) 


the effect of smoking in respect of a cause of death, 
smoking is seen to have the greatest effect in increasing 
the number of deaths from coronary heart disease and 
the least effect in increasing the number of deaths from 
lung cancer. If the simple difference of the crude death 
rates is used, essentially the same conclusions are reached 
as with the preceding net rates. On the other hand, when 
the mortality ratio is used, the rate for lung cancer 
appears very much higher than the rates for other 
causes of death. 

The authors propose the use of net rates as more 
logical and state that “‘ the finding of associations of 
smoking with classes of diseases in general, not only, 
or even chiefly, with lung cancer, is of central importance 
for an interpretation of the possible biologic significance 
of these statistical studies ”’. 


(J. Gurland) 
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COHEN, A. C., Jr. 


(University of Georgia, Athens, USA) 


4.3 (2.8) 


Estimation in the truncated Poisson distribution when zeros and some ones 


are missing—Jn English 


J. Amer. Statist. Ass. (1960) 55, 342-348 (9 references, 2 tables, 1 figure) 


The author considers the estimation of the parameter 
of a Poisson distribution for the special case in which 
all zero observations are missing as well as some of the 
ones. The maximum likelihood estimating equations 
for the Poisson parameter and for the proportion of 
missing ones are obtained. A table and a graphical 
aid for estimation of the Poisson parameter are included. 
Asymptotic variances and covariances of estimates are 
obtained by inversion of the information matrix. A 
function useful in computing the variance of the Poisson 
parameter is tabulated. 

As an example, the author suitably alters data of 
Varley on distribution of gall-cells in flower-heads of 
black knapweed. 


(S. Krane) 


2/651 


COHEN, A. C., Jr. 


(University of Georgia, Athens, USA) 


4.3 (2.8) 


Simplified estimators for the normal distributions when samples are singly 


censored or truncated—ZJn English 


Technometrics (1959) 1, 217-237 (14 references, 2 tables, 3 figures) 


The author presents a scheme to estimate the mean 
and variance of a normally distributed population when 
single truncation or censoring is involved. Only 
maximum likelihood estimation is used for the cases 
with restricted sampling in one tail. Other estimating 
techniques are described from the references along with 
conditions under which these may be preferable. The 
most important feature of the paper is that only one 
simple auxiliary function is required for each of the 
restricted sample cases covered. 

A table and a graph are presented of the single auxiliary 
estimation function for singly, either right or left, 
truncated samples. One table and two graphs for varying 
ratios of the number of censored observations to the 
total sample size are given of the single auxiliary estima- 
tion function for singly, right or left, censored samples. 
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The latter table and graphs are applicable for the cases 
in which either the point of censoring is fixed or the 
number of censored observations is fixed. Examples 
utilising the graphs and tables are presented for the 
six separate cases of truncated and censored sampling 
covered in the paper. 

In addition, the variances and covariances of the 
estimates are explained for truncated samples, for fixed 
point-censored samples and for fixed sample-size 
censored samples. Examples show how to utilise this 
information. 


(V. L. Anderson) 
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CORNISH, E. A. (Division of Math. Statist., CSIRO, Adelaide) 


4.8 (3.6) 


The simultaneous fiducial distribution of location parameters in a multivariate 


normal distribution—In English 


CSIRO Diy. Math. Statist. (1961) Tech. Paper No. 8 (8 references, 2 figures) 


In this paper the author proves directly the formula 
for the simultaneous fiducial distribution of the location 
parameters in a multivariate normal distribution, 
originally conjectured by Fisher [J. R. Statist. Soc. B. 
(1954) 16, 212-213]. 

If z’ = (z, ..., Z,) is distributed in a multivariate 
normal distribution with mean ¢ = (¢,, ..., C,) and 
variance-covariance matrix £, where this has a non- 
singular estimate §, based on n degrees of freedom and 
distributed independently of z, then the simultaneous 
fiducial distribution of the parameter ¢ is of the form 


FE@+p)1f{,, C-2's 6-2) 80+ ae 
(an)#?T (4n) n |s 
The author gives details for the fiducial method in 


at 
2a 


practice. Three illustrative examples are given based 
on rainfall, altitude, latitude and longitude data, 
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(Wayne State University and Stanford University) 


pressure and temperature observations at four stations, 
and shearing strength and density of kauri trees 
respectively. In conclusion the fiducial method is 
compared with the method of confidence intervals given 
by Hotelling [Ann. Math. Statist. (1931) 2, 360-378]. 


(G. Yeo) 


4.4 (2.5) 


Estimates of bounded relative error for the mean life of an exponential 


distribution—ZIJn English 


Technometrics (1961) 1, 107-109 (1 reference, 1 table) 


If the underlying distribution of life is exponential, an 
estimation procedure is presented for mean, 0, having 
a small relative error. More precisely, an estimate of 
0 is found which is within a certain percentage, 1006, 
of @ with preassigned confidence 1—«. In life testing 
this requires observing a sufficiently large number of 
failures, r. 

For the loss function which is zero when 1—dSa/0 
<1+6 and one elsewhere, the minimax estimate a, of 
@ is found. For small 6, a is approximately equal to 
0, the maximum likelihood estimate of 0. Then for 
small 56 and small «, the central limit theorem is used 
to determine how large r must be in order that 


Pr =" <i| 21i-a. 
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It is also shown how to compute, for any preassigned 
r, the confidence level of the assertion that 


where a is the minimax estimate of @. This is important 
when r is small and asymptotic results are no longer 


valid. 


Note. See also abstract No. 2/674, 5.2. 


(A. P. Berens) 
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EPSTEIN, B. (Wayne State University and Stanford University) 


4.4 (2.5) 


Estimation of the parameters of the two parameter exponential distribution from 


censored samples—In English 
Technometrics (1960) 2, 403-406 (2 references) 


The two-parameter exponential has the density function 
f(x; 4, A) = Aexp {—A(x—A)}, where x>A, and A 
is, for example, a minimum “ guaranteed ” life span. 
In a context of life-testing, n items are placed on test 
and testing is stopped after r items have failed. The 
r failure times are x;<x,<...<x,: there is no replace- 
ment of these items. Let 
k 


T, = ¥) x,+(n—k) x,. 
i=1 

Then x, and (7,—nx,) are independent and sufficient 
for estimating A and 0 = 1/2. 

Furthermore, 4 = x,—0/n and 6 = (T,—nx,)/(r—1) 
are unbiased and minimum variance for A and 0. 

Confidence limits on @ are obtained from chi-square 
tables because 2(r—1) 0/0 follows chi-square with 
2(r—1) degrees of freedom. Confidence limits on A 
are determined from F tables since n(x,—A)/0 and 
2 (r—1)0/0 are independent and chi-square with 2 and 
2 (r—1) degrees of freedom respectively; and thus this 
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ratio is equal to n(x,—A)/0 and follows F with 2 and 
2(r—1) degrees of freedom. If 4 actually equals A, 
the ratio is unity; therefore the tabulated limits of F 
give limits on n(x,—A)/0 and hence the «-limits on A 
are easily obtained. 

. The density function of z = n(x,—A)/@(r—1) is 
g (Zz) = (r—1)/(x+1)’, z=0. These can also be used to 
find confidence limits on A. 

The paper closes with a numerical example. 


(E. H. Lehman, Jr.) 


4.1 (4.3) 


On the attainment of Cramér-Rao and Bhattacharyya bounds for the variance of an 


estimate—In English 


Ann. Math. Statist. (1959) 30, 381-388 (5 references) 


The results of this paper may be applicable in problems 
of minimum variance unbiased estimation. Suppose it 
is desired to estimate a real parameter 0, or some function 
thereof, from some data XY having a (joint) density with 
a single unknown parameter 0. The Cramér-Rao 
bound or the Bhattacharyya bounds [see Cramér, 
Mathematical Methods of Statistics (1946) Princeton and 
Bhattacharyya, Sankhyd (1946) 8, 1-32] on the variance 
of an estimator are sometimes useful for verifying 
whether or not an unbiased estimator has uniform 
minimum variance. 

The author first proves, under common regularity 
conditions for these bounds, that the Cramér-Rao 
bound is achieved if and only if the logarithm of the 
density is of the form ¢ (x) g (0)+go (9)+/(x) where t 
is the estimator and g is not constant. He extends this 
result to prove that, if the density is of the Darmois- 
Koopman exponential family, and if the variance of 
the estimator attains the kth order Bhattacharyya bound 
but not the (k—1)th, then the logarithm of the density 
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is of the form given above and the estimator is a poly- 
nomial in ¢ of degree k. In fact, the variance of any 
kth degree polynomial in ¢ will achieve the kth bound 
so that any such polynomial is a uniform minimum 
variance unbiased estimator of its expectation. 

He points out that an unbiased estimator may achieve 
the Ath order bound without being a maximum likelihood 
estimator, contrary to the case of the first order bound; 
that is the Cramér-Rao bound. He treats the example 
of a single observation x for which x0~°(c>0) has a 
negative exponential distribution on the positive real 
numbers. If 1/cis an integer k, the Ath bound is achieved 
by x;,/k!, the only unbiased estimator of 6, whereas 
the maximum likelihood estimator is x,; if c is integral 
(>1), no polynomial is an unbiased estimator of 0, so 
that no unbiased estimator of 0 can achieve any Bhatta- 
charyya bound. In fact, by Laplace transform theory, 
xT (1+1/c) is the unique, and hence uniform 
minimum variance, unbiased estimator of 0. 


(W. J. Hall) 
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GODAMBE, V. P. (Indian Statistical Institute, Calcutta) 4.9 (8.2) 
An admissible estimate for any sampling design—In English 
Sankhya (1960) 22, 285-289 (5 references) 


For a finite population consisting of N individuals, on 
each of which a metric character X is defined, a funda- 


N 
mental problem is to estimate T= )}' X, on the basis 
a= 

of a sample. In an earlier paper [J. R. Statist. Soc. B 
(1955) 17, 269-278] the author gave the most general 
definitions of the terms “‘ sample ”’ and “ linear estimate ”’. 
With that terminology, he proves that the well-known 
unbiased estimator T=(X,/P,), .Horwitz & 
Thompson’s estimator [J. Amer. Statist. Ass. (1952) 
47, 668-670], where 4 runs through all distinct units 
in the observed sample and P, is the probability of the 
Ath unit being included in a random sample is an 
admissible estimator, in the sense that there exists no 
other unbiased estimator which is uniformly superior 
pony 

The author next proves that in the class of sampling 
designs with a fixed expected number of distinct units 
per sample, that is to say one which can be interpreted 
as a fixed expected cost of sampling, the estimator T 
has its maximum variance minimised, so that it possesses 
the minimax property. 

(T. V. Hanumantha Rao) 
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GOODMAN, L. A. (University of Chicago) 4.3 (4.1) 
A note on the estimation of variance—In English 
Sankhya (1960) 22, 221-228 (10 references, 2 tables) 


In this paper the author proves that if X is a random 
variable such that & log (X | 0) is a known constant A, 
then among all statistics of the form aX, the expression 
e 4X minimises & [A (0) (log aX—log 0)*]: it is also 
unbiased with respect to the loss function 

[log f (x)—log O]*. 

A similar result was proved by the author in an earlier 
paper [Ann. Math. Statist. (1953) 22, 114-117]. The loss 
function considered there was 1 (0) [f(x)—0]*. The 
conditions of the theorem are satisfied if the density 
function of X is of the form 0 * g (x | 0), (@>0). 

The results mentioned above are applied to obtain 
best-unbiased estimates for the variance of a normal 
population. The loss functions considered are 
(/f-—0)?/o?, (log /f—log a)’, (f—o7)/o*, (f—0°) 
and Pr (| f—o? | >ko’). 


(S. John) 
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HARTER, H. L. (Aeronautical Research Laboratory) 


The use of sample ranges and quasi-ranges in setting exact confidence bounds for the 


4.4 (3.8) 


population standard deviation o. I. The range of samples from a rectangular population 
—probability integral and percentage points; exact confidence bounds for o—In English 
ARL 31, Part I (1961) Aero. Res. Labs., Wright-Patterson Air Force Base. 


Xil+48 pp. (7 references, 3 tables) 


The author discusses point estimates and_ interval 
estimates of the population standard deviation o, based 
on the sample range and quasi-ranges. Harter [Wright 
Air Development Center Technical Report 58-200 (1958)] 
has tabulated the most efficient unbiased estimates of 
o based on one quasi-range for samples of size 
n = 2(1) 100 from normal, rectangular, and exponential 
populations and on linear combinations of two quasi- 
ranges for samples of size n = 4 (1) 100 from the normal 
population. Chu, Leone & Topp [Ann. Math. Statist. 
(1957) 28, 530-531] have proposed a procedure, using 
sample quasi-ranges, of setting confidence bounds for 
the population standard deviation. Leone, Rutenberg 
& Topp [Air Force Office of Scientific Research Report 
60-408 (1960)] have obtained approximate confidence 
bounds for the standard deviation of normal, exponential, 
and rectangular populations by first applying distribu- 
tion-free methods and then imposing the distribution. 
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The use of sample ranges and quasi-ranges in setting exact confidence bounds for the 


Though the basic procedure is sound, the use of distribu- 
tion-free methods yields approximate confidence intervals 
which are much longer than the exact confidence intervals 
which are given in the present report. 

In the case of a rectangular population, the efficient 
point estimate and the most effective interval estimates 
are those based on the sample range, so it is not necessary 
to consider estimates based on sample quasi-ranges. 
The coefficients of the sample range w in the exact 
confidence bounds for the population standard deviation 
o are found by taking the reciprocals of percentage 
points of the (standardised) range W = w/o. 

The following tables for the rectangular population 
are included: 

(i) An eight-decimal-place table of the probability 
integral of the range for W = 0-01 (0-01) 3-46 
(0-001) 3-464 and sample sizes n = 2 (1) 20 (2) 40 
(10) 100; 


continued 


4.4 (3.8) 
continued 


population standard deviation c. I. The range of samples from a rectangular population 
—probability integral and percentage points; exact confidence bounds for o—In English 


ARL 31, Part I (1961), Aero. Res. Labs., Wright-Patterson Air Force Base. 


xli+48 pp. (7 references, 3 tables) 


(ii) a six-decimal-place table of the percentage points 
of the range corresponding to cumulative prob- 
abilities P = 0:0001, 0-0005, 0-001, 0-005, 0-01, 
0-025, 0:05, 0-1 (0-1) 0-9, 0-95, 0-975, 0-99, 0-995, 
0-999, 0:9995, 0-9999 for the same values of n; 

(iii) a table, to seven significant figures or six decimal 
places, whichever is less accurate, of the co- 
efficients of the sample range w in the exact lower 
confidence bounds for o for the above values of 
P and n. 


Work is now in progress on Part II of this report, 
which will give exact confidence bounds, based on 
sample quasi-ranges, for the standard deviation of a 
normal population. 


(He Harter) 
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MURTHY, M. N. & NANJAMMA, N. S. (Indian Statistical Institute, Calcutta) 


4.1 (4.9) 


Almost unbiased ratio estimates based on interpenetrating subsample 


estimates—ZIn English 
Sankhya (1959) 21, 381-392 (4 references) 


In this paper, the efficiencies of different types of ratio 
estimates built up from n independent, interpenetrating 
subsamples have been compared with regard to their 
biases and mean-square errors. The main comparison 
isp betweense they twoe types Ry = (vib i) + y,)/ 
(4 +xX2+...4+x,) and R, = (¥y/x,+-..+y,/x,)/n; the 
y;,s and x,’s being unbiased estimates of the mean values 
of y and x obtained from the ith subsample. It is shown 
that the bias of R, is n times that of R, to the second 
degree of approximation. The mean-square error, M, 
of R,, is M,+n~? (n—1) A+(1—1/n) B? where B is a 
constant; A is another constant positive for bivariate 
normal distributions, and M, is the mean-square error 
of R, to the fourth degree of approximation. 

Using these results an almost unbiased estimate 
R, = (nR,—R,)/(n—1) is constructed and conditions 
are derived for the variance of R, to be less than the 
mean-square error of R, in terms of B, A and the correla- 
tion p between R, and R,. Similarly, for some suitable 
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values of p, the minimum value of n for which the 
gain in using R, is greater than that in using R, is 
positive and the value of n for which the gain is a 
maximum are tabulated. 

With the help of an empirical study conducted with 
n = 2, the validity of the theoretical results is demon- 
strated: generalising the results further, ratio estimates 
of the type R,,, obtained by grouping the n subsamples 
into m groups, are compared taking the approximation 
to any general 2kth degree. 


(G. Parthasarathy) 


(Indian Statist. Inst., Calcutta) 4.9 (4.3) 


Some sampling systems providing unbiased ratio estimators—-Jn English 


Sankhya (1959) 21, 299-314 (5 references) 


This paper deals with sampling designs that provide 
unbiased ratio estimators in the case of sampling from a 
finite population. For example, selecting the first unit 
with probability proportional to x, which is known for 
every unit, and the rest of the (n—1) units from the 
remaining (V—1) units with equal probability without 
replacement leads to the unbiased ratio estimator y/x. 
An estimate of its variance is derived and it is shown 
that this unbiased ratio estimator is more efficient than 
the corresponding biased one which can be obtained 
by simple random sampling of units from N units, 
if the correlation between y?/x and xX is not negative. 
In the case of a stratified one-stage design, (2 N,y;)/ 
(x N;x;) becomes unbiased if one unit out of all the 
units from all the strata is selected with probability 
proportional to x; the remaining (”;—1) units from 
that stratum which contains the selected unit and all 
the n; units from the jth stratum (j #i) are selected 
with equal probability without replacement from the 
respective strata. 

As a generalisation, a necessary and sufficient condi- 
tion has been derived for estimating a general class of 
parameters F, expressible as y f(@) where a is a 
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subset of the population; A is the class of sets a, and 
f(a) a set function. F is estimable if and only if each 
sample w contains at least one set « and each « is con- 
tained in at least one w. Then 


F=[)> f(@) ¢(w, «)+Pr (w) ] 
where $(w, «) is such that )’ $(w, «) =1 for all a, 


w-@ 
is an unbiased estimator F. An estimate of the variance 
of F is also given. By taking ¢ (w, «) = Pr(w |), the 
conditional probability of getting the sample w given 
that the test « has been selected first, the usual estimators 
are obtained in many cases. In order to obtain an 
unbiased estimate of F/G, where G is another parameter, 
it is sufficient first to select the set « with probability 
proportional to g(a) and the rest of the units with any 
probability scheme. In this case, the unbiased ratio 
estimator is given by 
R= ¥ f@Pr(w|a)/ > g@Pr(w] a). 

The estimate of its variance ‘is also derived. Applica- 
tions of these generalised results are made to estimate 
population mean, variance and the regression coefficients: 
the procedure is extended to two-phase sampling. 


(G. Parthasarathy) 
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NARAYANA, T. V.. CHORNEYKO, I. & SATHE, Y. S. (Univ. Alberta, Canada) 
Sufficient partitions for a class of coin tossing problems—In English 


4.7 (4.3) 


Biom. Zeit. (1960) 2, 269-275 (9 references, 1 figure) 


The authors consider estimation problems in a Markoff- 
chain of alternatives. Let the probability of a coin 
falling head be p, (be p,), if in the previous trial the 
outcome was tails (was heads). At the first trial the 
probability of a head is p,. Two different stop-rules 
are considered: first, the coin is tossed until the total 
number of heads exceeds the total number of tails by 
r, and second the coin is tossed n times. 

(i) If the first trial results in tails, let m be the total 
number of tails and let m be the number of changes 
from tails to heads or, if the first trial results in a head, 
let n again be the total number of tails, but m—1 be 
the number of changes from tails to heads. It is shown 
that (n, m) is a sufficient statistic and the minimal 
sufficient partition in the sample-space is constructed. 
The statistic (n, m) is complete, so there exists a unique 
unbiased estimator for p,, which has minimum variance. 
This estimator and also the estimator for p, are given. 
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The analysis of life test data—In English 


(ii) Let r be the number of heads in n tosses. We 
define R = +r, if the last trial resulted in a head and 
R= —-r, if the last trial resulted in tails. Let k be the 
number of changes from tails to heads. Then (R, k) 
is a sufficient statistic and it is shown, that it is also 
complete. Finally the authors give the minimum 
unbiased estimator for p,. 


(W. Vogel) 


4.3 (3.8) 


Technometrics (1959) 1, 9-19 (15 references, 5 tables, 3 figures) 


A common problem in life testing is the estimation of 
the probability that the life, L, exceeds a specified value 
Ly. Estimation procedures are given, without derivation, 
which are in current use. An introduction is given to 
the ideas behind order statistics, with illustrative examples 
from the normal and logistic distributions. A graphical 
method is then discussed for estimating the parameters 
of these distributions from singly censored samples. 

A second quick estimation procedure for estimating 
the mean and standard deviation from a small normal 
sample involves equating the lowest and highest censored 
values to their respective expectations and solving for 
the unknown parameters; this technique has good 
efficiency for very small sample size. 

Exact techniques, developed by Gupta [Biometrika 
(1952) 39, 260-273], and expanded by Sarhan and 
Greenburg [Ann. Math. Statist. (1956) 27, 427-451 and 
(1958) 29, 75-105], give the appropriate coefficients 
for best linear unbiased estimators of « and o, based 
on observations in a small (n<15) censored sample 
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from a normal population. An alternative technique, 
applicable for n<50 which is the present limit for 
available tables of normal order statistics, is presented 
which does not depend upon the tables of Gupta and 
other writers; and which employs only tables of the 
normal distribution and tables of the expected values 
of normal order statistics. The latter technique is also 
illustrated as applied to the logistic distribution. 

Finally, an approximate technique to handle large 
(n>50) samples from a normal population is discussed, 
in which the observations are grouped into frequency 
classes. The normal equivalent deviations for each 
midpoint is taken to represent the average expectation 
of the normal order statistic, after which the technique 
described for small (n<50) normal samples is applied. 


(H. E. McKean) 
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STEIN, C. (Stanford University) 


4.1 (4.3) 


The admissibility of Pitman’s estimator of a single location parameter—In English 
Ann. Math. Statist. (1959) 30, 970-979 (10 references) 


Pitman [Biometrika (1939) 30, 391-421] has presented 
an estimator for a location parameter, based on a 
random sample with a single unknown parameter, 
which has uniformly minimum mean-square error 
among all translation invariant estimators. This 
estimator has later been proved to be minimax [Girshick 
& L. J. Savage, Proc. Second Berk. Symp. (1951) 53-73], 
without any invariance restriction; and Karlin [Ann. 
Math. Statist. (1958) 29, 406-436], under fairly strong 
conditions on the form of the density function, has 
proved it to be admissible in terms of mean-square 
error. This paper presents an alternative, and less 
restrictive proof of this admissibility. 

The method is a generalisation of one due to Blyth 
and requires essentiaily that Pitman’s estimator have 
a third moment. Thus, for example, for samples of 
size at least seven, the sample median is an admissible 
estimator of the median of a one-parameter Cauchy 
distribution. 


(W. J. Hall) 
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STORMER, H. (Siemens & Halske, Munich) 


4.4 (2.5) 


On a test for the parameters of a life distribution—Jn German 
Metrika (1961) 4, 63-77 (1 reference, 1 table, | figure) 


The author considers the distribution function 
f@ = 1—-exp {-(@/T)} 

which occurs with life-studies. On account of the 
difficulty of obtaining estimates for « and 7 in the 
usual way, the following procedure is introduced. A 
number of items (NV) are observed; as soon as an item 
fails it is substituted by another item of equal age, so 
that the number of items under observation is constant 
all the time. Let f¢,, ..., 4, be the ages of the first n 
items which fail. By considering the joint distribution 
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Uy = tq, (Uy, U2) is a sufficient statistic. The maximum- 
likelihood method leads to the estimates & = n/u, and 
Dian ON ny al. 


., t,), it is seen that with u, = 
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Finally, the variables «/@ and (7/T)* are introduced 
which have a distribution function independent of « 
and 7. The marginal-distribution of the first variable 
is closely related to the y?-distribution and 2na/@ is 
distributed as 73 (n-1): The marginal-distribution of the 
second variable seems to be in no immediate relation 
to one of the usual distribution functions. 


(J. Pfanzagl) 
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THEIL, H., & SCHWEITZER, A. (Netherlands School of Economics, Rotterdam) 


4.3 (6.0) 


The best quadratic estimator of the residual variance in regression analysis—In English 


Statist. Neerlandica (1961) 15, 19-23 (2 references) 


This article discusses the estimation of the variance 
o* of the disturbances in regression analysis under 
classical assumptions: the “ descriptive ” variables take 
non-stochastic values and the disturbances are normal, 
independent and homoscedastic. 

It is known that the maximum likelihood estimate is 
found by determining the sum of the squares of the 
disturbances estimated by least-squares and dividing 
this sum by 7, the number of observations; also, that 
the best quadratic unbiased estimator of o* is obtained 
by dividing the sum of squares of the disturbances 
estimated by least-squares by T—A, where / represents 
the number of coefficients to be estimated; i.e. the 
number of characteristic variables (+i if a constant 
term is present). It is shown that the best quadratic 
estimator is obtained by dividing the sum of squares 
of the disturbances estimated by least-squares by 
T—A+2. 


(H. Theil) 
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TRYBULA, S. (Wroclaw) 
On some loss functions—IJn English 
Coll. Mat. (1960) 7, 297-305 


Two statistical games are considered in which minimax 
estimators are found to be unbiased. 

The first deals with the multinomial distribution and 
may be defined in terms of an urn scheme as follows: 
m = (mM, ..., m,) is the result of drawing with replace- 
ment 7 balls from an urn containing balls of k different 
colours, m,; representing the number of balls of the 
ith colour in the sample. A vector function f = (fi, 
..., f,), Where f; = f; (m), is an estimator of the unknown 
composition p = (p,, .-., p,) of the urn, and constitutes 
a strategy of the statistician. The prior distribution of 
p is the strategy of his opponent. The author shows 
that if the loss function is 


1 ed eos 


na (fi— Pi)" | Pp 


then the minimax estimator is f = (m,/n, ..., m,/n), 
the game is determined and uniform distribution over 
the (k—1)-dimensional range of p is the least favourable 
prior distribution. 
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4.0 (0.7) 


The author also derives a similar result for the hyper- 
geometric distribution. Here it is assumed that an 
urn contains a fixed number WN of balls of k different 
colours, and the statistician wants to estimate how 
many balls of particular colours there are in the urn, 
knowing the composition of a sample of size n drawn 
without replacement. The author finds a loss function 
for which the unbiased estimator (Nm,/n, ..., Nm,/n) 
turns out to be the minimax estimator of the content 
M = (M,, ..., M,) of the urn, the uniform distribution 
of M over its range to be the least favourable prior 
distribution, and the statistical game in question to be 
determined. 


(S. Zubrzycki) 
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WALSH, J. E. (System Development ae Santa mone Cal) 


Nonparametric linear estimation of common median of symmetrical opulations from roe 
symmetrically censored observations—In English p one 


Sankhyi (1960) 22, 295-300 (5 references) han Sap tt | Meee NAN 


The author is concerned with the problem of estimating 
the median of a symmetrical population using the 


underlying data consist of observations whose ordered available order statistic ae 
values are denoted by x(1)S%(2)S---S%m- In ‘the 505-511 eunls! EEL 


X(nt1-r are known. 
The author considers linear estimates of the form 
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WEISSBERG, A. & BEATTY, G. H. (U.S. Food & Drug Admin., Washington, 


4.5 (11.1) 


& Batelle Memorial Inst., Columbus, Ohio) 
Tables of tolerance-limit factors for normal distributions—In English 
Technometrics (1960) 2, 483-500 (15 references, 2 tables, 1 figure) 


The aim of this paper is to present tables permitting an 
experimenter to predict that at least a given percentage 
P, of future measurements from a normal universe will 
with a specified probability y, lie within a computed 
interval X+Ks, where X and s are conventionally 
defined and K (N, P, f, 7) = ru. 

Define r (N, P) in this way: 

Y (/./N+r)—¥ (1/,/N—r) = iP, where iN is the 
sample size and P (x) is the probability that a standard- 
ised normal variate does not exceed x. 

Define u(f, ) in this way: u* =f/x;,,, u20, where 
f, not necessarily equal to N—1 is the degrees of freedom 
for s and y;, is the abscissa on the y” density function 
with f degrees of freedom such that y of the distribution 
lies to the right. 

Newton’s method, 7;,, = 7r;—[/(r)/f' (7)], was used in 
order to compute r. Here 

ff) = //N+nr)—¥ (1/,/N—r)—-P = 0. 

Table r gives r for N = 1 (1) 100 (5) 200 (10) 300 (20) 
500 (100) 1000, 2000, 3000, 5000, 10,000, oo, and 
Pa )- 50. 0-75, 0°90; 0795, 0-99:.0:999. 
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WORMLEIGHTON, R. 


(University of Toronto, Ontario) 


Calculations from Thompson’s, Hald and Sinkbaek’s, 
and Pearson’s tables were required to compute uw. 
Table u gives u for f = 1 (1) 150 (2) 250 (5) 500 (10) 
800 (20) 1000 (1000) 10,000, oo, and y = 0:90, 0:95, 
0-99. 

Calculations of r and u were performed on the IBM 650 
at Battelle with an accuracy of at worst one unit in the 
fourth decimal place, except for u(f<20, y = 0-999) 
which is uncertain. 

Now since K = ru, given N, P, F, and y we may 
uniquely determine K: were we to list tables for 
K(N, P, f, y) directly, about 1660 times as many entries 
would be required. Bowker has given a table similar 
to this but requiring f = N—1. 

The authors point out an imperfection in their tables. 
The interval X +rus really corresponds to N, P, f, and 
y’. The exact value | y’—y| is unknown, but Wald & 
Wolfowitz have given some typical limits on y. The 
authors include the estimated number of significant 
decimal places in the product ru for selected combina- 
tions of N, P, f, and y. 


(E. H. Lehman, Jr.) 


4.7 (4.4) 


A useful generalisation of the Stein two-sample procedure—In English 
Ann. Math. Statist. (1960) 31, 217-221 (1 reference, 1 figure) 


A general two-sample procedure is proposed for estimat- 
ing the mean, py, of a normal population with unknown 
variance o”. Special cases discussed include: 


(i) fixed sample size, 

(ii) Stein’s two-sample method, in which the precision 
of the estimate is fixed [“‘ A two-sample test for a 
linear hypothesis whose power is independent 
of the variance’’, Ann. Math. Statist. (1945) 16, 
243-258], 

(iii) a modified Stein procedure with an upper bound 
on the total number of observations, 

(iv) a method in which the decisions on sample size 
and precision are postponed until after the first 
sample is obtained, and 

(v) a procedure which minimises the total cost, i.e. 
the cost of the observations plus the cost resulting 
from a lack of precision of the estimate. 


The procedure is as follows: let s? be the sample 
variance of a first sample of m9 observations. Let n 
be the total number of observations, let 2L be the 
length of a (1—«) confidence interval for yu, and let 
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t™ be the (1 —«) point of a r-variate with (nj) — 1) degrees 
of freedom. Consider the family of curves, indexed 
by s, 
n(L) = (st/L)?. 

(a) On each curve choose in advance, a single point: 
this set of points constitutes a “cut” across the family. 

(6) Take the first sample and calculate s: this deter- 
mines a particular curve of the family and, on the cut, 
a unique point, (n*, L*). 

(c) Take [n*—n,]+1 further observations, where [q] 
denotes the largest integer strictly less than q. 

(d) Calculate x, the mean of all the observations: 
then x +L* is a (1—«a) confidence interval for p. 

The proof is that given by Stein. Appropriate 
definitions of the cut lead to the special cases discussed. 


(R. Wormleighton) 
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ELTEREN, P. van (Catholic University, Nijmegen) 
Two restricted sequential sign tests—Jn English 


5.7 (1.8) 


Statist. Neerlandica (1961) 15, 91-96 (3 references, 3 figures) 


During the last few years Armitage, “ Restricted 
Sequential Procedures” [Biometrika (1957) 44, 9-26] 
has stimulated the use of sequential methods which 
limit the number of observations required to arrive at 
a decision. There also exist restricted versions of the 
Classic sequential schemes of Wald. Here a sequential 
sign test of Armitage with three decision possibilities 
is compared with the corresponding Wald test which is 
further discussed by de Boer. The Wald scheme is a 
sequential test with three possible decisions for testing 
an unknown probability. It is shown that the efficiency 
of the Wald test is dependent more on the order in 
which the results are obtained than the Armitage test 
so that, in many cases, the latter will be preferable. 


(P. van Elteren) 
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EPSTEIN, B. (Wayne State University and Stanford University) 


5.2 (2.5) 


Tests for the validity of the assumption that the underlying distribution 


of life is exponential. I.—Jn English 
Technometrics (1960), 2, 83-101 (31 references) 


This article deals with a number of graphical and 
analytical procedures for testing, on the basis of life 
test data, the assumption that the underlying distribution 
of life is exponential. Part I contains descriptions of 
the mathematical and graphical procedures, while Part II 
[Technometrics (1960) 2, 167-183] contains many worked 
examples illustrating the methods of Part I. 
Included are 12 tests, labelled by Epstein as follows: 


(1) a graphical procedure, 

(2) the y” test for goodness of fit, 

(3) a criterion based on the conditional distribution 
of total lives, 

(4) a test for abnormally early failures, 

(5) a test for an abnormally long first failure, 

(6) a test for whether the mean-life, or rate of 
failure in the first half of a life test differs 

’ significantly from the mean-life or rate of failure 

in the second half of a life test, 

(7) a test for whether or not the mean-life fluctuates 
during the life test, 
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(8) a special case of (7) when the group size is unity, 
(9) a test based on the maximum F distribution, 
(10) tests for abnormally long periods in which there 
is no failure, 
(11) a graphical procedure based on the Kolmogoroff- 
Smirnoff test, and 
(12) a test based on the conditional rate of failure. 


The description of these tests is followed by general 
comments, a short summary, an appendix in which some 
of the more mathematical arguments are found, and 
finally by a bibliography. 

Some of the tests are restricted in various ways which 
are clearly noted by the author. Several remarks are 
given following each test indicating variations on that 
particular test. 

Note. See also abstract No. 2/654, 4.4. 


(R. L. Conn) 
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KLOTZ, J. (University of California, Berkeley) 
Null distribution of the Hodges bivariate sign test—In English 
Ann. Math. Statist. (1959) 30, 1029-1033 (1 reference, 1 table) 


5.6 (5.8) 


For measurement vectors (x;, y;) (xj, y}), i= 1, ..., n, 
from two bivariate distributions, Hodges [Ann. Math. 
Statist. (1955) 26, 523-527] proposed a sign test for the 
hypothesis that the bivariate distributions are the same 
and that the vectors are independent. He considered 
the vectors of differences (x;—x;, y/—y,) and the 
maximum value (= M) of the number of projected 
vectors with a given sense over all possible directions. 
The hypothesis is then to be rejected if M is too large. 
The particular combinatorial problem of determining 
the number of sign patterns was then resolved by 
identifying it with the classical gambler’s ruin problem. 
Hodges gave expressions for the number of ways m, (n) 
of proceeding from the point (0, 0) to (k, n—k) and 
then hitting the line: y = x+A only at the nth step. 
Values of Pr {K<k} for the number of patterns with at 
most k negative signs uppermost were computed for 
n = 1 (1) 30 in case k<n/3. 
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PFANZAGL, J. (Universitat K6ln) 


The present paper extends the values of Pr {K<k} 
for n = 31 (1) 50 and provides values of P for k such 
that P just exceeds 0-10. Expressions for Pr {K = k} 
are also derived in the more general case: 


(2j)+1) h<n<(2j+3) h. 


(B. M. Bennett) 


5.0 (5.1) 


On the existence of uniformly most powerful tests—In German 


Metrika (1960) 3, 169-176 (4 references) 


Let {P,, 0€Q} be a class of probability measures on 
a measurable space (X, A). Assume that each Py, is 
absolutely continuous with respect to a o-finite measure 
pt. Denote the Radon-Nikodym derivatives of Py with 
respect to p by p(x). One of the probability measures 
of the class {P,} is marked as the null hypothesis Po, 
the remaining P, form the alternative hypotheses. The 
quotient of p, (x) divided by py (x) is called the likelihood 
ratio f, (x). It is well known that the monotony of 
these likelihood ratios is sufficient for the existence of a 
uniformly most powerful test for P) against all alter- 
natives P, for every level of significance «. 

In the present paper it is proved that this condition 
of monotony is also necessary. 

For every level of significance « (0S a<1) there exists 
a uniformly most powerful test statistic ¢ (x, ), defined 
for almost every x¢ X, for the null hypothesis Po 
against all alternatives Py. Then there exists a real 
valued, measurable function T(x), depending only on 
P, and ¢, and a class of real monotonic functions H, 
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with the property that the likelihood ratios f, (x) can 
be represented <in- the “form: -f, (x) =H, T(x) for P39, 
almost every xe X and all 0€Q. T(x) is a sufficient 
statistic. The relation to a theorem of Neyman & 
Pearson [Statist. Res. Mem. (1936) 1, 113] is mentioned. 


(W. Uhlmann) 


Remark of the Regional Editor: 

There is one weak point in the proof concerning the 
independence of T(x) from the alternative; a supple- 
mentary note will appear in Metrika (1961) 4, issue 
No. 2. 


(J. Pfanzagl) 


ROY, J. (Indian Statistical Institute, Calcutta) 


Non-null distribution of the likelihood ratio in analysis of dispersion—IJn English 


5.8 (3.9) 


Sankhya (1960) 22, 289-294 (8 references, 3 tables) 


It is well known that every problem in analysis of 
dispersion, that is the generalised analysis of variance, 
can be put in the following form: a random matrix Y 
of order (pxm) of uncorrelated p-dimensional column 
vectors is given; and it is known that all the column 
vectors have the same dispersion matrix. An unbiased 
estimate V of the dispersion matrix is available. In 
order to test whether &(Y) =0, we assume that the 
joint distribution of the elements of Y is multivariate 
normal and that the estimate of the common dispersion 
matrix of the column vectors is distributed independently 
of Y according to the Wishart’s law with n degrees of 
freedom. The appropriate likelihood ratio criterion is 
A=|W|/|W+YY’| where W=nV. The null 
distribution of this criterion has been investigated 
extensively, see for instance, [C. R. Rao, Advanced 
Statistical Methods in Biometric Research (1952) New 
York: Wiley]. 
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UZAWA, H. (Stanford University, California) 


The present paper obtains an asymptotic expansion 
of the non-null distribution of A in the case where 
& (Y) is of rank one. ‘The expansion is in terms of the 
distribution function of the chi-square variable. The 
author observes that the sum of the first few terms in 
the expansion gives satisfactory results provided the 
deviations in the parameters are small and that n is 
large. A table of the power function of the analysis of 
dispersion test at the 5 per cent. level of significance is 
given for p = 1 (1).4; m = 2,3 and n = 200. 


(S. John) 


5.6 (5.3) 


Locally most powerful rank tests for two-sample problems—Zn English 
Ann. Math. Statist. (1960) 31, 685-702 (7 references) 


WeteN Ge eeA, |) X,)445s--5 4, Delindependent random 
variables where, for 1<i<n,, X; is distributed according 
to a distribution function F, and for n,+1Sisn, X; 
is distributed according to a distribution function G. 
Let Z; be zero or one according to whether the ith 
smallest member of the combined sample is distributed 
according to F or G. The author discusses the problem 
of characterising tests of the hypothesis F = G which 
are based on Z,, ..., Z, and which are locally most 
powerful in the sense that there is a one-parameter 
family of alternatives {(Fo, Gp): OSOS0,} such that 
Fo = Go and the test is uniformly most powerful against 
{(F,, Go): OS OS6_} for some 0<6) <0. 

The form of any locally most powerful rank test is 

“reject the hypothesis if }) a;Z,;>c, and accept the 
i= 14 

hypothesis if the opposite inequality holds”. The 

constants a; are given explicitly. 

Criteria are given for four problems to enable one 
to decide if a test of the above form is locally most 
powerful against some one-parameter family of alter- 
natives. In order to test F = G against the alternative 
F #G, then any test of the above form is locally most 
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powerful. Every such test is locally most powerful 
against a one-parameter family having the property 
that F,Ga‘' (x) =(1—0) x +OH (x), where H is a dis- 
tribution function on the closed unit interval. 

To test F = G against the alternative F>G, a rank 
test is locally most powerful if and only if it is of the 
stated form and, in addition, (Co, ¢;, ..., C,-2) has all 
non-negative Hankel determinants, where 


a= ([u("7 7) Gi )@n-a 


If the problem is to test F = G, symmetric, against 
the alternative F # G, symmetric with the same median, 
a rank test is locally most powerful if and only if it is 
of the stated form and, in addition 

>i (751) ues--a) =0 
s=jti\ J 
for j = 1, 2,...,. To test F = G against the alternative 


FdG>4, a rank test is locally most powerful if and 


only if it has the form mentioned with 
[(n+1)/2] 


) ((n+1)/2—s) (Qn41-s—45)20. 
ed (D. R. Truax) 


WIJSMAN, R. A. (University of Illinois, Urbana) 


5.7 (5.1) 


A montonicity property of the sequential probability ratio test—In English 
Ann. Math. Statist. (1960) 31, 677-684 (6 references) 


Based on the inequality 
aP, (a<Yy<b)SP, (a< Yy<b)<bP, (a< Yy<b) 


where Yy is the sequential probability ratio at the nth 
observation, the indices 1, 2, refer to the two hypotheses 
under consideration, and a and b are any real numbers, 
“it is shown that, if, in a sequential probability ratio 
test, the upper stopping bound is increased and the lower 
stopping bound decreased, and if the new test is not 
equivalent to the old one, then at least one of the error 
probabilities is decreased’. This result implies a previous 
result of Weiss [Ann. Math. Statist. (1956) 27, 1178-1181] 
concerning the monotonicity of the error probabilities 
as a function of the stopping bounds and generalises 
the result of T. W. Anderson & Friedman [Contributions 
to Probability and Statistics. Essays in Honor of Harold 
Hotelling (1960) Stanford: University Press] concerning 
the optimum properties of the sequential probability 
ratio test. : 

In the words of the author, ‘The monotonicity 
property is a consequence of the following stronger 
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ZAJTA, A. (Agricultural University, Budapest) 
On the Lehmann test—Jn Hungarian 


result. Let the old and new tests be given by the 
stopping bounds (B’, A’) and (B, A), respectively, with 
B<B’<A'<A; let (a, «3) and («,, 2) be the error 
probabilities and Aw; = «;—«; the changes in the error 
probabilities; then the vector (Aq,, Ac.) is restricted 
to a cone consisting of the third quadrant, plus the 
part of the second quadrant where — Aa,/Aa,<B, plus 
the part of the fourth quadrant where — Ax,/Aa,> A. 
Another consequence of this result is that (a, 2) 
cannot lie in the closed triangle with vertices (a}, «3), 
(0, 1) and (1, 0)”. 

“* Finally, the following monotonicity property follows: 
if the lower stopping bound is fixed and the upper 
stopping bound increased, then «,/(1—a,) decreases 
monotonically. The same holds for «,/(1—«,) if the 
upper stopping bound is held fixed and the lower stopping 
bound decreased.” 


(C. A. Bennett) 


5.3 (5.6) 


Publ. Math. Inst. Hung. Acad. Sci. (1960) 5, 447-495 (3 references, 1 table) 


Lehmann, see ‘‘ Consistency and unbiasedness of certain ~ 


non-parametric tests [Ann. Math. Statist. (1951) 22, 
165-180] and Rényi [Magy. Tud. Akad. Alk. Mat. Int. 
K6ézl. (1953) 2, 243-265], has modified the Wilcoxon 
test into a form which is consistent against all alternative 
hypotheses. Referent, ‘On two modifications of the 
Wilcoxon test” [Publ]. Math. Inst. Hung. Acad. Sci. 
(1959) 4, 313-319], has shown the identity of these 
modifications to be a simple linear relationship between 
their statistics. 

In this paper a new proof of this connection is given 
and a new expression of the statistic is established. 
This new expression is, in certain respects, simpler than 
the earlier ones as it is only a linear function of the rank 
numbers. 

The variance of the statistic is computed in the general 
case of alternatives. As a special case the variance in 
the case of null hypothesis is also obtained. 


(E. Csaki) 
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BASMANN, R. L. (Hanford Laboratories Operation, G.E.C., Richland, Washington) 


6.6 (6.1) 


An expository note on estimation of simultaneous structural equations—In English 
Biometrics (1960) 16, 464-480 (27 references, 2 tables) 


The purpose of this paper is to make known to biologists 
and biometricians several methods of generalised linear 
estimation of structural parameters not yet widely 
known outside econometrics; in particular, the least- 
variance difference method. The technique of least- 
variance difference estimation is formulated exactly as 
the generalisation by Aitken [Proc. Roy. Soc. Edinb. 
(1934-35) 55, 42-48] of Gauss’s method of treating 
independent observations of unequal precision to the 
case of interdependent observations. 

The author begins with a discussion of some logical 
features of a simple causal path model of the type 
presented by Turner & Stevens [Biometrics (1959) 15, 
236-258: abstracted in this journal No. 1/258, 6.6]. 
The model is expressed functionally in two equivalent 
ways, reduced form and structural form. The hypothesis 
that the structural equations are linear in the unknown 
parameters and some tests that might be made are 
considered. The author then discusses three of the 
several ways one may obtain the least-squares estimates 
of the structural parameters which are needed, these 
methods differing only in the degree to which estimates 
of the coefficients involved are constrained to satisfy 
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described are conceptual only; 


(Hanford Labs. Operation, G.E.C., Richland, Washington) 


the hypothetical restrictions deduced from the structural 
model. It is emphasised that the estimation processes 
the computational 
processes, however, are equivalent to the conceptual 
substitution procedure used. The three methods con- 
sidered in this paper lead to the full-information 
estimates, the least-variance ratio estimates using the 
limited-information single equation method, and what 
the author calls the least-variance difference estimates. 

In continuing, the author describes the heuristic 
considerations that led to his initial formulation of the 
least-variance difference method, these considerations 
being inspired by the classical Gauss-Fisher least-squares 
problem. This formulation involves nothing more than 
a slight generalisation of the problem that motivated 
the Gauss-Legendre formulation of the method of 
least-squares itself and could serve as an introduction 
to the broader subject of this article. 

A thorough discussion of the statistical properties 
of the estimates yielded by the three methods follows. 
The paper ends with a numerical example of the least- 
variance difference method using artificial data and a 
test of the hypothesised relation by two structural 
equations. (S. R. Knox) 


6.6 (4.8) 


On the asymptotic distribution of generalised linear estimators—Jn English 


Econometrica (1960) 28, 97-107 (22 references) 


This paper gives a sufficiently general proof of the 
asymptotic normality of the distribution of the general- 
ised classical linear estimators of a single equation 
in a system of structural equations: see also T. W. 
Anderson & Rubin [Ann. Math. Statist. (1949) 20, 
46-63 and (1950) 21, 570-582]. The author here makes 
the assumption that the random structural disturbances 
and the reduced-form disturbances are not normally 
distributed. The generalised classical linear estimators: 
for example the minimum variance-ratio estimators are 
not maximum-likelihood estimators. The effect of this is 
to require proof of the property of asymptotic normality 
under non-normal assumptions. 

The second section of the paper states the stochastic 
and structural assumptions: use is made of matrix 
methods here. The next section deals with linear 
decompositions of the general classical linear estimators. 
After dealing with the asymptotic normality of two of 
the estimators in section 4 the final section considers 
a large sample test of identifiability. 


(W. R. Buckland) 
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On a contingency-table with missing data—In German 


6.9 (4.3) 


Biom. Zeit. (1960) 2, 236-243 (3 references, 3 figures) 


Let f;; be the observed frequencies of the contingency- 
table, where the first index j denotes the row and the 
second index i the column of the table. The total 
number of columns is denoted by k and the number of 
rows by /. The frequencies of a column are supposed 
to be proportional up to random deviations; that is, 
say homogeneity exists between the rows, and may 
represent the marginal probability p;. Further, it is 
supposed that in the first row only m frequencies, in 
the second only m+n, etc. are available, where 
m<m+n<...<k. The problem is to evaluate the 
maximum likelihood estimates of the marginal prob- 
abilities p;. 

The author calculates these estimates and points out 
that they are proportional to the column sums of the 
available frequencies. The unbiasedness is proved 
for the case of k = 3 columns and / = 2 rows only, 
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providing the last frequency in the first row (= f;3) is 
missing. For the same special case asymptotic formulae 
for the variances and covariances of the estimates are 
obtained by inversion of the information matrix 
associated with the likelihood function. 


(B. Schneider) 


6.1 (5.2) 


Tests of equality between sets of coefficients in two linear regressions—IJn English 


Econometrica (1960) 28, 591-604 (10 references) 


In this paper the author presents systematic tests which 
are involved when, having estimated a linear regression 
with p coefficients, it is necessary to test whether m 
additional observations belong to the same regression. 
The author relates the prediction interval for m = 1 
with the analysis of covariance for m>p within the 
framework of the general linear hypothesis for any m. 
These results are then extended in order to test the 
equality between subsets of coefficients. 

Firstly, the author discusses various applications of 
normal linear regression in the representation of an 
economic relationship and points out some of the 
complications and the methods used to overcome them. 
His first test concerns the hypothesis that both samples 
belong to the same regression and uses the well-known 
prediction interval, see Mood [Introduction to the 
Theory of Statistics (1950) New York: McGraw-Hill, 
pp. 304-305]. This can be utilised when the number m 
of observations in the second sample equals one; the 
analysis of covariance test, see Kendall, M. G. [The 
Advanced Theory of Statistics (1946) 2, London: 
Griffin, p. 242], may be used when m>p. The author 
presents two tests for the case where 2SmSp. 
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The first test is based on a prediction interval for the 
mean of m observations and the second test is an 
F-test (Section 3). Later the author explains the relation- 
ship between this F test, the prediction interval and the 
analysis of covariance and proceeds to test the equality 
between regression coefficient subsets in both the 
regressions. Two examples are given from the field of 
economic relationships concerned with the temporal 
stability of a statistical demand function for car owner- 
ship and of a statistical demand function for new cars. 
The author summarises his results as follows: “To 
test the equality between sets of coefficients in two 
linear regressions, we obtain the sum of the squares 
of the residuals assuming their equality, and the sum of 
squares without assuming their equality. The ratio 
of the difference between these two sums to the latter 
sum, after adjustment for the corresponding degrees of 
freedom, will be distributed as the F-ratio under the 
null-hypothesis. This latter sum of squares will be 
computed only from the first sum of 1 observations when 
the second sample is not large enough for computing a 
second regression.” 


(W. R. Buckland) 


CHRIST, C. F. (University of Chicago) 


6.6 (11.7) 


Simultaneous equation estimation: any verdict yet ?—Jn English 


Econometrica (1960) 28, 835-845 (7 references) 


This paper is one of four papers presented by various 
authors in (1958) at the meeting of the Econometric 
Society. It is a revised version of the original paper 
and here the position taken is slightly more in favour 
of simultaneous equations methods than that taken 
previously; [see also abstracts No. 2/688, 6.6, No. 
2/690, 6.6 and No. 2/692, 6.6]. The relative merits of 
ordinary least-squares estimation and several sim- 
ultaneous equations methods are discussed and the 
importance of the idea of a system of simultaneous 
structural economic equations is stressed. The available 
evidence suggests that on occasions while economical 
least-squares methods may give results which are as 
good or better than the simultaneous equations method 
it is preferable to use the latter. 

After giving some standard definitions and remarking 
that Monte Carlo experiments are more often applicable 
as there is yet little small-sample theory of the properties 
of the different estimators, the author begins by discussing 
over-identified linear stochastic models. A model is 
stated and methods for obtaining its reduced form 
given. The estimation methods discussed for the equa- 
tions of the model, termed “ structural equations ”’ are: 
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(i) ordinary least-squares, (ii) two stage least-squares, 
(iii) limited information, (iv) full information maximum 
likelihood. Estimation methods are also proposed for a 
reduced form equation. 

A discussion of the literature concerning the equations 
and their properties follows and other related evidence 
is commented on at some length. The evidence provided 
by Monte Carlo experiments is stated and various com- 
parisons made. The author then reports on three 
Monte Carlo studies by Wagner [Econometrica (1958) 
26, 117-133]: Basman [A paper presented at the 
Econometric Society Meeting (1958)] and Summers [A 
paper presented at the Econometric Society Meeting 
(1958)] and comments on work by Nagar [Econometrica 
(1959) 27, 575-595: abstracted in this journal No. 2/300, 
4.8, and Econometrica (1960) 28, 573-590: abstracted 
in this journal No. 2/777, 11.7] and by Theil [Statistical 
Estimation of Simultaneous Economic Relationships (1958) 
Amsterdam: North Holland Publishing Co.] and states 
that Nagar has done further experiments with Wagner’s 
(1958) models. In summing-up the author says it is not 
yet clear that the least-squares method for structural 
equations is dead or should be discarded. 

(W. R. Buckland) 


6.2 (6.9) 


Publ. Math. Inst. Hung. Acad. Sci. (1960) 5, 311-323 (8 references, 1 figure) 


The bivariate stochastic connections are investigated 
with the aid of the theory of Hilbert spaces. Let 
(Q, S, P) be a probability field, L? the Hilbert space of 
the random variables € with finite variance, S- the 
smallest o-algebra for which ¢€ is measurable and 
ie = 1(O7S,, P); further 


i = {0 af PAF (s)<coh 


Li, denotes the subspace of Lz corresponding to zero 
expectation. Let A-¢ be the orthogonal projection of 
C to Lz that is, As¢ = M(¢|é). Then the measures 
of the stochastic connection, the correlation coefficient 
R(é,, €>), and the correlation ratio ©, (¢) may be 
written as R(G, 62) =(66 ): O.(0 = []Acl* 
where ¢* is the standardised form of ¢. For the maximal 
correlation S and for the contingency C the following 
relations hold: 
S(é,n) = || Az ll =IIAg lls CC 1) = Ill Ag Ill = II Ay I 
provided that the domains of A; and A, are restricted 
to L,, and L,. 
Let A, be the operator, corresponding to A,, mapping 
L2,, into L3,,; H, and H, are the density functions of ¢ 
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and y, respectively, and H (x, y) is their joint density 
function. Similarly the operator A, corresponding to 
A; can be defined. A, and A, are integral operators 
if and only if the measure generated by H is absolutely 
continuous with respect to the product of the measures 
P, and P, generated by H, and H, respectively. With 
the aid of a theorem, the above expression of the con- 
tingency can be proved to be equivalent to the definition 
given by Rényi, “ New version of the probabilistic general- 
isation of the large sieve ” [Acta Math. Acad. Sci. Hung. 
(1959) 10, 441-451: abstracted in this journal No 1/252, 
6.0]. 

Replacing the operator A, by A = A,A, instead of 
the measure P a symmetric measure P is obtained. 
H (x, y) being the distribution function corresponding 
to P we have 


In addition the square of the maximal correlation of 
H (x, y) equals the maximal correlation of H (x, y). 


Computations will be often facilitated by using this 
connection. 


” f(x) 9 (y) dH (x, y) = (Ag Ai f(s 9 (0). 


— © 


(P. Bartfai) 
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CSAKI, P. & FISCHER, J. (Math. Inst., Hungarian Academy of Sciences, Budapest) 


6.2 (6.0) 


Contributions to the problem of maximal correlation—ZJn English 
Publ. Math. Inst. Hung. Acad. Sci. (1960) 5, 325-337 (5 references, 2 figures) 


Preserving the notations of authors’ preceding paper, 
see “‘ On bivariate stochastic connection ” [Publ. Math. 
Inst. Hung. Acad. Sci. (1960) 5, 311-323; abstracted in 
this journal No. 2/686, 6.2] let the domain of A, be Li, 
then the functions fel. geL;, form a pair of eigen- 
functions with eigenvalue 4, if A.g = Af and A, f = dg. 
If f and g are standardised variates then the following 
three statements are equivalent: 


(i) f and g form a pair of eigenfunctions, 

Gi) ©-(9) = 9,(f) = (49), 

(ui) f and g are linearly correlated and ©; (g) = ©. (g) 

O7 (Ff). = ©, (Ff). 

In the remaining part of Section one the case 
S(é, 7) = 1 is considered and it is proved that the 
linearly correlated standardised variates — and y are 
homoscedastically correlated if and only if €*—1 and 
n’—1 form a pair of eigenfunctions with eigenvalue 


(e7)7: 
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Simultaneous equations: any verdict yet?—Jn English 


Econometrica (1960) 28, 846-854 (25 references) 


This paper is the second in the series presented at the 
Symposium [see also abstracts No. 2/685, 6.6, No. 2/690, 
6.6 and No. 2/692, 6.6]. The author discusses his 
reasons for believing that the experience with applica- 
tions of existing simultaneous equation procedures 
does not yet furnish a reliable proof of their practical 
usefulness. 

He gives an example of a situation which is not yet 
covered by these studies and gives extensive references 
as to the work already completed. While deploring 
such statements as “‘ Least-squares methods are obsolete” 
he remarks that at present such statements are un- 
justified and can present a biased view of the basis upon 
which to choose methods and models. 

He discusses both changes in the type of model and 
work being done on a new estimation procedure 
[Econometrica (1959) 27, 302] and stresses the import- 
ance for the future of the choice of the best procedure 
to employ in various circumstances. The author points 
out that even minor generalisations of traditional models 
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The computation of the maximal correlation is 
facilitated for some cases if {¢,} and {,} are linearly 
independent systems in Lz, and in L;,, respectively, 


n 


then the functions f,= )) ond, and g,= >) a W, 


are eigenfunctions if and only if, A,¢, and Az, can be 
written as linear combinations of ,, Wo, ..., W, and 
#1, Po, +++» b, respectively; the eigenvalue can easily 
be calculated from the coefficients. If the eigenfunctions 
are polynomials then the choice of ¢, = &"—M (é") 
is suitable. 

The second section of this paper gives examples for 
the determination of the maximal correlation. 


(P. Bartfai) 


6.6 (11.7) 


have led to severe mathematical complications. This 
has lead to difficulties in the development of procedures 
and made almost non-existant finite sample distribution 
and test theory. The case in which the Monte Carlo 
method fills an important gap is stated generally. The 
importance of testing predictions is stressed and the way 
in which the Monte Carlo method can be used in this 
context is discussed. 

Finally, the author remarks that he would not wish 
to give a verdict on simultaneous equations at this 
time but would like to see theoretical and empirical 
investigations continued. 


(W. R. Buckland) 


JOHNSON, E. A. & BROWN, B. W. (School of Public Health, Univ. Minnesota, USA) 


6.8 (4.3) 


The Spearman estimator for serial dilution assays—ZIn English 
Biometrics (1961) 17, 79-88 (9 references, 1 table, 2 figures) 


To estimate the density of a specified organism in a 
given suspension, a common method is to take n doses 
of constant volume from each of a series of successive 
dilutions of the original suspension and to observe the 
proportion of doses which show growth of the organism. 
The dilutions may be so chosen that the logarithms of 
the resulting concentrations are equally spaced, and 
in this case Spearman has proposed an estimator for 
the logarithm of the density of organisms in the original 
suspension which is a linear function of the original 
log-concentration, the difference between successive log- 
concentrations, and the observed proportion of doses 
showing growth. 

Assuming a Poisson distribution of organisms per 
dose in the original suspension and assuming such a 
range of concentrations and such a dosage that it is 
very probable that growth will occur in a dose from 
the original suspension and very unlikely in a dose from 
the lowest concentration, the authors -show that the 
Spearman estimator has negligible bias and find its 
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approximate variance which is of order 1/n. Based on 
these results and the fact that the estimator is asymp- 
totically normal, confidence limits for the log-density 
are proposed using the normal approximation. The 
authors also find the approximation to the mean and 
variance of the estimator of the density derived from 
Spearman’s estimator of the log-density. 

Two graphs are provided to facilitate finding the 
point estimate and interval estimate of the density. 


(C. C. Thigpen) 


6.6 (4.1) 


Single equation versus equation system methods of estimation in economics—ZJn English 


Econometrica (1960) 28, 866-871 (3 references) 


This paper, the fourth of the symposium (see also 
abstracts No. 2/685, 4.1, No. 2/688, 4.1 and No. 2/692, 
4.1) is a discussion by the author on the history of the 
various developments in the methods of estimation and 
includes his comments on various points such as bias, 
misuse of the choice of method to be used for estimation 
and the way in which the use of the computers of various 
sorts has influenced the field of estimation. He points 
out that limited-information estimates can now be 
obtained in a matter of minutes and high speed pro- 
grammes are being made available for the full maximum- 
likelihood estimation. Recent interesting Monte Carlo 
results have also been made possible by the advent of 
high-speed computers. 

When comparing the two methods of estimation the 
author remarks that a comparison coefficient by co- 
efficient of structural estimates is often not a suitable 
approach and he sees the need for comparisons of 
summary statistics of a whole system. Reduced-form 
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coefficients may differ widely even though component — 
structural coefficients appear closely related. The 
efficiency properties of single-equation least-squares 
estimates do not hold under the transformation from 
structural to reduced-form coefficients. Finally the 
author discusses the applicability of least-squares. 


(W. R. Buckland) 


KLOEK, T. & MENNES, L. B. M. (Netherlands School of Economics) 


6.6 (6.4) 


Simultaneous equations estimation based on principal components of 


predetermined variables—Jn English 


Econometrica (1960) 28, 45-61 (7 references, 6 tables) 


Estimation of economic equation systems is usually 
based on rather short time series: if such a system 
happens to be large this gives rise to difficulties in all 
cases where the estimation method used requires 
estimates of the moments of reduced-form disturbances. 
These difficulties often arise in cases of limited-in- 
formation, maximum-likelihood and two-stage least 
squares estimations. In this paper the author demon- 
strates a practical solution of solving these difficulties 
by means of the use of a small number of principal 
components of predetermined variables which are used 
to represent the complete set of the predetermined 
variables. 

In the first part the author discusses over-identifica- 


tion and undersized samples and shows the normally 


adopted procedures. An account of some experiments 


with the introduction of principal components of pre- 
least-squares 
residuals is given. The author says that although these 


determined variables and of certain 


approaches do not avoid all disadvantages they do 
tackle one procedure with a certain amount of success. 
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The third section gives the common part of the 
alternative estimation procedure and the subsequent 
sections discuss the use of four variants. These are 
applied numerically in the last section of the paper. 

One of the six tables used to illustrate this paper gives 
the alternative of Klein’s model I, [Klein, L. R. & 
Goldberger, A. S. ‘An econometric model of the 
United States from 1929-1952 ” (1955) Amsterdam]. 


(W. R. Buckland) 


6.6 (4.8) 


‘ Under-identification, structural estimations, and forecasting—ZJn English * 


Econometrica (1960) 28, 855-865 (13 references) 


In a system of simultaneous linear equations a structural 
relationship is identifiable only if the following condition 
is met: 

Rank [TM, 44] =G*-1; K**2G°—1 where [II,, «+] 
is the matrix of the reduced-form of the regression 
coefficients of the jointly dependent variables included 
in the relationship given (G*) on the predetermined 
variables which are included in the structure but excluded 
from the given number of relationships, K**. This 
condition appears unimportant but is necessary for 
identification: it is satisfied in all existing statistical 
estimates concerning structural relationships. Both the 
original supposition underlying the simultaneous equa- 
tion approach and, in the empirical models, the universal 
fulfillment of the necessary condition of identification, 
present an interesting situation. Economic variables 
are considered as mutually dependent but the degree 
of simultaneity is only recognised in so far as it does 
not prevent identification of the structural coefficients. 
Literature already available answers that we are often 
required to exclude a sufficient number of variables 
from the structural relationship either by theory or by 
a priori information. When considering economic 
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theory, this paper suggests, one requires to include a 
much greater number of variables than is usual in 
existing structural models for economics. 

The complexity of the modern economic structure 
makes the true structural relationship more likely to 
be under-identified rather than over-identified. Deter- 
mined efforts to obtain “‘ significant ”’ structural estimates 
may well be the cause of the fact that all existing empirical 
structural relationships are over-identified. The respons- 
ibility for the over-identification as a result of over- 
simplification is due more to statistical difficulties than 
to economic theory or a priori information. 

The author then proceeds to discuss and demonstrate 
the misleading nature of the type of so-called structural 
estimate. 

Finally the paper deals with the rather doubtful 
prospects of over-identifying an otherwise under- 
identified structural relationship by resorting to more 
division of the time periods or by disaggregating the 
variables. He does not think that the results of either of 
these procedures are of significant use. (See also abstracts 
No. 2/685, 6.6, No. 2/688, 6.6 and No. 2/690, 6.6.) 


(W. R. Buckland) 


MANN, H. B. (Ohio State University, Columbus, Ohio) 


The algebra of a linear hypothesis—In English 
Ann. Math. Statist. (1960) 31, 1-15 (3 references) 


James introduced the concept of the relationship algebra 
of an experimental design [Ann. Math. Statist. (1957) 
28, 993-1002], and demonstrated that the important 
properties of an experimental design, including the 
analysis of variance appropriate to it, can be revealed 
by analysing the structure of its relationship algebra. 
Essentially, this algebra is generated by the basic relation- 
ships that exist between the experimental units of the 
design. 

As an illustration, James presented a detailed analysis 
of the relationship algebra of a balanced incomplete 
blocks design. Mann has considered a more generalised 
algebra in which the relationship algebra of an experi- 
mental design is a special case. He illustrates his methods 
by decomposing the algebra of an s-dimensional cubic 
lattice into its principal components. 


(B. S. Pasternack) 
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6.1 (9.1) 


Programming Fisher’s exact method of comparing two percentages—ZJn English 
Technometrics (1960) 2, 103-107 (1 reference, 1 figure) 


By re-writing the expression for Fisher’s exact method 
of comparing two percentages, a high-speed computer 
programme is developed for computing these exact, to 
a degree of decimal accuracy, probabilities. The differ- 
ence in two percentages is expressed first in terms of 
the usual 2 by 2 table with entries: a, b, c and d where 
a is the min (a, b, c, d). The probability of interest is 
the sum of this observed percentage difference or per- 
centage more diverse from these such as (a—1, b+], 
c+1, d—1) etc. Fisher proves that the probability of 
the observed results when the hypothesis of no real 
difference in population percentage is true is given by: 


_ (a+b) !(e+d) \(a+c)!(b+d)!: 
~— albic!d!(at+b+c+d)! 


Robertson shows that this may be written as: 


Py 


i | 


(a+b)—-x 


P, = J] ——— 
* <20(atb+c+d)—x 


a=t (at+ec—y “  (ct+d)-z 


Vaya. z= 0 (b+c4-d)—z 


This expression lends itself to a three-loop programme 
given in the paper where x is increased from 0 to a—1, 
then y from 0 to a—1 and finally z from 0 to c—1. The 


6.9 (11.5) 


whole programme must be repeated for more diverse | 


percentage differences than the one observed, but one 


can include a preset significance level such that the — 


programme will stop when this value is exceeded. 


(C. R. Hicks) 


SANDEE, J. (Central Planning Bureau, The Hague) 
A development model for India—Jn Dutch 
Statist. Neerlandica (1961) 15, 189-210 (4 tables) 


A model for the economic development of India in 
the 1960-70 period is described from the statistical 
point of view. The model contains 13 sectors, each 
with their specific products for which balance-equations 
are set up. Input-output and capital-output relations 
connect the sectors. Agricultural output is linked to 
the use of fertiliser, irrigation and advice. Consumption 
in 1970 is maximised within thirty-five constraints. 
Finally, the author discusses the effects of statistical 
uncertainties and possible improvements of the model. 


(J. Sandee) 
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6.0 (6.6) 


6.1 (5.1) 


Tests for regression coefficients when errors are correlated—In English 
Ann. Math. Statist. (1960) 31, 929-938 (7 references) 


The author is continuing some earlier work in the 
study of regression problems in which the residual 
errors are not independent, see “‘ Covariances of 
least-squares estimates when residuals are correlated ”’ 
[Ann. Math. Statist. (1958) 29, 1251-1256; abstracted in 
this journal No. 1/92, 6.1]. In the present paper he 
investigates the results of using the usual test statistics 
for regression coefficients, derived for independent 
residuals, when this condition is not fulfilled. He points 
out that it may be possible to transform the problem 
into the independent case, but such a move is not 
always wise or feasible. 

With a specific covariance matrix for the errors, 
namely one in which elements more than one removed 
from the main diagonal are negligible, and with the 
assumption that the residual sum-of-squares and vector 
of parameter estimates may be considered independent, 
he obtains the distributions of certain test statistics, 
and relates these to their common counterparts. From 
here it is possible for him to point out that the usual 
statistics may be effectively used in many areas where 
the residuals are correlated. 
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The author gives applications of the theory to the 
cases of testing a single sample mean, the difference 
between the means of two samples, and the coefficient 
in a linear trend and in regression on trigonometric 
functions. A concluding statement is given indicating 
the meaning of the work presented and pointing out 
directions for possible extension. 


(D. H. Shaffer) 


SCHNEEBERGER, H. (Wiesbaden) 


6.1 (7.0) 


Theory and applications of tests on the hypothesis of linearity in total 


regression equations. IIl—Jn German 
Math. Tech. Wirtschaft (1960) 4, 175-179 


The second part of the author’s paper on tests of 
linearity contains in the first place the demonstration of 
the second part of the first theorem (abstracted in this 
journal No. 2/107, 6.1). It is proved by the statement 
of &s{>&sz assuming H, : &y;; = «; +f; (x;;—X), that 
is to say, not all the straight lines received for T>-type 
groups and/or not all the means received for T,-type 
groups coincide with the straight line 6y = «+f (x—x) 
generally valid on Ho. 

It is shown that the standard technique to test the 
linearity in total regression equations follows from the 
first theorem in the special case of m, = O and therefore 
m =m). Jtis also shown that theorem | is a generalisa- 
tion of a fundamental theorem in analysis of variance, 
where Hy : &y; = « for each j. 

In chapter two of the paper, theorem 1 is generalised 
to theorem 1’ by weakening the third assumption to 
é (v,;-&y))? = € y—6y)"/9;; g;>0; originally in 
theorem 1: &(y—@€y)* = const for each x. Now, in 
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Theory and applications of tests on the hypothesis of linearity in total 


regression equations. IIl—Jn German 
Math. Tech. Wirtschaft (1960) 4, 175-179 


If H, cannot be rejected, central confidence intervals 
may be specified for « and B by « = atts/,/ ())9;1;) 
j i 

with ne 

X =) 9; > xij! dg jn; and so = v9; Y Wy Yin) . 

j J I i 

The author states that he intends to publish a further 
generalisation of his theorems in which a single in- 
dependent variable x will be replaced by a set of k 
independent variables x,, x2, .--, X,- 


(R. K. Bauer) 
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the statistic F of theorem 1, s{ is to be replaced by 
(51) =) 93D iyi")? (m, + 2m, —2) and s3 by 
A] L 


(53)? = ye gj yy, (Vij)? /n- m,—2m). 


All other conditions and results remain. Theorem 1’ 
is a very theoretical one, for the g,’s will never be known 
in practice, and their estimates will always be biased. 
As a by-product of theorem 1’ the author obtains a 
second theorem to test the identity of m regression 
coefficients «; and f£, simultaneously. He interprets 
the m = m, T,-type groups (m,; = 0) as stochastically 
independent sets of x values and tests the simultaneous 
hypothesis Hj: «a, =a and f£;=f for each j by 


F = (s{)7/(s3)” with (s{)* = » 95» (is vii*)?/2 (m—1) 
J i 
and (5%)? = yy gj Ne (vi; —Yj*)?/(n—2m), F-distributed 


j i 
with f; = 2(m—1) and f, = n—2m degrees of freedom. 


continued 


6.1 (7.0) 
continued 


ABT, K. (U.S. Naval Weapons Lab., Dahlgren, Virginia) 


7.6 (6.1) 


Analysis of covariance and analysis of differences. I1—Jn French 
Metrika (1960) 3, 177-211 (7 references, 1 table, 3 figures) 


The first part of this study, see [Metrika (1960) 
3, 26-45 and 95-116: abstracted in this journal No. 
2/111, 7.6] dealt with the analysis of covariance with one 
independent variable. The second part discusses the 
relations between the ‘analysis of covariance and the 
analysis of variance of differences. 

Referring to linear regression analysis it is first shown 
that the generalised differences d = y—Box, where By 
is postulated, can be properly used for measuring the 
increment, adjusted if By #1, of a variable if y is 
linearly dependent on x and if By equals the true value 
of the regression coefficient of y on x. Here x and jy, 
respectively, stand for the values of the variable before 
and after a “treatment ”’ being applied to the experi- 
mental unit on which d is measured or observed. The 
test of the null hypothesis & (6) = Bo, therefore, is 
referred to as “‘ Test for Admissibility of Differences ”’. 

Then, for the case of t groups of paired observations 
‘before and after treatment”, the power of the test 
for the null hypothesis “‘ no group effects ’”’ in the one- 
way analysis of variance of generalised differences is 
found to be at least equal to that of the analysis of 


\ 
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covariance if again &(b) = By; b this time being the 
regression coefficient for error “‘ within-groups ”’. The 
F-test in the analysis of differences is even more powerful 
if besides & (6) = By appreciable differences between the 
group means in x exist. This applies especially to the 
case of ¢ = 2 groups. 

For computational purposes the F values of the two 
types of analyses themselves are compared. For the 
special but most important case of ‘‘ ordinary” differ- 
ences (Bf) = 1) the analysis of covariance can be by- 
passed if in fact the F-test in the comparatively simple 
analysis of differences is of equal or higher power as 
can be concluded from certain procedures described in 
the paper. Numerical examples are given to illustrate 
the findings. 

The results found for the one-way classification are 
then generalised for the cases of two and more ways of 
classification of the data. 


(K. Abt) 


Confidence bounds connected with ANOVA and MANOVA for balanced and 
partially balanced incomplete block designs—Jn English 


Ann. Math. Statist. (1960) 31, 741-748 (5 references) 


Simultaneous confidence bounds on particular sets of 
parametric functions associated with a strongly testable 
linear hypothesis for either ANOVA (analysis of 
variance—univariate) or MANOVA (multivariate analysis 
of variance) have been obtained by S. N. Roy & J. 
- Roy [/nst. Statist., Univ. North Carolina Mimeo. Series 
No. 208, (1958)]. These parametric functions can be 
generally regarded as measures of deviation from the 
“total”? hypothesis and its various components, and 
are such that for any specific problem one of them can 
be appropriately designated the “total” parametric 
function while the others are “ partial”? parametric 
functions of various orders. The “total” function, 
firstly in the univariate case, is related to the non- 
centrality parameter of the usual F-test of the “ total ” 
hypothesis in ANOVA, and secondly, in the multi- 
variate case, is the largest characteristic root of a certain 
parametric matrix which is related to another parametric 
matrix whose non-zero characteristic roots occur as a 


2/700 


set of non-centrality parameters in the power function 
for the test, regardless of which of the standard tests is 
used, of the total hypothesis in MANOVA. This | 
statement also applies to “ partial’ functions of various 
orders and their related “‘ partial ” hypothesis. 

In this paper for both ANOVA and MANOVA, 
the author considers the hypothesis of equality of 
treatment effects and vector equality in the multivariate 
case. He then proceeds to develop explicit algebraic 
expressions for the total and partial parametric functions 
that are attached to the simultaneous confidence state- 
ments subsumed under the hypothesis of equality of 
treatment effects, both for ANOVA and MANOVA. 
The author later indicates the procedure for obtaining 
algebraic expressions for the confidence bounds on each 
such parametric function, without a derivation of these 
expressions in an explicit form. 


(B. S. Pasternack) 


7.1 (9.1) 7 
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COX, D. R. (Birkbeck College, London) 


8.8 (1.8) 


Serial sampling acceptance schemes derived from Bayes’ theorem—ZIn English 
Technometrics (1960) 2, 353-360 (4 references, 3 tables) 


Since consecutive batches of a sampling scheme may be 
correlated, the author proposes an acceptance scheme 
wherein one examines adjacent samples before sentencing 
a “lot”. A stochastic process is set up representing 
the system and Bayes theorem is applied to obtain a 
sentencing rule. 

The model is as follows: consider batches ..., B,_,, 
B,, By+1, ---» with counted defectives in random sample 
eases gs isos, Wet 6 (X,) =I, The assump- 
tion that m, either has the value a or b (a<b) is made, 
and if m, = a the batch is good; if m, = b the batch 
is bad. Runs of good and bad batches are distributed 
independently with mean run length 1/t, and 1/t,. It 
follows that the prior odds that a batch is bad is ii fires 
if t,+t, = 1 then the sequence is truly random. Con- 
sidering only B,_, and B,; 


Pr [imi Sh a] > ty ¢l x Cr i). 
Pr [m,-1 tg b| =) tq (1—1)/G,+%); 


Pr[m,-; =a, m, =b] =Pr[m,-, =}, m, =a] 
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If before sentencing B, we consider (X,_,, Xy-1441> 
..., X,,) the author’s rule is to accept B, if, in view of 
the above probabilities and the assumed prior distribu- 
tion Pr [m, = a] is sufficiently large. If W, and W, 
are respective losses from accepting a bad batch and 
rejecting a good, then we accept if Pr [m, = a,]<W,/W,. 
Each scheme is then characterised by a rejection rule, 
an expected loss, and by probability of an incorrect 
decision given a prior distribution. 


(E. H. Lehman, Jr.) 


8.8 (8.7) 


Some statistical aspects of the economics of analytical testing—Jn English 
Technometrics (1959) 1, 49-61 (6 references, 3 figures) 


General methods are developed for determining the 
number of analytical tests, n, to be made on a batch of 
chemical products in order to ‘minimise the total cost, 
testing cost plus the cost due to losses by incorrect 
decisions on batches. The author extends and applies 
the methods of others to routine laboratory testing of 
batch products. The assumption of a linear increasing 
testing cost and quadratically decreasing losses is made 
throughout. 

Two situations are explained: firstly, when the 
quality of a product is measured on a continuous scale 
the problem is to find a point on this scale, actually on 
the scale of averages based on n repeat tests, where the 
risk of passing B grade (bad) material is fixed and then 
the total cost is minimised. If T is the total cost per 
batch, the author minimises T(= Pc+nb) subject to 
the restriction that p(—a; X, n) = %, where P is the 
proportion of A grade (good) material rejected by the 
plan, c is the cost per batch for reprocessing, b is the 
cost per test of testing, X is the critical point beyond 
which a batch is rejected, and —a is the quality which 
should be rejected. An example is given where the 
distribution of quality and errors are both assumed to 
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be Gaussian: the proportion of “good” batches 
rejected is determined from bivariate normal tables, as 
quality and errors are both assumed normal. . 

The second situation discussed is the case where the 
amount of intermediate product to be added “in 
process ’’ for maintaining a given final quality is tested. 
The same problem of optimisation is posed; but the 
loss curve is assumed to be a quadratic function of x, 
where x is the per cent. error from the optimum amount 
that should be added: or, L = ax”. In the appendix 
it is proved that the average loss due to errors in the 
amount added = aa/n, where ao is the standard error of 
testing. Thus the total cost per batch is: T = (aa6/n)+nc, 
here c is the cost per test. Minimising, we find 


(Me Gipnf (a/c). If prior information is available on 
batch to batch variation (¢,), he shows that 
Pd po 
060 
T=nct+a "+ 
do ar no, 


and for a minimum: 


a a 
n=do =a) rei MBe 
c ot 


which is a smaller number of tests than when we have 
no prior information. (C. R. Hicks) 
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GOLDSMITH, P. L. & WHITFIELD, H. (Imperial College, London) 


8.9 (5.7) 


Average run lengths in cumulative sum chart quality control schemes—IJn English 
Technometrics (1961) 3, 11-20 (5 references, 2 tables, 8 figures) 


This paper’s primary aim is to give curves for the average 
run length before action for cumulative sum chart 
schemes. The cumulative sum is the algebraic sum of 
the deviations from some desired process mean, which is 
plotted against time. A V-mask with angle 206 is 
placed with vertex d time-units to the right of the last 
point and with axis horizontal. If any points lie outside 
the V, action is taken. 

Two process levels 9 when a process is on target, and 
Hy = Mo+ko are to be specified, then the average run 
lengths are to be specified for each, perhaps Ly = 500 
(X’s), and some small amount L, for p,. 

The average run lengths were carried out by Monte 
Carlo simulation for k = 0, 0-5, 1, 2, 3 with d’s 1, 2, 
5 and 8 and suitable values of 0. A comparison is then 
made between the present cumulative sum procedures 
and Shewhart charts for individual X’s, using only 
limits at +3-096, without warning limits and runs. 
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Average run lengths for cumulative sum charts are 
considerably smaller for low k’s than for this particular 
X-chart usage. Empirical formulae are given for Ly 
and L, in terms of d and 6. Curves are also given for 
the average run length in various cases where there is 
serial correlation in the data. 


(I. W. Burr) 


8.8 (1.6) 


An asymptotic distribution for an occupancy problem with statistical 


applications—ZIn English 


Technometrics (1961) 3, 79-89 (1 reference, 3 tables, 1 figure) 


In this paper the authors consider a problem of N items 
of which D are defective. Many times the items are 
divided into k lots of equal size s, and it is desired that 
the number of defectives in each lot be less than a given 
number, x. It is shown that the number of lots of x 
- defectives or more, m (x), is asymptotically multivariate 
normal. Thus Pr {m(x)<mpo} can be determined. 

The asymptotic results are derived through an 
occupancy argument, and the adequacy of the asymptotic 
theory is discussed. Applications to acceptance sampling 
are indicated. 


(P. H. Randolph) 
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HALD, A. (University of Copenhagen) 


8.8 (2.4) 


The compound hypergeometric distribution and a system of single sampling 
inspection plans based on prior distributions and costs—JIn English 


Technometrics (1960) 2, 275-340 (24 references, 18 tables, 12 figures) 


This long and full paper is divided into fourteen sections. 
Section one is introductory. Sections two to five define 
and discuss conventional terms ‘‘ Operating Character- 
istic’ and ‘‘ Acceptable Quality Level’? and describe 
and criticise four Quality Control systems, namely the 
Statistical Research Group of Columbia, U.S. Military 
Standard 105, Philips Standard Sampling and Dodge- 
Romig. The author’s criticism on all these systems is 
the arbitrariness of certain assumptions, for instance, 
always choosing the producer’s risk to be some fixed 
percentage, ignoring cost of accepted defectives, or not 
precisely defining the prior distribution. 

Section six presents two models based on cost 
minimisation: the average cost per item submitted in 
model one consists of three terms: 

(1) sampling inspection cost, | 

(2) expected loss due to accepted defectives, 

(3) cost of a rejected lot times probability of rejection. 


This model does not consider manufacturing costs: in 
model two an adjustment is made to include manufactur- 
ing costs. 
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The compound hypergeometric distribution and a system of single sampling 


Section seven derives the mean, variance, and co- 
variance of the two variables x, the number of defectives 
in the sample, and y = X—x, the remaining number of 
defectives in the uninspected portion of the lot; given 
that x is hypergeometric with parameters N, the lot 
size, n the sample size, and X the number of defectives 
in the lot: and that X itself is subject to some prior 
density function fy (X). 

The formulae which are lengthy for o2, 05 and o,, 
are given, depending upon the original fy(X). This 
density function is called normal, hypernormal or sub- 
normal as dy = 0, >0, <0 in this formula: 


ox = Npg (1+6y), 6y=—1, where p = 1-g = & (X)/N. 


Section eight shows how if fy (X) is hypergeometric, 
binomial, Polya, or any weighted average of these 
where the weights do not depend upon WN and_X, then 
x has the same distribution as XY with n substituted for N. 
That is, fy (X) is “ reproducible’ by random sample 
selection. 

Section nine defines the cost K(n, c) of a single 


continued 


8.7 (2.4) 
continued 


inspection plans based on prior distributions and costs—IJn English 
Technometrics (1960) 2, 275-340 (24 references, 18 tables, 12 figures) 


sampling plan as a function of n the sample size, and 
c the acceptance number. 


K (n, c) =n (k,—k,)+ Nk, 
BUN =; 1) y ia (x | Xk, In (x) 


where k, = cost of inspection per item, k, = cost per 
item of rejected lots, p, (x | X) = binomial probability 
of exactly x defectives in a sample of size n given there 
are X defectives in lot of size N, and g,(x) = marginal 
density function of x = 


Gn (x) aa ey In (X) Pn (x | X). 


An optimum sampling plan chooses n and c to 
minimise K given the other parameters. This involves 
study of several secondarily defined parameters such 
as the average fraction defective in the non-inspected 
part of accepted lots 


ce 


P,(C) = ss Dixy Ga) | BL Go); 


x= 


and the quantity 


Y (n, c) = 5) In (x) [k,— Pn (x)]. 
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The optimum plan is shown to be the solution with 
respect to n and c of the two inequalities: 
Pn (C)Sk,SPp (e+ 1) 
and 
(N—n—1)[y a+1, c)—y (, c)]<k,—-k, +7 (1, ©) 
<(N-n+ Ny (n, Cy lee i, c)]. 

Exact solutions depend upon the prior distribution 
fy (X) of X the number of defectives in a lot of size N. 

The following three sections give optimal solution for 
K(n, c) in the cases where fy (X) is rectangular, Polya, 
and binomial respectively. The Polya density function 
is similar to the binomial except that p is not constant: 
instead p = p(V, uw, y) = (pt Vy)/[1+(V+ uy] where 
V = number of defectives at any time, » = number of 
effectives at that time, and y is a parameter which may 
have any real value. If it is zero, the’density function 
is binomial. 

Tables are given to show optimum n, c, and K (n, c) 
for various forms of the prior distribution and for 
various values of the parameters. The paper concludes 
with a summary and a bibliography. 


(E. H. Lehman, Jr.) ~ 


I ——————————eeee 


rd 
bya kay ches ny 


Ras 


Lae 


JAMBUNATHAN, M. V. (University of Mysore, India) 


8.3 (4.3) 


A note on the efficiency of double sampling for stratification—Jn English 


Sankhya (1960) 22, 367-370 (3 references) 


To estimate the total value of a character y in a popula- 
tion double sampling may be employed if an auxiliary 
character x, which can be more economically measured, 
is available. In this method a large sample, called the 
first-stage sample, serves as a basis for stratification of 
the population according to the value of x, and at the 
second stage, stratified sampling is adopted for measuring 
the main character y and hence estimating the population 
total of y. 

Neyman, who introduced this method [J. Amer. 
Statist. Ass. (1938) 33, 101-116] also considered the 
problem of optimum allocation of the sampling units 
to the large and small samples, under the assumption 
that the cost of a single observation on y is the same in 
all strata. 

_ In an earlier book [Some aspects of sampling (1953) 
Bangalore City: India Book Company] the author 
considered the more general case of the cost of single 
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measurement varying from stratum to stratum and 
obtained approximate expressions for optimum alloca- 
tion. In this note he obtains an expression for the 
variance of the estimate under that optimum allocation. 
He compares it with the variance in the case of a simple 
random sample and obtains a condition which causes 
the double-sampling method to be more efficient. 


(T. V. Hanumantha Rao) 


8.8 (1.8) 


Probability of misclassifications. I. Variables—In German 
Qualitatskontr. Operat. Res. (1960) 5, 109-112 (4 references, 5 figures) 


Given a variable x, a population, for example an 
inspection lot, may be classified into three classes of 
quality in the following way: the population shall be 
of type A if w<x,, Bif x,Sp<x,, Cif u2x,. If the 
sample result x +ts/Nn contains the boundary x, or x, 
to which class should the population be attached? 
Here n denotes the sample size, X the sample mean, s 
the sample-standard deviation, and ¢ is the factor 
according to the level of significance chosen. 
Theoretically, the problem may be solved in a very 
simple manner by increasing n to n; = [ts/(x-—x) 3 
i= 1, 2, respectively. Practically, this method is 
impossible if the unknown p lies near x; or x,. The 
alternative method of solution is to classify by the 
original sample and to tolerate a certain probability of 
misclassifications (Pmj,): often, a compromise between 
the theoretical and the practical way may be chosen. 
Under the assumption of normality the author 
calculates the probability of misclassification by a 
Bayesian argument but without stating this. A diagram 
shows this probability for different sample sizes and 
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position of the sample mean. The assumption of 
uniform distribution of the true means is used for 
calculating the mean p,,,, dependent on the fraction 
paar where x, and x, are boundaries of the 
classes of quality. With reference to the mean p,,;, the 
author discusses the possibilities of determining the 
necessary sample size for a given maximum Of Pynjs, 
and of improving a given classification scheme. He 
also refers to the consequences which may result from 
classifying in a lower category in order to be on the 
safe side. 


(R. K. Bauer) 
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OLDS, E. G. 
Power characteristics of control charts—In English 


(Carnegie Institute of Technology, Pittsburgh) 


8.9 (5.1) 


Annual Cony. Trans., Amer. Soc. Qual. Contr. (1960) 195-208 (20 references, 6 tables) 


This article is a summary of several papers on power 
characteristics of control charts. Numerous statements 
are abstracted from these papers. Segments of some 
tables from these papers are presented to illustrate the 
effects on power of the number of samples, sample 
size, non-centrality, and width of the control limits. 
The author considers the case where a_ previous 
standard is not given. Hence, it is necessary to collect 
several samples, estimate the mean and standard devia- 
tion, and compute the control limits. The process is 
then defined as being in a state of control if all of the 
sample points fall within the control limits. A process 
is in a state of control during a time interval if all of 
the x’s are independent and f(x) is identical for every x. 

For counts, the control limits are ¢ +a,/ © where ¢ 
is the average number of counts from k samples. The 
probability of rejecting the null-hypothesis of a state of 
control when in fact the process is in a state of control 
is w. In general, exact computations of a are lengthy. 
An upper bound, «,, for « was found such that «,—>« 
as k->oo. Two tables of «, are presented which show 
the effects of k, a, and the true average number of 
counts. 
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Cumulative sum charts—Jn English 
Technometrics (1961) 3, 1-9 (21 references) 


This paper traces the development of process inspection 
schemes from Shewhart’s original work to the present 
techniques using cumulative sums. This includes 
control limits, warning limits, analysis of errors of both 
kinds and attendant costs, the study of average length 
of run before action is taken, and the extension suggested 
by sequential analysis to cumulative sums. The variable 
plotted on a cumulative sum chart is the algebraic sum 
of the deviations from some desired or typical process 
mean. Then, subjectively, if the points show a tendency 
down or up, that is, a change of slope, a process change 
is indicated. 
Objective criteria are briefly discussed : 


(i) take action whenever the points climb to a'height 
of h above the lowest point or conversely fall 
to h below the minimum, and 

(ii) placing a V-shaped mask with angle 20, with 
vertex at d time units ahead of the last point 
and axis horizontal, opening to the left, then 
take action if any points lie outside the V. 
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A single slippage in a process occurs when during 
all k sub-intervals the true means are equal to c; except 
for a single sub-interval with c,. The probability that 
at least the sample point in the slippage sub-interval 
is out of control (probability of detecting the single 
slippage) was studied in relation to k, a, c,, and c}. 
A few values are given for the power of detecting a 
double slippage, two intervals with true means of c}. 

For x-charts, the components of variance model, 
Xij = Hite;;, was considered where pu; is normally 
and independently distributed (0, 0707), ¢,; is similarly 
distributed (0, o7), and ¢;; and p; are independent. 
Control limits are determined by X +aA,r/3 where x is 
the overall mean and 7 is the average range or by 
X+as/,/n where s is the standard deviation and n is 
the number of samples per sub-interval. Type I and 
II errors were investigated for the null hypothesis that 

= 0 against the alternative that 0>0 for various 
values of k, n, 0, and a. The analysis of variance test 
is more powerful than the control chart test for these 
alternatives. The power of the control chart test was 
compared with the analysis of variance test for other 


alternatives. 
(D. W. Gaylor) 


8.9 (5.7) 


It is suggested that criteria for determining h or 0 
and d involve average run length before action is 
indicated, with long runs if process is satisfactory and 
short runs if “* off” by a specified amount. Comparisons 
of average run length with Shewhart charts are made 
for charts for individual X’s. The approach requires 
knowledge of process o and stability of o. Further 
problems needing solution are suggested. 


(I. W. Burr) 
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PLESZCZYNSKA, Elzbieta (Inst. Math., Polish Academy of Sciences, Warsaw) 


8.8 (6.2) 


Screening in statistical testing for correlated characteristics—In Polish 


Zast. Mat. (1960) 5, 47-59 (2 tables, 8 figures) 


In this paper the authoress considers lots of items whose 
quality depends decisively on a characteristic Y, which 
is difficult to test. It is known, however, that Y is 
correlated with another characteristic ¥, not so im- 
portant but easy to test. A further assumption is made 
that the pair of characteristics ¥, Y has two-dimensional 
normal distribution with known correlation coefficient 
p and known standard deviations, but with unknown 
means. To be concrete in subsequent description let 
it be assumed moreover that p<0O and that “ good” 
items are those with Y less than a given number yo. 

In reference to an earlier discussion by Oderfeld 
[“* Statistical testing for correlated characteristics ”’ 
Zastosowania Mat. (1959) 4, 255-264: abstracted in 
this journal No. 463, 8.8], the authoress investigates 
the statistical effects of screening, that is to say the effect 
of rejecting items with X smaller than a given x) and 
computes the dependence between per cent. defective 
in a lot before and after screening. 
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Two acceptance procedures involving screening are 
discussed, the one, where the lot is rejected or screened 
before acceptance according to the result of estimating 
mean value of Y on the base of a small pilot sample, 
and the second, where the parameters of screening 
depend on the result of that pilot sample. 


(S. Zubrzycki) 


8.6 (4.1) 


On the relative efficiencies of BAN estimates based on doubly truncated and 


censored samples—ZJn English 


Proc. Nat. Inst. Sci. India (1958) 24, 366-376 (6 references) 


The author supposes that in drawing samples from a 
statistical population distributed according to a prob- 
ability law f(x; 0), observations which are less than x9 
or greater than x, are not measured; but, the number 
of observations of each kind is recorded. Such a 
sample, we should expect, would supply less information 
regarding the parameter 0 than a sample where all 
individuals in the sample are measured. The author 
establishes a certain inequality and states that it proves 
the result mentioned above. Cases where the number 
of observations with values less than x or greater than 
x, are not separately recorded are also considered. 
The results are extended to cover distribution laws 
with more than one parameter. 


(S. John) 
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8.4 (4.3) 


On the estimation of the population size by capture-recapture method—ZJn English 
Bull. Calcutta Statist. Ass. (1960) 9, 93-110 (10 references) 


In this paper some methods of estimating population 
size by “‘ capture-recapture ’’ method are discussed. A 
sample of size n, is drawn out of an unknown N and 
the units are “released” back into the population 
with a suitable marking made on each of them. Next 
a sample of size n is drawn and the number r of them 
bearing the marks previously made is noted. The 
problem posed by the author is to estimate N. 

He considers two methods of drawing the second 
sample of size n: firstly by direct sampling, drawing 
a fixed number » of units; and secondly by inverse 
sampling, drawing units until a fixed number 7 of 
previously marked units appear. The cases where 
sampling is with replacement and also without replace- 
ment are both discussed. Since the maximum-likelihood 
estimators do not possess moments of any positive order 
in the case of direct sampling, whether it be with or 
without replacement, modified estimators due to Bailey 
[Biometrika (1951) 38, 293-306] are considered. It is 
shown that Bailey’s estimators have (1) bias which 
~ decreases more rapidly than that of the maximum- 
likelihood estimator, (2) finite unconditional moments 
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SEVERO, N. C. & NEWMAN, A. E. (Statist. Eng. Lab., Nat. Bur. Stand., Washington) 


of all orders, and (3) for sampling without replacement, 
conditional expected squared loss is not greater than 
that of the maximum likelihood estimator; the condition 
being that r is greater than a fixed positive number c. 
On the other hand, for sampling with replacement, the 
variance of Bailey’s estimator is larger than that of 
the maximum-likelihood estimator, though both tend 
to the information limit as the sample size increases. 

In the case of inverse sampling, the maximum-likeli- 
hood estimator is the minimum-variance-unbiased 
estimator when sampling is with replacement. For 
sampling without replacement, the maximum-likelihood 
estimator is biased, while Bailey’s estimator is unbiased, 
and has a smaller variance than that of the above 
estimator; though the two variances are equal, asymp- 
totically. Next, sequential procedures are discussed. 
A comparison of the information, in Fisher’s sense, 
for the sequential case with that of nonsequential case 


- shows the sequential procedure to be superior. 


(T. V. Hanumantha Rao) 


8.1 (4.9) 


A statistical chain-ratio method for estimating relative volumes of mail to given 


destinations—IJn English 


J. Res., Nat. Bur. Stand. C (1960) 64, 37-47 (9 references, 4 tables, 4 figures) 


A sampling method, called the “‘ chain ratio” method, is 
applied in estimating the distribution of mail by destina- 
tion. The details and results of the application of this 
sampling method to outgoing first-class letter mail in a 
study conducted in San Francisco are given: Los Angeles 
and Baltimore are other locations studied. 

The basic idea involved in the estimation formulae 
consists of multiplying together a chain of ratios. Thus, 
if the records at the post offices are available, or sampling 
plans are devised, to give the following unbiased estimates 


(i) the ratio of primary (the first stage of distribution 
of outgoing mail) to the total volume (7,/T), and 
(ii) the ratio of mail to the primary destination to 
the primary volume, (D,/7,), 
then the product of these two ratios gives the estimate 
of the ratio of mail to a primary destination to the total 
volume, (D,/T), or 


(D,/T) = (Dp|Tp) T/T). 
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This concept, with modifications, is extended to obtain 
estimates of the ratios to the final destinations through 


the primary, secondary, and tertiary stages of 
distributions. 
Variances and coefficients of variations for the 


estimators are given. 
also described. 


The method of collecting data is 


(H. H. Ku) 
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STRAUCH, H. (Niirnberg) 


8.9 (8.8) 


Continuous sampling inspection from the belt—Jn German 
Qualitatskontr. Operat. Res. (1961) 6, 21-25 (3 references, 3 tables, 9 figures) 


The problem of a combined sample and systematic 
inspection technique by attributes, suitable for belts or 
other continuous productions, has been solved theoretic- 
ally by Dodge [Ann. Math. Statist. (1943) 14, 264-279]. 
The defective units on the belt are assumed to be 
randomly distributed. To start the inspection, piece 
by piece is to be tested until i successive units are found 
to be acceptable; type I part of the series. From now 
on only every 1/fth unit is to be tested up to the next 
defective one; type II part of the series. Then another 
type I part of the series starts and so on. All detected 
defective units are replaced by non-defective ones. 

Let u be the mean length of an equal number of 
units of type I sections and let v be the mean length 
of type II sections. Then the mean proportion of all 
units to be tested is F = (u+fv)/(u+v). Furthermore, 
let p be the proportion of defectives in the production. 
After the inspection, the mean proportion of defectives 
is given by D=p(l—F). The author calculates 
u = (1—q')/pq' and v = 1/fp withg = 1—p. D assumes 
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SUKHATME, B. V. & KOSHAL, R. S. (Inst. Agric. Res. Statist., Delhi & FAO, Cairo) 


A contribution to double sampling—ZJn English 


a maximum D,,,,. A special sample design will be 
obtained by the following procedure. A certain value 
Of Dyyax iS preassigned. Then D,,,, is assumed for 
P (Dax) = (1 +iDpax)/(i+ 1) and it holds 
Fog! (Dee) oes ED pated (Daa)' ays 

By these formulae numerous pairs of values i and f 
are obtained which solve the problem. A graph given 
by Dodge shortens this procedure. 

Another numerical example shows that u is decreasing 
with decreasing i and increasing f. Therefore an 
approximate value for p, and the efficiency of the 
inspection team needs to be estimated before a maximum 
i and a minimum f will be chosen. This and some other 
results in sample design follow from a broad discussion 
of the formulae given. Finally the author refers to 
some generalisations of the inspection technique pro- 
posed for cases of more than one kind of defect and of 
grouped lots. 


(R. K. Bauer) 


8.3 (4.1) 


J. Indian Soc. Agric. Statist. (1959) 11, 128-144 (2 references) 


It is well known that the efficiency of an estimate of the 
population total of a characteristic Y can be increased, 
for a given cost, if information regarding an auxiliary 
variable XY is available or can be ascertained more 
economically than Y. One such technique is double 
sampling, where values of X are obtained for each 
individual in the sample, but values of Y are obtained 
only for individuals in a sub-sample. From the sub- 
sample, the regression of Y on X or the ratio of total Y 
to total X is first estimated. The estimate of total Y from 
the whole sample is obtained and used in conjunction 
with the regression or ratio based on the subsample 
to estimate total Y. 

Cochran [Sampling Techniques (1953) New York: 
Wiley] gave the formulae for the estimator, its variance 
and an estimator of that variance for the case of single- 
stage double sampling. This paper extends these results 
to the case of multi-stage designs and then to stratified 
multi-stage designs. As in the case of single-stage 
sampling, the estimators given are biased. A direct 
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analogue of the conditions under which the auxiliary 
variable X increases the precision of the estimates is 
derived. 

A numerical illustration is given. The authors state 
that they have extended these results to the case when 
the selection probabilities are unequal, and propose 
to publish them in a subsequent paper. 


(T. V. Hanumantha Rao) 
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TIKKIWAL, B. D. (Karnatak University, Dharwar, India) 


8.3 (4.1) 


An examination of the effect of matched sampling on the efficiency of estimators in the theory of 


successive sampling—Jn English 


J. Indian Soc. Agric. Statist. (1958) 10, 16-21 (8 references, | table) 


It was shown by the author in an earlier paper that when 
one or more characters of a population are studied on 
Successive occasions by a scheme of matched sampling, 
the efficiency of estimators of the time-conditioned 
parameters on any occasion can be increased by a scheme 
of matched sampling. The author shows in this paper 
that this efficiency goes on increasing with the occasion 
till it reaches a limiting value. 

The proof is given in two steps. Firstly, it is noted 
that the problem imposes certain restrictions on the 
sample size, etc., which reduces the multivariate case 
to the univariate one. Secondly it is observed that the 
variance of an estimator on any occasion is a certain 
multiple of the weight associated with the estimator on 
that occasion. It is shown that the weights, which are 
non-increasing with the occasion, tend to a limit as 
occasions increase. The author further examines, 
theoretically as well as by tabulation, the values of the 
weights for different values of the correlation between 
consecutive occasions and of the replacement fraction; 
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WETHERILL, G. B. (Birkbeck College, London) 


also the rapidity with which the weights converge to a 
limiting value. It is noted that the rate of convergence 
decreases with the increase in absolute value of the 
correlation between two consecutive occasions. While 
proving the non-increasing behaviour of the weights 
in the second step, the author obtains a lower limit to. 
the weights which provides an easy computational check 
on the weights. 


(S. G. Prabhu Ajgaonkar) 


8.8 (1.8) 


Some remarks on the Bayesian solution of the single sample inspection plan—ZJn English 
Technometrics (1960) 2, 341-352 (2 references, 4 tables, 1 figure) 


A prior distribution of the number X of defectives 
among K batches of size N; is assumed binomial 
_b(X; Nj, p) where p has the value p; with probability 
a;3 > @;=1, py>p2--->P,. Let W,; be the cost of 
i=1 

accepting a batch when p; is true and W,; be the cost of 
rejecting a batch when p; is true. Let po be the value of 
p at the breakeven point when W,; = W,; = 0: then 
W,; =0 if p;SPo and W,; = 0 if p;2Po- Equations 
for an optimum solution take a simple form if we assume 
pilq; = (p’)"* where p’ and « are suitable positive 
constants. 

To optimise n and c the “ neutral line ” is first deter- 
mined; that is, the locus of points (n, c) such that 
W,1; = W2;. This line has the equation i 


k 


& a:(”) piat* Was Mi) =0. 


i=1 
Next we let c be a function of n, determined by the 
neutral line and proceed to choose an optimal n. Let 
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the expected loss for a sample of size n be R(n). Then 
we seek an n such that AR = R(n+1)—R(n) changes 
from negative to positive at the optimum z. It is 


shown that 
k 


AR =} 2, a; (") Pidi © (W2:-W,:) (2p.-1) +1. 


rh 

In order to find the optimal m one must assume the 
as and W’s known and then perform repeated calcula- 
tions of several fairly lengthy equations, using several 
reasonably guessed values of n. It is shown that the 
optimum is “flat”? and hence quite large changes in 
some of the parameters have little effect on R(n). 
That is, the test is robust. The most critical parameter 
iS Po, the breakeven quality. Some numerical examples 
are given. 


(E. H. Lehman, Jr.) 
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ZUBRZYCKA, Ludmilla (Inst. Math., Polish Academy of Sciences, Wroctaw) 


8.0 (1.4) 


On the distribution of sampling points in a plane—IJn Polish 
Zastosowania Mat. (1960) 5, 161-171 (3 tables, 3 figures) 


If we want to estimate the integral of a certain function 
J (p) over a given region G by observing the values of 
J (p) in n points selected from G, then it is advisable to 
distribute the n points as uniformly as possible over 
G; for example by choosing them as vertices of a 
triangular net. The uniformity of such a sample of 
points is then disturbed in a crude manner if we add 
a few new points to our sample. It is not so if we 
select the points at random; but then we must expect 
random clustering of points and this is not desirable. 

Mrs Zubrzycka tries to find a compromise between 
the competitive requirements of uniformity of a sample 
for a given n and the preservation of this uniformity if 
a sample is enlarged point by point in that she defines 
a sequence of points p, in the unit-square such that 
each segment of this sequence is distributed more 
uniformly than, on the average, the same number of 
random points. As it is known, the sequence of rests 
modulo 1 of nz, denoted below by ¢,, where 


z = (\/5—1)/2 = 0-61803... 
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has remarkably good equi-partition properties over the 
unit interval O<r<1. 

Advantage is taken of this by defining p, as the image 
of ft, by a continuous Peano curve transforming a 
closed unit-interval into a closed unit-square. In 
table I the co-ordinates of points p,, n = 1, 2, ..., 248, 
are given. These points are also shown in a figure. 

The uniformity of the distribution of the points p, 
has been examined by comparing the variance of the 
observed numbers of points p, in equal disjoint squares, 
into which the unit square has been divided, to their 
mean. It has also been examined with the aid of a 
dendritic method so that the fractions of vertices of 
particular orders have been compared with the corres- 
ponding fractions for random points. 

The note terminates with some indications as to the 
practical use of the points p,,. 


(S. Zubrzycki) 


BOX, G. E. P. & HUNTER, J. S. (Statist. Tech. Research Group, Princeton Univ.) 


9.3 (11.4) 


Condensed calculations for evolutionary operation programmes—In English 
Technometrics (1959) 1, 77-95 (6 references, 1 table, 14 figures) 


This paper presents a detailed description of a simple 
method for calculating estimates of the “ effects’ and 
the error standard deviation for a 2” factorial design 
with centre points. 

Since an evolutionary operation programme is 
normally run by plant personnel and the calculations 
done by them, it is desirable to have a straight-forward 
method which will tend to reduce arithmetical errors. 

Examples of an information board and work sheets 
are given along with detailed numerical examples of 
their use for 2? and 2? designs. The 2° design is blocked 
and the calculation procedure is such as to remove the 
block effect. The range is used as an estimate of the 
standard deviation and table of correction factors is 
included. 

There is a discussion of the rationale of the com- 
putation method presented. 


(P. J. Cislak) 
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9.1 (9.2) 
Blacksburg, Va.) 


Intra- and inter-block analysis for factorials in incomplete block designs—In English 


Biometrics (1960) 16, 566-581 (12 references) 


The authors’ main objective in this paper is to extend 
the recent work of Kramer & Bradley [Ann. Math. 
Statist. (1957) 28, 349-361 and Biometrics (1957) 13, 
197-224] to the use of factorials in the several suitable 
classes of two-associate class, partially balanced in- 
complete block designs. Their secondary objective is 
to complete the earlier work through utilisation of 
inter-block information in combined intra- and inter- 
block estimators. They consider specifically the recovery 
of inter-block information in group divisible designs, 
both intra- and combined intra- and inter-block analyses 
for latin-square, sub-type L, designs and  latin- 
square, sub-type L, designs and then discuss other 
partially balanced incomplete block designs with two 
associate classes. In each situation tests of significance 
and estimators of factorial effects are obtained along 
with efficiencies of contrasts. 

Only basic results are given in this paper. The 
authors work essentially from treatment estimators for 
the basic designs and show how they yield the necessary 
information for the analyses for factorials. Usual 
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concepts of experimental design and factorials apply, 


and it is emphasised that attention in application must © 


be: given to choice and blocking of experimental units, 
applicability of models, estimation of weights in the 


recovery of inter-block information, and presentation — 


of results including tables and figures supplementary 
to basic analysis of variance tables. Methods of analysis 
follow usual patterns and numericai examples are not 
included. 

The material in this paper greatly increases the 
numbers of incomplete block designs available for use 
with factorial treatment combinations. Results for 
balanced incomplete block designs are well known and 
may also be obtained by setting 1, = 4, in formulas 
presented here for group divisible designs. 


(R. E. Walpole) 
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BROOKS, S. H. & MICKEY, M. R. (CEIR, Inc., Los Angeles, Cal.) 
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9.3 (4.7) 


Optimum estimation of gradient direction in steepest ascent experiments—IJn English 


Biometrics (1961) 17, 48-56 (4 references, 2 tables) 


Suppose R is a response surface functionally dependent 
on ” concomitant variables in some region ®. Let 
o be any subregion of ® such that in ¢, R is remote 
_ from its maximum. The method of steepest ascent, 

proposed by Box & Wilson [J. R. Statist. Soc. B (1951) 
13, 1-45] involves the step-wise progression from an 
initial set of trials in ¢ to a region $9 of high response. 
The path of ascent is determined by successive inferences 
as to the maximum gradient direction based on sets of 
trials conducted sequentially en route. 

It is desirable to assign ¢ trials per set in such a way as 
to maximise the gain towards ¢, per unit effort. In 
the absence of experimental error, ¢ = n+1 trials per 
set constitutes the optimal allocation: however, if 
experimental error is present, then the estimate of 
gradient direction becomes uncertain and the problem 
of optimum allocation of trials reappears. 
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The authors present a proof which indicates that the 
maximum gain toward $, per unit effort is achieved 
when ¢ = n+1 trials per set; regardless of the magnitude 
of either the experimental variation or the gradient in 
b. Proof stems from a derivation of & (cos 0) where 0 
is the angle of discrepancy between the true gradient 
and an estimated gradient resulting from a set of f¢ 
trials. Defining (1/t) & (cos) to be the expected 
improvement in location per unit effort, then (1/t) & (cos@) 
is shown to be maximised for t = n+1. 

Tabular values are given both for & (cos @) and for 
the decrease in efficiency associated with t>n+1 trials 
per set. 


(J. E. Dunn) 


9.1 (0.6) 


Proper spaces related to triangular partially balanced incomplete block designs—Jn English 
Ann. Math. Statist. (1960) 31, 498-501 (3 references) 


The incidence matrix for an incomplete block design, 
N, is the matrix with v rows, v is the number of treat- 
ments, and 5 columns, b being the number of blocks, 
where the typical element 7;;, is unity if the ith treatment 
occurs in the jth block, and is otherwise zero. The non- 
“negative symmetric matrix Q@ = NN’ of order v has 
elements g;,, where q;; is equal to the number of 
replicates of the ith treatment, and q;;, i #j, is equal 
to the number of blocks in which the ith and the jth 
treatment occur together. 

Knowledge of proper values and spaces of @Q is of 
interest in finding necessary and sufficient conditions 
for the existence of designs with given sets of parameters. 
Knowledge of the proper values of Q for several cases, 
including the triangular designs as given by Connor & 
Clatworthy [Ann. Math. Statist. (1954) 25, 100-112], 
has been utilised in the derivation of necessary conditions 
for the parameters of such designs. 


2/724 


In this note the author explicitly exhibits the proper | 


spaces of Q for triangular designs. Such information 
led Ogawa [Ann. Math. Statist. (1959) 30, 1063-1071: 
abstracted in this journal No. 1/701, 9.1] to derive other 
conditions for the existence of such designs. 


(B. S. Pasternack) 
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DANIEL, C. (116 Pinehurst Ave., New York 33, New York) 


Parallel fractional replicates—Jn English 
Technometrics (1960) 2, 263-268 (2 tables) 


This paper discusses 2’~‘ fractional replicates designed 
to study the effects of factors upon two response variables, 
y, and y, for certain situations in which y, and y, are 
affected by different subsets of the factors. Suppose 
for example that in a 2* experiment factor A affects 
y, only, factors B and C affect both y, and y,, and 
factor D affects y, only. This situation is described 
by the author as “influence pattern”? A—BC—D or 
1—2—1. In general, the influence pattern p,—p,—p, 
indicates that p, factors affect y, only, p, factors affect 
both y, and y,, and p; factors affect y, only. For p,—p3, 
the influence pattern is called a ‘“‘ symmetric influence 
pattern ”’. 

Some 8-run and 16-run fractional replicate plans are 
discussed for a few symmetric influence patterns. Two 
tables of symmetric pattern plans are also presented, 
_ one table for 8-run plans and the other for 16-run plans. 
For each influence pattern tabled the following are 
given: 

(i) the generators of the alias subgroups, plus other 
important subgroup members in parentheses, 

(ii) the generators of the principal blocks. 
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9.2 (4.8) 


All of the plans tabled have the following properties: 


(1) All confoundings of main effects with two-factor 
interactions are zero. 

(2) For 16-run plans, two-factor interactions are not 
confounded with each other. 

(3) For the 8-run plans one of the following conditions 
is satisfied: 
(a) two-factor interactions are not confounded. 
(b) two-factor interactions are confounded in 

pairs. 


(Helen Bozivich) 


9.3 (9.2) 


Third order rotatable designs in three dimensions—In English 
Ann. Math. Statist. (1960) 31, 865-874 (5 references, 3 tables) 


In this paper the author discusses the construction of 
third-order rotatable designs in three dimensions, that 
is to say, designs involving three factors. He presents 
the conditions for the existence of third order designs 
and the two conditions which are necessary for non- 
singularity. It is shown that a third-order rotatable 
design is singular if and only if all of its points, excluding 
centre points, lie on a k dimensional sphere centred at 
the origin. 

The four previously known third-order designs are 
listed in their original geometric framework and are 
also translated into the mathematical structure used 
by the author in constructing several infinite classes of 
third-order designs. Sequential designs, that is to say, 
those formed from two second-order rotatable designs 
in three dimensions, are the author’s primary con- 
sideration. An example for one infinite class is con- 
sidered in detail. It is shown that two of the previously 
known designs are members of this class and that their 
design parameters are located at the extremes possible 


for this class. 


The author states that he intends to present, in a 
future report, some specific single designs which are 
to be selected from the infinite classes developed in 
this paper and will be presented in a form in which 
they can be used conveniently by experimenters. See also 
abstract No. 2/727, 9.3. 


(J. W. Wilkinson) 
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DRAPER, N. R. (Math. Res. Centre, U.S. Army, Madison, Wisconsin) | 973972) 
A third order rotatable design in four dimensions—In English 
Ann. Math. Statist. (1960) 31, 875-877 (2 references, 1 table) 


In this paper the author presents a third-order rotatable 
design in four dimensions, that is to say, a design for 
four factors, which requires 96 points plus No centre 
_ points where m9 may be zero. 

This design is a combination of four second-order 
rotatable arrangements, each containing 24 points. | 
This permits the use of four convenient blocks of equal 
size or, alternatively, sequential performance in several | 
ways. The notation and definitions of a paper by the 
author in the same issue are used in obtaining the 
design; see abstract No. 2/726, 9.3. 

It is mentioned that previously only one third-order 
rotatable design in four dimensions was known and 
that design required 128 points plus centre points. 


(J. W. Wilkinson) 
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GRAYBILL, F. A. & SESHADRI, VY. (Oklahoma State University, Stillwater) 9.1 (7.3) 


On the unbiasedness of Yates’ method of estimation using interblock 
information—ZJn English 


Ann. Math. Statist. (1960) 31, 786-787 (2 references) 


This refers to the results in a previous article [Ann. 
Math. Statist. (1959) 30, 799-805: abstracted in this 
journal No. 1/474, 9.1]. A proof is given of the un- 
biasedness of Yates’ procedure of using the estimator 
‘based on recovery of interblock information if the 
estimate of the block variance component is positive, 
and the randomised complete blocks estimator if this 
estimate is negative. 


(R. L. Anderson) 
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9.1 (0.6) 


On the uniqueness of the triangular association scheme—ZJn English 
Ann. Math. Statist. (1960) 31, 492-497 (6 references) 


It has been shown by Connor [Ann. Math. Statist. 
(1958) 29, 262-266] that the relations that exist among 
the parameters of the triangular association scheme 
themselves imply the scheme if n=9. 

Shrikhande [Ann. Math. Statist. (1959) 30, 39-47: 
abstracted in this journal No. 1/292, 9.1] extended 
this argument by demonstrating that this result is also 
true for nS<6: the problem has no meaning for n<4. 

The aim of the author in this paper is to show that 
the result is valid ifm = 7, but that it is invalid ifn = 8. 
For the case m = 8, the author proves by the con- 
struction of a counter-example that the relations among 
the parameters of the triangular association scheme do 
not necessarily imply the triangular association scheme. 

The results of this paper have also been obtained 
through the use of different methods by Chang [Science 
Record (1959) New Series, 3, 604-613]. Chang has also 
shown there are exactly three counter-examples when 
n = 8, see [Science Record (1960) New Series 4, 12-18]. 


(B. S. Pasternack) 
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KISHEN, K. (Department of Agriculture, Lucknow) 


9.2 (9.1) 


A note on the construction of the (21°, 2'') and other associated confounded designs 
keeping up to second interactions confounded—ZJn English 


J. Ind. Soc. Agric. Statist. (1959) 11, 180-185 (3 references, 2 tables) 


Let m,(r, s) be the maximum value of m for which it 
is possible to have a symmetrical factorial design 
(s”, s") ; that is, a design with m factors each at s levels, 
accommodated in blocks of size s’, so that degrees of 
freedom belonging to a main effect or interaction 
involving t or less factors are not confounded. Bose 
[Sankyd (1947) 8, 107-166] and Rao [J. R. Statist. Soc. B 
(1947) 9, 128-140] showed that m3(r, 2) = per and 
from this it follows that in a factorial experiment in 
which each factor is at two levels, and block size is 2°, 
the maximum number of factors that can be accom- 
modated so that no main effect or first order interaction 
is confounded is 16. 

In this note the author, using Bose’s method, con- 
structs the (2'°, 21!) design, by selecting 16 points in 
PG (4, 2), no three of which are collinear. The intra 
block sub-group of the 32 combinations in the key 
block is given in a table. 


From the (2'°, 211) design, the author gets (2'*, 2*°) © 


(2'4, 2°), ... (27, 27), (2°, 2) designs by omitting columns 
successively. The numbers of degrees of freedom 
confounded for the various higher order interactions in 
these designs are given in a further table. 


(C. S. Ramakrishnan) 


MIKHAIL, W. F. (University of North Carolina, Chapel Hill) 


9.1 (-.-) 


An inequality for balanced incomplete block designs—In English 
Ann. Math. Statist. (1960) 31, 520-522 (2 references) 


In this paper the author discusses a balanced incomplete 
blocks design with v treatments in b blocks of size 
k<v, each treatment appearing in r blocks. It is proven 
, that if v = nk, n an integer, b>v+r—1. 

This is a weaker condition than the usual one that 
the b blocks can be divided into r sets of n blocks each, 
so that in each set every treatment occurs exactly once. 


(R. L. Anderson) 
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PAVATE, M. V. (Central Tobacco Res. Inst., Rajahmundry, India) 
Combined analysis of balanced incomplete block designs with some commori 


treatments—ZJn English 


Biometrics (1961) 17, 111-119 (3 references, 3 tables) 


The purpose of this paper is to suggest a simplified 
method to obtain adjusted treatment components for 
the combined analysis when the individual experiments 
are laid out in balanced incomplete block designs. 
Gomes & Guimares [Biometrics (1958) 14, 521-526: 
abstracted in this journal No. 1/113, 9.7] have considered 
the case when the individual experiments are laid out 
in randomised complete block designs. 

The author presents briefly the notation used. C 
denotes the number of treatments common to all 
experiments. The adjusted treatment sum-of-squares 
for the combined analysis is found by two methods. 
The general method suggested by Rao [J. Amer. Statist. 
Ass. (1947) 42, 541-561] is outlined. This method treats 
the combined data as that of a new incomplete block 
experiment. Hence the problem of combined analysis 
reduces to the solution of a set of adjusted normal 
equations. 

The simplified method, which the author introduces, 
makes use of the individual analysis already performed 


Pe? |732 a 


discussed. The combined analysis is by the author’s i 


9.1 (7.1) 


on the balanced incomplete block experiments. This 

method is presented in some detail. Estimation of — 

treatment estimates and variances of treatment differ- — 

ences is also discussed. Four special cases of combined 

analysis of these balanced incomplete block experiments : 

are discussed. 

(a) Experiments with only one common treatment _ 

which might be a control. In this case C = 1. _ 

(b) Several experiments with identical parameter — 

values and treatments. In this case C =v 

where v is the number of treatments in an — 

individual block. 

(c) & (d) Cases where balanced incomplete block designs _ 

reduce to randomised complete block designs — | 

with C common treatments. These cases are | 
considered in detail by Gomes & Guimares. 

Finally the author presents an example dealing with 

two experiments: the individual analysis of both is 


simplified method. 
(J. J. Bartko) 


RAGHAVARAO, D. (University of Bombay) 


9.1 (0.4) 


A generalisation of group divisible designs—In English 
Ann. Math. Statist. (1960) 31, 756-771 (9 references, 10 design plans) 


Group divisible designs were first introduced by Bose 
& Connor [Ann. Math. Statist. (1952) 23, 367-383], 
and were later extended to m-associate classes by 
P. M. Roy [Science and Culture (1953) 19, 210-211]. 
Roy called these designs “ Hierarchial Group Divisible 
designs with m-associate classes’. Further work in 
this direction was not pursued in the literature until 
the publication of this paper in which the author studies 
these designs systematically. 

In the second section the author gives a compact 
definition of a group divisible m-associate design in 
which the parameters of the design are presented in a 
slightly different form than that of Roy. In Section three 
the author demonstrates the uniqueness of the associa- 
tion scheme from the parameters of the design. The 
designs are divided into (m+1) classes in Section four, 
and in Section five some interesting combinatorial 
properties of these designs are discussed. The following 
section offers a presentation of some known results 
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RAGHAVARAO, D. (University of Bombay) 
Some aspects of weighing designs—ZJn English 


concerning the Legendre symbol, the Hilbert norm 
residue symbol and the Hasse-Minkowski invariant. 
The final section of this paper describes necessary con- 
ditions for the existence of symmetrical regular group 
divisible m-associate designs, and numerical illustra- 
tions of these designs are given in the Appendix. 


(B. S. Pasternack) 


9.8 (9.1) 


Ann. Math. Statist. (1960) 31, 878-884 (9 references) 


The author supposes we have WN objects to be weighed 
in N weighings with a balance having no bias. Let the 
design matrix X (Nx N) = (x;,;) where x;; = 1(—1) if 
the jth object is placed in the left or alternatively the 
right pan in the ith weighing, and x;; = 0 if the jth 
object is not weighed in the ith weighing. In a previous 
paper by the author “‘ Some optimum weighing designs ”’ 
[Ann. Math. Statist. (1959) 30, 295-303: abstracted in 
this journal No. 1/288, 9.8] it was shown that the 
design matrices X = Py and X = Sy are the most 
efficient weighing designs under Kishen’s definition of 
efficiency, see “On the design of experiments for 
weighing and making other types of measurements ig 
[Ann. Math. Statist. (1945) 16, 294-300]. This is related 
to the magnitude of the trace of the matrix Core 
when N is odd and N = 2 (mod 4), respectively, subject 
to the conditions that the variances of the estimated 
weights are equal and the estimated weights are equally 
correlated. Here PyPy = (N—1)IytEyy and SySy = 
(N—1) Iy where Iy is the identity matrix of order N 
and Eyy is an Nth order matrix with positive unit 
elements everywhere. = 

In this paper it is shown, subject to the same conditions 
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as stated above, that the P, matrices are the best 
weighing designs under the definitions of Mood, see 
“On Hotelling’s weighing problem” [Ann. Math. 
Statist. (1946) 17, 432-446] and Ehrenfeld, “On the 
efficiencies of experimental design ” [Ann. Math. Statist. 
(1955) 26, 247-255], when N is odd; while the Sy matrices 
are the best weighing designs under the definition of — 
Ehrenfeld when N=2(mod 4). Under Mood’s — 
definition of efficiency, the best weighing design X, 
when N = 2(mod 4), is shown to be that for which — 
X'X = (N—2)Iyt+2Eyx. The best weighing design X _ 
according to Mood’s definition is the one for which the | 
determinant of X’X is maximum. The best weighing 
design according to Ehrenfeld’s definition is the one for — 
which the smallest characteristic root of X’X is maximum. — 

It is shown that the necessary condition for the 
existence of a weighing design X such that X’X¥ = © 
(N—2) Iyt2Eyy is that N = {44+(3f7+4)73/3, where — 
f is an integer. For N<200, this implies the existence 
of such weighing designs only for N = 6 and N = 66. 
The author also shows the impossibility of S, matrices — 
of order N = 22, 34, 58 and 78. 


(J. W. Wilkinson) 


pe 


RAO, P. V. (College of Science, Nagpur, India) 9.1 (0.6) 
The dual of a balanced incomplete block design—In English 
Ann. Math. Statist. (1960) 31, 779-785 (9 references) 


Shrikhande [Biometrics (1952) 8, 66-71] and P. M. Roy 
[Sankhyad (1954) 14, 39-52] have shown that certain 
balanced incomplete block designs can be dualised 
to give partially balanced incomplete block designs 
with exactly two associate classes. J. Roy & Laha 
[Sankhyda (1956) 17, 115-132] have obtained a necessary 
and sufficient condition for the dual of a balanced 
incomplete block design to be a partially balanced 
incomplete block design with two associate classes. 

In this paper, a general result regarding the dual of a 
balanced incomplete block design is established, and 
the results of Shrikhande and P. M. Roy are obtained 
as particular cases. An illustration to show the use 
of the result when the dual is not a 2-associate partially 
balanced incomplete block design is also given. 


(P. V. Rao) 
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SCHWARZ, G. (Columbia University, New York) 9.2 (7.4) | 
A class of factorial designs with unequal cell-frequencies—Jn English ; 
Ann. Math. Statist. (1960) 31, 749-755 (8 references) 


A class of multifactorial designs are defined and analysed. 
The designs considered each have a total number of 
observations that cannot be divided equally among the 
cells of the designs; however, by distributing the 
observations in a manner that is in a certain sense 
symmetrical, the equations that determine the least- 
squares estimates of the linear parameters become 
explicitly solvable. 

The first case which the author treats in this paper 
is the case of two non-interacting factors with arbitrary 
numbers of levels. In the n-factor case, he restricts 
himself to the discussion of factors having an equal 
number of levels. After defining the designs, the estimates 
are computed. Some general discussions of the sym- 
metries and algebraic properties involved conclude the 


paper. 
(G. Schwarz) 
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9.1 (0.6) 


A matrix substitution method of constructing partially balanced designs—In English 


Ann. Math. Statist. (1960) 31, 34-42 (6 references) 


Given an incomplete blocks design for v treatments 
and 6 blocks with the following incidence matrix, 
N = (n;;), where n;; is the number of times the ith 
treatment occurs in the jth block. The author is con- 


cerned only with designs for which n;; = 0 or 1 and 
b v 

>) my=r<b and > nj =k<v. 

a4 . v= 1 

If each treatment appears with every other treatment 

in 4 blocks, that is to say )' n,n; =4 for 1 4 i, the 


design is balanced; if certain eamenn appear together 
in 4, blocks, others together in 1, blocks, ..., others in 
Am blocks, the design is partially balanced with m 
associate classes. 

Vartak [Ann. Math. Statist. (1955) 26, 420-438] 
considered the construction of experimental designs 
utilising kronecker products of incidence matrices. 
His method essentially consisted of the replacement of 
two elements, 0 and 1, by two matrices. Shah, in an 
earlier paper [Ann. Math. Statist. (1959) 30, 48-54: 
abstracted in this journal No. 1/291, 9.1] generalised 
this approach by using the incidence matrices of balanced 
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SHAH, B. V. (University of Bombay) 
Balanced factorial experiments—Jn English 


incomplete block designs for substitution. In the 
present paper the author continues his generalisation 
by considering the incidence matrices of partially 
balanced incomplete block designs and factorial experi- 
mental designs as replacement elements. 

The author introduces fundamental notions concern- 
ing canonical matrices with reference to the partially 
balanced incomplete block designs in Sections 2 and 3 
followed by a discussion of associable designs and their 
properties in Section 4. In Section 5 balanced matrices 
are defined and in Section 6 a method is presented for 
the construction of designs by substitution of the 
incidence matrices of a balanced matrix. Finally, the 
application of this method to the construction of factorial 
experimental designs is given in Section 7. 

The author indicates that these methods may be of 
considerable importance in the construction of con- 
founded asymmetrical factorial designs when there are 
many factors. 


(B. S. Pasternack) 


D240) 


Ann. Math. Statist. (1960) 31, 502-514 (15 references) 


Factorial experiments are usually of such a nature as 
to prohibit all possible treatment combinations from 
appearing in a block: it is necessary, therefore, to 
consider methods of confounding. The author in an 
earlier paper [Ann. Math. Statist. (1958) 29, 766-779] 
considered the problem of achieving ‘“‘ complete balance ” 
over various interactions in factorial experiments. 
“Complete balance” is achieved over an interaction, 
if and only if all the normalised contrasts belonging 
to the same interaction are estimated with the same 
variance. In the present paper, the author considers 
a class of factorial experiments which are defined as 
balanced factorial experiments. An experiment is 
called a balanced factorial experiment, if the following 
conditions are satisfied: 
(i) Each of the treatments is replicated the same 
number of times. 
(ii) Each of the blocks has the same number of plots. 
(iii) Estimates of contrasts belonging to different 
interactions are uncorrelated with each other. 
(iv) “ Compiete balance” is achieved over each of 
the interactions. 
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In the fifth section of this paper Shah proves several 
theorems which indicate a method of analysis for 
balanced factorial experiments which is particularly 
useful when estimates of interactions at different levels 
are required. A method for constructing these experi- 
ments is' given in the final section. 


(B. S. Pasternack) 
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9.1 (0.6) 


Relations among the blocks of the kronecker product of designs—In English 
Ann. Math. Statist. (1960) 31, 772-778 (7 references) 


For certain incomplete block designs interesting relations 
among blocks have been found. Fisher [Ann. Eugenics 
(1940) 10, 52-75] has shown that in the case of a 
symmetrical balanced incomplete block design with 
parameters v = b, r = k, A, any two blocks have exactly 
4 treatments in common. For an affine ‘resolvable 
balanced incomplete block design with parameters 

v = nk =n? {(n—1) t+1}, b = mr = 2 {n?t4+nt+hi, 
4 = nt+1, Bose [Sankhyd (1942) 6, 105-110] has shown 
the blocks can be divided into sets of n blocks, such 
that each set is a complete replication and any two 
blocks have either k?/v = (nt—t+1) or nil treatments in 
common according to whether they belong to different 
groups or the same group. Similar results have been 
obtained by Connor [Amn. Math. Statist. (1952) 23, 
57-71] and Bose & Connor [Amn. Math. Statist. (1952) 
23, 367-383]. 

In this paper the author confines his attention to 
partially balanced incomplete block designs having two 
or three associate classes with the objective of obtaining 
analogous relationships among blocks of balanced 
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incomplete block designs, as those described above, 
for the blocks of their kronecker product. The inter- 
connection of the results he obtains and the theorems 
on inversion of designs, P. M. Roy [Bull. Cal. Math. 
Soc. (1954) 46, 47-58] are also discussed. 

In one main theorem, the author proves that if for 
each of two designs with incidence matrices N, and N>, 
there exists a pair of blocks with, respectively, M, and 
M, treatments in common, the design with incidence 
matrix based on the kronecker product N, x N, will 
have a pair of blocks with M,M, treatments in common. 

Other theorems and corollaries show how to classify 
the method of the classification of blocks according to 
the number of treatments that they have in common 
with some other block. 


(B. S. Pasternack) 


BALOGH, T. (Kossuth University, Debrecen) 


10.8 (0.6) 


Generalisation of a theorem of U. Grenander—In English 
Second Hungarian Mathematical Congress (1961) (4 references) 


The extension of the theory of stochastic processes to. 


the case of matrix-valued variables made it necessary 
to determine the minimum of the spur of the so-called 
generalised Toeplitz-form 


20 
=) Po (2) 6) PS (2) dxs7 2 =e 
T Jo 


under the condition P, («) = E, where P,, (z) is a trigono- 
metric polynomial matrix, and f (x) is a positive definite 
Hermitian functional matrix. 

In the case of ordinary Toeplitz-forms the above 
problem has been solved by Szegé. The method by 
which we obtain our result is essentially a generalisation 
to the case of matrices of the proof of Szegé. If we 
apply the theorem so obtained to matrix-valued 
stochastic processes, we obtain a natural generalisation 
of the theorem of Grenander. 


(T. Balogh) 
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BARLOW, R. E. & HUNTER, L. C. (Sylvania Electric Defense Lab., 
Mountain View, Cal.) 


System efficiency and reliability—In English 


10.1 (10.2) 


Technometrics (1960) 2, 43-54 (7 references, 5 figures) 


This paper presents a model useful in determining the. 


reliability and efficiency of large, complex systems. 
Repair is an integral part of the model and the assump- 
tion of component independence need not be made. 

A system operating in time can be in any one of many 
states. At a given time the system will be in one of 
the states. Thus the state of the system is essentially 
a random variable and as a function of time it is a 
stochastic process [X(t), OSt<oo]. Usually it is 
advantageous for the system to be in a certain class of 
states. The designation of these states is made by the 
definition of a gain function g defined on the state-space 
of the system. Since the state-space of the system as a 
function of a time is a stochastic process, the gain 
function is composed with the process; that is, g [X (¢)]. 
The reliability of the system at time ft, R(t), is defined 
as the expected value of the gain function. For example, 
if g[X(]=1 when the system is in a favourable 
state and g [X (t)] = 0 otherwise, the reliability is the 
probability that the system is in a favourable state at 
time f. 
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The final factor in the model deals with the possibility 
that the system is more likely to fail at one time than 
another. Thus the environment may provide a distribu- 
tion function F(t) on the time axis. Efficiency is defined 
as the expected value of the reliability function. If the 
system is operating under constant threat, efficiency is 
the time average of the reliability function. The limiting 
efficiency is defined in terms of the time average. 

The model presented is very general and permits the 
use of probabilistic techniques developed for stochastic 
processes. However, actual computations of system 
efficiency and reliability are easily accomplished only 
for systems with finite state spaces and for those which 
can be described by Markoff processes. 


(A. P. Berens) 
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BHAT, B. R. (University of Western Australia, Nedlands) 


10.2 (4.1) 


Maximum likelihood equation for positively regular Markoff chains—Jn English 


Sankhya (1960) 22, 339-344 (8 references) 


In this paper the author considers the problem of 
maximum likelihood estimation of an unknown para- 
meter entering in the specification of the transition 
probability matrix of a simple stationary and positively 
regular Markoff chain. The main results can be stated 
as follows: let 2(@) denote the transition probability 
matrix when the true value of the parameter is @ and 
p, the estimate of m based on the empirical distribution 
of n successive observations. Let T be any Fisher- 
consistent estimate of § based on p. The variance of 
T then satisfies the well-known Cramér-Rao information 
inequality under certain regularity conditions. If 6* 
denotes an estimator of 0 which makes the likelihood 
a maximum and 


|9-@0| >6 ah 

for every 6>0 and all i, 0) being the true value of the 
parameter, then 6* is a consistent estimate of 6. Under 
certain conditions (which are not stated in the paper) 


the uniqueness of the root of the likelihood equation, 
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BLOM, G. (University of Lund, Sweden) 


at which the maximum is attained, and the asymptotic 
normality of its distribution are asserted. All these 
results are generalisations of those of Rao [Sankhya 
(1957) 18, 139-148] in the case of independent . 
observations. 


(K. R. Parthasarathy) 


LOT) 


Hierarchical birth and death processes. I. Theory—Jn English 


Biometrika (1960) 47, 235-244 (9 references) 


A Markoff chain is defined with certain transition rates, 
m—1 indices for “‘ births ’’ and m indices for “‘ deaths ”’. 
Thus transitions occur only between states with one 
more or less indices, and only one index at a time is 
either deleted or added. Each index is restricted to 
the same finite set of values 1 to N. 

Steady state probabilities are found for the finite 
(mn) “last come, first served ” process in which only 
the index at the right-hand end can change. It is shown 
that the same solutions are obtained in the “ first come, 
- first served ” process in which additions to the indices 
are at the right-hand end and deletions at the left-hand 
end, provided that g, = q2 = .--=4n- 

A sufficient set of conditions is given for the general 
infinite processes to tend to a steady state. The steady 
state probabilities are shown to be the same as those 


= 
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derived earlier, for Type A processes in which 


Olea D000 me AY EY 
and for Type B processes in which a, _1,, = CmDm, y 
where ¢) > 0 and? v= 1)..5 40, 
In a second paper the author discusses some practical 
applications of this theory. See abstract No. 2/745. 
10.1. 


(A. S. Miller) 
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BLOM, G. (University of Lund, Sweden) 10.1 (-.-) 
Hierarchical birth and death processes. II. Applications—Jn English 
Biometrika (1960) 47, 245-251 (5 references) 


The theory of hierarchical birth-and-death processes is 
discussed in Part I of the paper. See abstract No. 2/744, 
10.1. 

_ It is shown in this second part of the paper that this 
theory may be applied to the machine interference 
problem with different rates of breakdown and different 
service times for different machines and with either one 
or several operators. 

The theory is also applied to a telephone traffic 
problem. A group of channels is used by N subscribers. 
Subscribers make calls at random but at different rates 
and the durations of their calls have (i) exponential 
distributions with varying parameters, (ii) mixed 
exponential distributions. If all m channels are occupied 
no queueing is allowed, the subscriber makes another 
attempt later, 


(A. S. Miller) 
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CIUCU, G. (University “ C. I. Parhon”’, Bucharest) 10.1 (1.5) 


I. Some properties of chains with complete connections 
II. The law of large numbers for chains with complete connections—IJn Rumanian 


Com. Acad. R.P.R. (1960) 10, 561-563, 643-645 (7 references) 


In these papers concerning chains with complete con- 

nections the author considers a chain with a denumerable 

set of states Q, defined by the transition probabilities 
Bocas) ce ||. Q% 2)=0Ge—N), 


jeEe—-N 
eee 10.0 cAV OY = il OY, 
j=l 

Qf = Q(1SjSs/). 
A chain of this kind is of type (G), if there exists a function 
Pp, (A), defined for all Ae OO ari acl Demet) 
tad two constants L, A>0 such that 

| BAC: A®)—P? (A) | <Lexp (—A,/n), 


where P‘”, (c; A) is the n-step transition probability. 
In the first note it is shown that for chains with 
complete connections of type (G) a central limit theorem 
and the theorem of the iterated logarithm hold true; 
in the second note the law of large numbers is considered. 


(R. Theodorescu) 
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COX, D. R. (Birkbeck College, London) 


10.1 (2.7) 


On the number of renewals in a random interval—Jn English 


Biometrika (1960) 47, 449-452 (7 references) 


A renewal process is defined by the author in this paper 
as a sequence (X,, X>, ...) of independent identically 
distributed random variables, taking only positive values, 
and is assumed here to have an absolutely continuous 
distribution. If X; is interpreted as the lifetime of the 
ith component to be used, renewals occur at times 
X,, X, +X, ..... The random variable N, is defined as 
the number of renewals (0, 7). 

There is an extensive literature both on the exact 
properties of N, and on asymptotic results which hold 


as too. The exact results involve the inversion of a 


Laplace transform, so that simple explicit results are 
obtained only for especially simple distributions of the 
ix. 

‘Suppose now that T is a positive random variable 
distributed independently of the Y;: 
number of renewals in (0, 7). It is pointed out that 
when T has a IT distribution, the properties of N can be 
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let N be the 


obtained without the inversion of a Laplace transform. 
Hence, particuarly when the index of the I distribution 
is a small integer, simple explicit results can be obtained. 
A number of special cases are worked out. 


(D. R. Cox) 


10.4 (0.1) 


Deterministic customer impatience in the queueing system. GI/M/I—IJn English 


Biometrika (1960) 47, 45-52 (5 references) 


This paper gives the solution of an integral equation 
occurring in the analysis of a particular queueing system 
when the common distribution function [V (x)] of the 
identically distributed random variables [V,,] takes the 


f Cor Nes 
ae V(x) = ve [ A(yje dy 


x 
where x>0. The properties of the queueing system are 
as follows: 

(i) Customers arrive at time instants (¢,) after fp = 0 
so that the t,=¢%,—-t,-1, (m=1, 2, ...) are 
independent’ and identically distributed non- 
negative random variables with common distribu- 
tion function A (x) and finite expectation 


[o8) 
Co | xdA (x). 
¥ . 

(ii) Customers are served in order of arrival and 
service is immediate if a customer arrives and 
finds the server free. 

(iii) A customer will wait for service only for a time 
not exceeding W at which time he departs never 
to return. ; 

(iv) The service times of customers who receive 
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x 


service are independently and identically dis- 
tributed random variables which are independent 
of {t,} and with a common distribution function 
Bia el ee torn). 
The case where {t,} is a Poisson process was studied 
by Barrer [Operat. Res. (1957) 5, 650-656] who obtained 
the limiting distribution of queue-size at an arbitrary 
instant of time. 
The second section of this present paper contains 
the solution of the integral equation 


G(x) = ie G(x—y) dV (y), 


and this is illustrated for the case where {t,} is an 
Erlang (£,) process with mean k/A. The following 
section relates the analysis to the case of W = oo and 
them titanic: untensity = 1-01 /(ro) <i. Where. 7 = "1p. a 
solution does not exist. The case of />1, W = o is 
the limiting-time distribution of a queue with Poisson 
input process with parameter m and service-time dis- 
tribution A(x). The last two sections of this paper deal 
with the limiting distribution of queue size at the 
instants a service commences and the case of the many 
server queue. (Florence N. David) 
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DANSKIN, J. M. (Rutgers University, New Jersey, USA) 


The theory of reconnaissance—In English 


10.5 (0.7) 


ARL Technical Report 60-337 (1960) Aero. Res. Labs., Wright-Patterson Air Force Base. 


iv+20 pp. (2 references) 


In a previous paper [WADC Technical Note 59-409: 
abstracted in this journal No. 2/381, 10.5] with the 
same title the author considered the one-sided recon- 
naissance problem, in which the reconnoiterer seeks to 
maximise the information or minimise the confusion 
obtained as a result of the expenditure of a given effort; 
it being assumed that the side being reconnoitered 
remains passive. This problem was formulated as a 
problem in information theory. This report concerns 
the two-sided reconnaissance problem, in which the 
side being reconnoitered seeks to minimise the informa- 
tion or maximise the confusion obtained by the recon- 
noiterer, while maintaining at least a certain minimum 
acceptable threat with a fixed budget. This problem is 
formulated as a zero-sum, two-person game: for example 
consider a defender trying to protect ballistic missile 
installations against reconnaissance. When planning 
at a period substantially preceding the period in which 
he will be reconnoitered, he is limited by a budget but 
not by equipment, and may choose to divert money 
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(Rutgers University, New Jersey, USA) 


to, for instance, concealment or decoys. If the period 
in which the reconnaissance is to take place is at hand, 
he will be more restricted, for example perhaps to 
moving a limited number of missiles from site to site 
in order to make empty sites decoys. Concealment 
will be fixed, since this involves large capital expense. 
The latter situation, equipment limitation, is a special 
case of the former budget limitation. A solution is 
given for the special but not the general case, though 
it is proved that there exists a solution in mixed strategies 
for the latter. 

A model is given for evaluating the effectiveness of a 
human photointerpreter, or a photo-interpreting system; 
this model separates the confusion due to the photo- 
interpreter (human or system) from that due to the 
photographs. The expected confusion to an interpreter 
using the system R is denoted by Cap, the confusion 
inherent in the photographic system used by Cp. 
Formulas are given for Cp and Cp; it is proved that 
CrSCp. The minimum achievable confusion is Cp: 


continued 


10.5 (0.7) 
continued 


ARL Technical Report 60-337 (1960), Aero. Res. Labs., Wright-Patterson Air Force Base. 


iv-+20 pp. (2 references) 


the measure of effectiveness of the photo-interpreting 
system is defined as E = Cr—Cp, the additional con- 
fusion introduced by that system over that inherent in 
the photographs. ; 

An error in WADC Technical Note 59-409, in the 
statement (Section 11, page 30) about the solution for 
the concave-convex case, is pointed out and corrected: 
the results are not invalidated. 


(H. L. Harter) 
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FISZ, M. (Math. Inst., Polish Acad. Sci., Warsaw) 


10.0 (-.-) 


Remarks on the sample functions of some stochastic processes—ZIn English 


Bull. Acad. Polon Sci. III (1960) 8, 355-358 


The author considers the real separable stochastic 
process {x,, t€1)}, where Jy is a closed finite interval. 
Denote further by x, the increment of x, in the interval / 
and by | /|—the length of 7. Consider. the following 
functions: 


a(I) = Pr(x, #0), 5, e) =Pr (x; | >8), 


A() = | a(1), BU, 8) = | b(, €), Q(t) =lima(l)/| 7| 
T F 

as J contracts to a fixed point te/. A(J/) and B(I, «) are 

the Burkill integrals of a(/) and b(/, «), respectively. 

The author characterises the sample functions of the 

process {x,, t € Jj} under different assumptions concern- 

ing above-mentioned functions. 


(L. Kubik) 
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GHOSAL, A. 
Emptiness in the finite dam—ZJn English 


(Central Fuel Research Institute, Dhanbad, India) 


10.4 (10.1) 


Ann. Math. Statist. (1960) 31, 803-808 (8 references) 


This paper discusses the general problem of emptiness 
in the finite dam and considers the probability that, 
starting with an arbitrary storage, the dam dries up 
before it fills completely. Both discrete and continuous 
inputs are considered. 

In the model considered, the storage Z, defined for 
few) ali, 2,0... 1S) CiVeTD by: 


Z,+X,-—m if m<Z,+X,<k 
Zs = <0 if Z,+X,<m 
k—m if Z,+X,2k 


where the YX, are identically distributed independent 
inputs, which flow into the dam in the intervals (t, t+ 1). 
Defining V; as the conditional probability that, starting 
with storage i, the dam becomes empty before it fills 
completely an exact solution for V; is obtained under 
the assumption of a geometric input, Pr (X, = j) = ab’; 
b=1-a, j=0,1,2,.... 
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The continuous analogue of V;, V(yv) is defined and 
exact solutions for V(y) are obtained under the assump- 
tion of an exponential input, g (X) = wexp(—-/X), 
X>0 and a gamma input, 


Bey ages 
g (X) ia 
It is shown that if H (y) is the stationary continuous 
distribution function of the dam content, then V (y) = 
H (k—y) for the continuous inputs considered. This 
relation also holds for a discrete input. 


Expi(= x) Xx), Ae 0: 


(P. L. Meyer) 


HAIGHT, F. A. (University of California, Los Angeles) 


Queueing with balking. II—Jn English 
Biometrika (1960) 47, 285-296 (8 references) 


In a previous paper [Biometrika (1957) 44, 360-369] 
the author considered the queueing problem where the 
rate of arrivals may be greater than the service rate. 
However, arrivals join the queue only if its length is K 
or less otherwise they balk, i.e. go away and do not return. 
The values of K chosen by successive arrivals are regarded 
as being random samples from the balking distribution. 
The inter-arrival time and service time distributions are 
negative exponential and there is a single server. 

In his paper the author considers queue lengths and 
waiting-time distributions under certain specific condi- 
tions and also the probability distribution for the time 
gap between an arrival (whether joining or not) and the 
previous joiner. The corresponding quantities for gaps 
between a joiner and the preceding joiner are also 
found. 
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IOSIFESCU, M. & THEODORESCU, R. (Central Statistical Office, and 


10.4 (2.5) 


Finally a numerical example is given based on a 
simulated queue using a Poisson balking rule. The 
observed numbers of negative exponential phases in 
inter-arrival and inter-joining gaps are given and geo- 
metric, negative-binomial and Borel distributions have 
been fitted. 


(C. Burrows) 


10.1 (1.0) 
Inst. Math., Bucharest) 


On certain linear chains with complete connections—I/n Russian 


Rey. Math. Pures Appl. (1961) 6, 167-170. 


The problem of the ergodic structure of the chains 
with complete connections was first considered by 
Fortet in his thesis for the alternative linear case. In 
the present paper the authors provide adequate condi- 
tions under which a linear chain with complete con- 
nections with a finite set of states, containing more than 
two states, can be considered as an alternative linear 
chain. For details on these processes see the recent 
monograph by Ciucu & Theodorescu [Processes with 
complete connections (1960) Bucharest: Ed. Acad.R.P.R. ]. 


(G. Simboan) 
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LAMPERTI, J. (Stanford University) 


The first-passage moments and the invariant measure of a Markoff chain—Jn English 
Ann. Math. Statist. (1960) 31, 515-517 (3 references) 


10,422) 


The author considers an irreducible recurrent Markoff 
chain and develops a simple expression for the (k+1)th 
factorial recurrence-time moment for state 0 in terms 
of the factorial first-passage time moments for all 
positive states and the invariant measure of the chain, 
Pvelid Vor k=0.015 20.0 This generalises the well- 
known formula for k = 0 and the formula for k = 1 
due to Chung [Trans. Amer. Math. Sec. (1954) 76, 


397-419]. 
(W. J. Hall) 
DIGS 
LESLIE, P. H. & GOWER, J. C. (Bur. Animal Pop., Oxford and Rothamsted Exp. Sta.) 10.9 (-.-) 


The properties of a stochastic model for the predator-prey type of 
interaction between two species—JIn English 


Biometrika (1960) 47, 219-234 (9 references, 4 tables, 2 figures) 


A model is considered in which the birth rates (per Carlo realisations of the process. Chances of extinction 
individual, per unit time) of two species, prey and _ of either species are also discussed. 
predator, are constant while the death rate of the prey A second model is then considered in which only a 
is approximately linear in both the numbers of prey proportion k of the prey population is exposed to risk. 
and predator and the death rate of the predator is This model is treated similarly to the first. 
approximately linear in the ratio of predator to prey. 
The model is in fact stated, for computational purposes, (A. S. Miller) 
in terms of the means and variances of the numbers of 
each species (assumed to be normally distributed before 
rounding off to the nearest non-negative integers) present 
at time ¢+1, given the numbers of each species at time f. 

Approximate theoretical variances and covariances 
for the fluctuations of the process in the neighbourhood 
of the stationary state are given and are shown to 
compare quite favourably with the results from Monte 
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LUKASZEWICZ, J. (Wroclaw University, Wroclaw) 10.4 (10.1) 
A remark on a theorem by Khintchine from the theory of random streams—Jn Russian 
Coll. Mat. (1960) 7, 285-287 


The author notices that a lemma in Khintchine’s book 
[Mathematical Methods in the T, heory of Queueing (1960) 
London: Griffin] stating that for every Palm stream, 
that is for every stationary orderly stream with a limited 
after-effect, and for every r = 1, 2, ..., we have 


ne Wr+i (Ui, (u) =0, 


where W, (u) is the probability of having at least r calls 
in a time segment of length u, is true only under the 
additional assumption that yw, (u)>0 for u>0. 
Moreover, he gives a generalised version of this 
lemma which is true for all Palm streams without 
exception and enables us to complete the proof of a 
theorem about probability distribution of distances 
between consecutive calls in Palm streams in its full 
generality; as it is stated in Khintchine’s book. 


7 (S. Zubrzycki) 
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MEERENDONK, H. W. van der & BOS, G. G. J. (N.V. Philips, Eindhoven) 10.4 (-.-) 
Rules for the classification of seasonal production—In Dutch 
Statist. Neerlandica (1961) 15, 177-187 (2 tables, 1 graph) 


This paper deals with the analysis of a problem often 
met when determining purchasing policies. From a 
number of different kinds of goods, which are to be 
bought, the total quantity is restricted by storing capacity 
or financial limits. The problem “‘ how much of each 
kind will be bought ”’ is solved for a discrete as well as 

_ for a continuous distribution of demand for the various 
goods. 


(H. W. van der Meerendonk) 


NELSON, E. & VARBERG, D. (University of North Dakota & Hamline University, Minn.) 


10.1 (1.2) 


Expectations of functionals on a stochastic process—In English 
Ann. Math. Statist. (1960) 31, 574-578 (6 references) 


Let Xt}, 0St< 00 be the separable and centred 
infinitely divisible stochastic process which gives the 
sum of N(t) independent Bernoulli random variables 


_ with equally likely values +1; where N (t) is an in- 


* 


dependent Poisson process. The characteristic function 


of this process is then 


& {exp [iéx {t}]} = exp [at (cos €—1)]. 
Let V be a non-negative real function and set 


W (u) = | “ee ex| —u | Vix(o} a m| dt 


where /(?, n) equals 1 or 0 according as DL Ot 
xn. 

The main theorem of this paper states that Y,, satisfies 
the difference equation . 

Pry —(2/a)[statuV (n)] Areas i a —(2/a)d,. 0» 


the quantity 6, 9 being the kronecker delta. The 
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PARTHASARATHY, K. R. (Indian Statistical Institute, Calcutta) 


method and notation follow those of Kac “On some 
connections between probability theory and differential 
and integral equations” [Proc. Second Berkeley Sym- 
posium Math. Statist. & Probability (1951) 189-215]. 

Two examples are then given, firstly V(x) is zero 
when a<x<b and is unity otherwise; and secondly 
Ve) x22 


(R. Pyke) 


10.6 (10.2) 


On the estimation of the spectrum of a stationary stochastic process—In English 
_ Ann. Math. Statist. (1960) 31, 568-573 (7 references) 


Let x1, Xz, ..., Xy be N observations of a discrete para- 
meter weakly stationary stochastic process; a second 
order process with stationary mean and covariance 
functions. setic. x, = O'and p= 6 x,x,4,,-., Let 


Ray nyt tT ON | oI xy)/(N—| n » 
be the sample covariances, and consider the following 
estimates of the spectral distribution F: 


R(N) 


DM Cl) is Ok, ema a 


k = —R(N) 
where the summand for k = 0 is interpreted to be 
Ay, y (A+n) and the a,, y are constants chosen to satisfy 
firstly, a, y>1 as Noo, secondly ay = 4 ,,y and 
thirdly, Fy is a distribution function. 

Investigations into the asymptotic behaviour of such 
‘weighted periodogram” estimates of the spectral 
density function have been carried out under various 
assumptions by Grenander & Rosenblatt [Statistical 
Analyses of Stationary Time Series (1957) New York: 
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Wiley] and by Parzen, [see, for example, ““On asymp- 
totically efficient consistent estimates of the spectral 
density function of stationary time series,” J. R. Statist. 
Soc. B (1958) 20, 303-322: abstracted in this journal 
No. 1/138, 10.6]. 

“In this paper, the problem of estimating the spectral 
distribution as well as the spectral density, if it exists, 
of a weakly stationary process is solved under the sole 
assumption that the sample covariances converge 
almost surely and in mean to the true covariances. 
. . . The existence of estimates which converge uniformly 
strongly to the spectral density of the process is proved 
under the assumption that the density has an absolutely 
convergent Fourier series,” states the author. He 
applies a theorem of Bochner to obtain a uniform bound 
on the bias of these estimates in the case of a continuous 
spectral distribution. 


(R. Pyke) 
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RAO, M. M. (Carnegie Institute of Technology, Pittsburg) 


Estimation by periodogram—ZIn English 
Trab. Estadist. (1960) 11, 123-137 (7 references) 


This paper is devoted to establishing the relationship 
between the methods of estimation of the period and 
amplitude of a harmonic model by periodogram and 
maximum-likelihood or least-squares. 

Let T be a set of integers and x,, te T, be a process 
which satisfies the relation x, = acos wt+b sin ott &,, 
where a, 6 and o are real constants to be estimated 
from the observations on x,, and &p- are ““errors ””: 
Assume: 

(i) —p<4a, 


b<p, 0<w<p. 


(ii) The random variables ¢,, t¢7, have zero mean’ 


and are uncorrelated with a finite positive 
variance independent of f. 
(iii) T<J, where J is a finite fixed interval. 
let X,(t= 1, ..., n) be the observations on the 
process and define 4, (0), 5,(@) and c,(@). Then the 
Schuster method is to choose that value of 0 as an 
estimator @ of w for which | c, (0)|* is a maximum 
and the corresponding values of d, (0) and 6, (0) as 
the estimators of a and b in (i). 

Two sequences {Y,} and {Z,' of random variables 
are said to be stochastically equivalent if for any e>0 
hm Pe Y= Z| > 2) =. 0. 
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THEODORESCU, R. (Institute of Mathematics, Bucharest) 
On chains with complete connections—/n Roumanian 


10.2 (-.-) 


Let ¢, be normal and suppose the assumptions (i) and 
(ii) hold. Then the periodogram estimators d, (0), 6, (8) 
and (0) = @, of the Schuster method are stochastically 
equivalent to the maximum-likelihood estimators 4, (0) 
5, (0) and @,,. 

The above proposition holds if maximum-likelihood 
is replaced by least-squares when normality is dropped 
and the independent random variables ¢, are assumed 
to have two finite moments with continuous distributions. 

A matrix is given which provides an approximation, 
or more precisely the lower bound, for the matrix of 
the moments of the maximum likelihood method; it 


is also an approximation for n large for the matrix of 


the moments of d, 6 and @ given by the Schuster method. 

The lower Cramér-Rao bound for the variance of 
the estimator of the period is given and lastly the 
following is proved: if (i) to (iii) hold and the e, are 
normal with mean zero and finite variance and the 
parameters (w, a, b, o”) belong to a bounded non- 
degenerate set of the 4-dimensional Euclidean space, 
then the maximum likelihood estimators, and hence 
also the Schuster estimators, are consistent. 


(R. Pro) 


10.1 (1.0) 


Studii cercetari matematica (Cluj) (1960) 11, 195-215 (49 references) 


The notion of a chain with complete connections called 
“chains of infinite order” by other authors, has been 
introduced into probability theory by Onicescu & 
Mihoc in their work “Sur les chaines de variables 
statistiques” [C. R. Acad. Sci. Paris (1935) 200, 511- 
512]. This notion, which generalises that of a Markoff 
chain, has been later extended by Doeblin and Fortet; 
the systematic study of these chains is due especially to 
the Roumanian school of probabilists. 

The material concerning the theory of chains with 


complete connections may be found scattered among © 


an entire series of memoirs and monographs. The 
aim of the present expository paper is to put into 
perspective the central results of this theory, emphasising 
the Roumanian contributions [for details, see the 
recent monograph ‘Processes with complete con- 
nections” by Ciucu & Theodorescu. See also abstract 


No. 2/746]. 
(M. Tosifescu) 
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THEODORESCU, R. (Institute of Mathematics, Budapest) 
An equation occurring in the study of stochastic processes—In Russian 
Rey. Math. Pures Appl. (1960) 5, 719-722 (4 references) 


10.1 (1.0) 


In this short paper a special functional equation occurring 

in the theory of processes with complete connections 

is studied. . , 
Under adequate conditions the solution exists, is 

unique and represents a probability. For details on 

these processes, see the recent monograph ‘“‘ Processes 

with complete connections” by Ciucu & Theodorescu. 

See also abstract No. 2/746. 


(M. Tosifescu) 
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THOMASIAN, A. J. (University of California, Berkeley) 10.5 (-.-) 
An elementary proof of the AEP of information theory—IJn English 
Ann. Math. Statist. (1960) 31, 452-456 (5 references) 


The author ‘develops some limit relationships among 
certain characteristics of the sequence of random 
variables—(1/n) log p for an arbitrary, not necessarily 
ergodic or stationary, information source. These permit 
an elementary combinatorial proof of the asymptotic 
equipartition property that the sequence from a station- 
ary ergodic source converges in probability to a constant. 
This was introduced and first proved by McMillan, 
“The basic theorems of information theory” [Ann. 
Math. Statist. (1953) 24, 196-219]. 


(W. J. Hall) 
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WALKER, A. M. (Statistical Laboratory, University of Cambridge) 
Some consequences of superimposed error in time series analysis—In English 


10.7 (10.2) 


Biometrika (1960) 47, 33-43 (7 references, 2 tables) 


In this paper a study is made of some of the problems 
that occur when there is a superimposed error in a time 
series. It is assumed that the errors are additive and 
form an independent stationary series which is in- 
dependent of the original time series. The observations 
are equally spaced and consist of the sum of a stationary 
linear autoregressive process of order p and a completely 
random series. 

Assuming that each of these two processes, and 
therefore the sum, are normal then the behaviour of 
the observed series is completely determined by its 
autocovariance function. It is easily shown that the 
effect of the superimposed error is to reduce all the 
autocorrelations in the same ratio. 

Two methods of estimating the parameters of the 
sequence are considered. The first (a) uses the first 
p+2 serial covariance while the second (b) uses a 
different set of p+2 consecutive serial covariances. 
The questions posed are: 


(i) What is the best set of covariances to be used 
in the second method judging by the (asymptotic) 
efficiency of the estimates? 
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WORKING, H. (Stanford University, California) 


(ii) How does the efficiency of the second method 
compare with the first? 

(iii) How far are the previous answers affected by 
non-normality of the processes? 


The calculations are very tedious and a detailed 
examination is given only for the case p= 1. Asymptotic 
efficiencies of the estimates are found and tabulated. 
Finally the effect for non-normality is considered when 
it is shown that the non-normality of the error term 
produces an effect which can, in principle, be corrected 
but leads to difficult algebra even for the case p = 1. 


(C. Burrows) 


10.7 (6.4) 


Note on the correlation of first differences of averages in a random chain—In English 


Econometrica (1960) 28, 916-918 


In this note the author, after pointing out the importance 
of bearing in mind the fact that with serial correlation 
the use of averages can introduce correlations which 

were not present in the original series, considers the 
effect of averaging successive groups of items contained 
inarandom chain. These are the most primitive type of 
stochastic series which have a close resemblance to 
either stock prices or to certain commodity prices. 
The author uses an example from his paper “ A random- 
difference series for use in the analysis of time series ” 
[J. Amer. Statist. Ass. (1934) 29] to show that with the 
number of items in each successive segment only 
moderately large, the variance of the first differences 
between averages over successive segments of a random 
chain is approximately two-thirds of the variance of 
the first differences between correspondingly positioned 
terms in the chain. 
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When the number of items is fairly small, the expected 
first-order serial correlation of the ‘first differences 
between average terms in the random chain approximates 
to &(r,) = +4. Serial correlation of order higher 
than the first remain zero for first differences of averages 
of successive groups of terms in the random chain. 


(W. R. Buckland) 
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ZUBRZYCKI, S. (Math. Inst., Polish Academy of Sciences, Wroclaw) 


Concerning plane sampling—ZJn English 


10.2 (8.1) 


Second Hungarian Mathematical Congress (1961) (6 references) 


Let y(p) be a plane stationary continuous stochastic 
process depending on a point p in a plane. How to 
choose n points in a given domain D in order to minimise 
the mean-square error of estimating the mean value 
n (D) of the process y(p) over D by the average of the 
values of the process y (p) observed in the chosen points 
is the question posed by the author. It is the correlation 
function of the process y(p) which determines the answer 
to that question. In practice, however, the correlation 
function is usually not known, and even if we knew 
it, the computations would be prohibitive. This suggests 
searching for optimal methods of sampling in some 
suitably narrowed class of sampling methods with the 
advantage of having optimal solution for a possibly large 
class of processes yp. 

This paper is a report on some recent investigations in 
this field: see also Zubrzycki [“‘ Remarks on random 
stratified and systematic sampling in a plane,” Colloq. 
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Math. (1958) 6, 251-264]; Dalenius-Hajek-Zubrzycki 
[On plane sampling and some related geometrical 
problems,” in the press] and Ludmilla Zubrzycka [‘‘ On 
the distribution of sampling points in a plane’, Zasto- 
sowania Mat. (1960) 5, 161-171: abstracted in this 
journal No. 2/719, 8.0]. 


(S. Zubrzycki) 
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BANKOVI, G. (Math. Inst., Hungarian Academy of Sciences, Budapest) 


11.7 (1.4) 


Evaluation of integrals by Monte Carlo methods based on the one-dimensional 


random space filling—Jn English 


Publ. Math. Inst. Hung. Acad. Sci. (1960) 5, 339-352 (8 references) 


The paper deals with the approximative evaluation of 


integrals of special type 


I {g} = i exp | g (u) du de 


by sampling methods. This procedure is based on the 
one dimensional random space-filling problem solved 
by Rényi [Publ. Math. Inst. Hung. Acad. Sci. (1958) 
3, 109-127: abstracted in this journal No. 1/18, 1.4]. 
This consists of placing at random disjoint intervals 
on a longer interval, and determining the mean value 
of the number of intervals which can be placed. It leads 
to a similar type of integral. 

The author generalised this problem by using intervals 
whose length is variable, weighting these and changing 
the placing procedures. 


(Ilona Palasti) 
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BAZLEY, N. W. & DAVIS, P. J. (National Bureau of Standards, Washington, D.C.) 


11.7 (10.1) 


Accuracy of Monte Carlo methods in computing finite Markoff chains—Jn English 
J. Res., Nat. Bur. Stand. B (1960) 64, 211-215 (5 references, 3 tables, 2 figures) 


Since Monte Carlo methods are generally applied to 
problems that are too difficult for other types of analysis 
there is seldom opportunity to compare their results 
with ‘‘ exact solutions”’. This paper reports results of 
an experiment performed on the children’s game “‘ Chutes 
and Ladders”, a game which can be interpreted as a 
finite Markoff chain with 82 states and with given 
transition probabilities. Statistics, such as the average 
length of play, are computed on the IBM 704 from 2'* 
simulated plays of the game. These Monte Carlo results 
are then compared with the “ exact solution ” obtained 
by powering the matrix of transition probabilities. 

The random element in the game of Chutes and 
Ladders arises from the toss of a dice for each move. 
This was simulated by employing the pseudo-random 
number generator of Taussky & Todd. The unit 
interval was divided into six equal subintervals and 
the result of the jth dice throw was identified with the 
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sub-interval in which the pseudo-random number f; 
fell. The average length of play after 16,384 games was 
39,420 throws. The exact answer is 39,224. 

The author concluded that the results of Monte 
Carlo calculations are in agreement with the theoretical 
results obtained from the matrix calculations and the 
convergence according to the 1 Ix N law. 

Note: In equation (5.2) “a” should read “ o? 
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(H. H. Ku) 
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DEMING, Lola S. (Statist. Eng. Lab., National Bureau of Standards, Washington) 


11.9 (6.0) 


Selected bibliography of statistical literature, 1930 to 1957. 


I. Correlation and regression theory—In English 
J. Res., Nat. Bur. Stand. B (1960) 64, 55-68 


This is the first in a series of bibliographies dealing with 
various specific subjects in the field of statistics. The 
series is intended to be of service by listing the most 
important contributions in statistical theory and method, 
as judged by two prominent reviewing journals [Zentra- 
blatt fiir Mathematik for the years (1930) to (1939), 
and Mathematical Reviews from (1940) onward]. The 
scheme of classification follows that of the Mathematical 
Reviews Annual Index. 

The titles are given exactly as in the reviewing journal, 
and are listed under the first-named author only in the 
case of joint authorship. Joint authorship is denoted 
by a symbol preceding the surname and the authors 
are cross-referenced to the first-named author. 

In addition to the name of publication, volume, 
initial page number, and date of publication, the volume 
and page numbers are given for the two reviewing 
journals in which the abstract appears, 
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DEMING, Lola S. (Statist. Eng. Lab., National Bureau of Standards, Washington) 


There are approximately 520 entries in this biblio- 
graphy on correlation and regression theory. See also 
Il. “‘ Time Series ’’; abstracted in this journal No. 2/772, 
11.9: III. “ Limit Theorems”; abstracted in this journal 
No. 2/773, 11.9 and IV. ‘‘ Markoff chains and stochastic 
processes ”’; abstracted in this journal No. 2/774, 11.9. 


(H: Bako 


11.9 (10.0) 


Selected Bibliography of statistical literature, 1930 to 1957. 


II. Time Series—Jn English 
J. Res., Nat. Bur. Stand. B (1960) 64, 69-76 


This is the second in a series of bibliographies dealing 
with various specific subjects in the field of statistics. 
References and titles of important contributions to the 
study of time series have been taken from a wide variety 
of technical journals published in the many languages 
and countries which have been actively engaged in 
statistical analysis. For a description of the bibliography 
in general see I. ‘‘ Correlation and Regression Theory ”, 
abstracted in this journal No. 2/771, 11.9. 

There are approximately 280 entries in this bibliography 
on time series. See also III. ‘‘ Limit Theorems ’”’; 
abstracted in this journal No. 2/773, 11.9 and IV. 
“ Markoff Chains and Stochastic Processes’’; abstracted 
in this journal No. 2/774, 11.9. 


(H. H. Ku) 
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DEMING, Lola S. (Statist. Eng. Lab., National Bureau of Standards, Washington) 11.9 (1.5) 


Selected Bibliography of statistical literature, 1930 to 1957. 
II. Limit theorems—-Jn English 


J. Res., Nat. Bur. Stand. B (1960) 64, 175-192 


This is the third in a series of bibliographies dealing 
with various specific subjects in the field of statistics. 
References and titles of important contributions con- 
cerning limiting distributions have been taken from 
technical journals published throughout the world since 
1930. 

There are approximately 670 entries in this biblio- 
graphy. The two preceding ones are: I. “Correlation and 
Regression Theory”; abstracted in this journal No, 
2/771, 11.9: II. “ Time Series”: abstracted in this 
journal No. 2/772, 11.9 and there is also IV. “ Markoff 
Chains and Stochastic Processes’ ; abstracted in this 
journal No. 2/774, 11.9. 


(H. H. Ku) 
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DEMING, Lola S. & GUPTA, D. (Statist. Eng. Lab., Nat. Bur. Stand. _ 11.9 (10.0) 
& Catholic Univ., Washington) 


Selected Bibliography of statistical literature, 1930 to 1957. 
IV. Markoff chains and stochastic processes—Jn English 
J. Res., Nat. Bur. Stand. B (1961) 65, 61-93 


This is the fourth in a series of bibliographies dealing 
with various specific subjects in the field of statistics. 
Titles and references of important contributions to the | 
study of Markoff chains and stochastic processes have 
been taken from technical journals published since 1930 
in the many countries that have been actively engaged 
in statistical analysis. | 

There are approximately 1250 entries in this biblio- | 
graphy. See also I. “Correlation and Regression 
Theory ”; abstracted in this journal No. 2/771, 11.9: IL. 
“* Time Series”’; abstracted in this journal No. 2/772, 11.9 
and III. ‘‘ Limit Theorems ’’; abstracted in this journal 
No. 2/773, 11.9. 


(H. H. Ku) 
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EVANS, Marilyn J. & HILL, G. W. (Div. Math. Statist., CSIRO, Adelaide & Melbourne) 
Tables of z = arc tanh (r)—Jn English 


Diy. Math. Statist. Tech. Paper No. 7 (1961) 13 pages (2 references, 5 tables) 


11.1 (3.1) 


This technical paper supplements Table VII of Fisher 
& Yates [Statistical Tables for Biological, Agricultural 
and Medical Research, 4th Edn. (1953) Edinburgh: 


Oliver & Boyd]. Tables are given for the function 
z = arc tanh (r) = $ {In (1+r)—In(1—-nr)} 


required to obtain Fisher’s z-statistic from known values 


Ola 


There is a brief introduction to the tables, followed 


by: 


Table 1 (a): 
Table 1 (5): 
Table 2 (a): 
Table 2 (d): 
Table 2 (c): 


values of z forr = 0 (0-001) 1, 
values of z forr = 0-9 (0-0001) 1, 
values of z for r? = 0 (0-0001) -05, 
values of z for r? = 0 (0-001) 1, 
values of z for r? = 0-9 (0-0001) 1. 


(J. Gani) 
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HAIGHT, F. A. (University of California, Berkeley) 


11.9 (2.0) 


Index to the distributions of mathematical statistics—Jn English 


J. Res., Nat. Bur. Stand. B (1961) 65, 23-60 


A fairly complete index of references to results on 
statistical distributions published before January 1958 
is presented. The distributions covered are characterised 
as normal, type III, binomial, discrete, distributions 
over (a, b), distributions over (a, 0), distributions over 
(—co, ©), miscellaneous univariate, miscellaneous 
bivariate, and miscellaneous multivariate. The number 
of entries varies from one or two for less well-known 
distributions to several hundred for the normal distribu- 
tions. This index should serve to eliminate unnecessary 
derivation of results already in the literature. 

The material given under each distribution consists 
of a number of entries, organised under a standard 
order of headings for well-known distributions with 
sufficient number of entries, for example, the normal: 


A. Functions and Parameters, 
B. Derived distributions. 

C. Estimation. 

D. Testing statistical hypotheses. 
E. Miscellaneous. 


Direct and indirect references are used in this index: 
direct references include twenty-six coded journals and 


eighteen coded books, and uncoded literature. The 
indirect references give reviews by Mathematical Reviews 
and Zentralblatt fiir Mathematik—usually for papers 
that appear in obscure sources. A distribution is 
included if its density, or probability function is a 
known, explicit function, and an entry must exhibit a 
property of the distribution in question. Exceptions 
include certain densities specified only implicitly or in 
terms of their cumulative probability function, or 
characteristic function. Historical information, impor- 
tant applications, and bibliographies are sometimes given. 
An alphabetical index to the index to distributions is 
also included. 

An example of the entries is: 

Wilk’s distribution of dispersion determinant, etc.: 
(e) 3:26 where the reference given is Sankhya, Vol. 3, 
D2; 


(H. H. Ku) 
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NAGAR, A. D. (Netherlands School of Economics, Rotterdam) 


11.7 (6.6) 


A Monte Carlo study of alternative simultaneous equation estimators—In English 
Econometrica (1960) 28, 573-590 (4 references, 8 tables) 


In this paper the author discusses a Monte Carlo 
approach to four methods of estimating the small 
sample properties of parameters of simultaneous 
equations. The methods discussed are: 


(1) least squares 
(ii) two-stage least squares 
(iii) unbiased second-moment 
(iv) minimum second-moment. » 


After a reference to a paper by Wagner “ A Monte 
Carlo study of estimates of simultaneous linear structural 
equations” [Econometrica (1958) 26, 117-133], the 
author states that his paper is a further examination 
of Wagner’s models. This examination is conducted 
in the light of the alternative procedures now available 
in order to carry out a comparison. The present author 
is able to disregard the methods of estimation already 
discussed by Wagner. After stating that the two 
equations used by the previous author are “‘ extreme ”’, 
in a certain sense, the present author hopes that a 
further examination may give valuable information 
regarding these “‘ extreme ’”’ type equations. This paper 
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THEIL, H. (Netherlands School of Economics and Stanford Univ.) 


considers both the over-identified and just-identified 
equation in the model. 

The second section deals with the alternative estimation 
procedures; Wagner himself having considered the 
limited-information estimates of the over-identified 
equations. The procedures and their properties are 
described but not proved. 

After a section which gives further consideration to 
Wagner’s models the author gives numerical results 
for two models in terms of one hundred samples of size 
twenty. The four methods of estimation are considered, 
all belonging to the class, and the final section of this 
paper gives results which include observations on 
sampling variance, bias and asymptotic standard error. 
The author gives eight tables of various points of com- 
parison between the methods. 


(W. R. Buckland) 


11.8 (6.4) 


Best linear index numbers of prices and quantities—Jn English 
Econometrica (1960) 25, 464-480 (5 references, 1 figure) 


This paper deals with the statistical type of approach 
to the subject of index numbers. It differs from the 
usual approach in as much as it is not necessarily only 
concerned with one single pair of periods. The author 
considers index numbers for an arbitrary number of 
periods which are determined by the prices and quantities 
of all separate periods. This procedure is of course 
necessarily ex post as it can only be used when the 
prices and quantities for the whole year are available; 
it is also suitable for computation of prices and quantity 
indices needed for an econometric model. Another 
situation to which it is applicable is the comparison of 
prices and quantities between different geographical 
units and the symbol ¢ in the matrix used in the section 
on best linear price index numbers can also be used to 
denote either time units or geographical units, although 
the exposition is largely confined to variations over 
time. 

The primary sections of this paper deal with the 
geometric and algebraic aspects of this approach which 
is a “linear”? one in the sense that the price index 


vector is a linear function of the matrix of the individual 
prices. It is closely related to principal component 
analysis from which it differs by virtue of the fact that 
it considers two sets of variables, both prices and 
quantities. Sections two and three deal with prices only 
and are followed by a completely analogous approach 
to quantities in section four. 

In the fifth section price and quantities are combined 
and the following section deals with a filtering index 
which serves to measure the statistical quality of the 
price and quantity index vectors. Partial index numbers 
and problems of aggregation follow in the next section 
and in conclusion the author comments on low filtering 
indexes and the weighting of procedures to deal with 
years which are considered “ outlying ”’. 


(W. R. Buckland) 
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THEIL, H. (Netherlands School of Economics, Rotterdam) 
Best linear index numbers of prices and quantities—In English 


11.8 (3.4) 


Econometrica (1960) 28, 464-480 (5 references, 1 figure) 


This paper is statistical in that it is concerned with 
the optimum specification of a central tendency for 
price and quantity ratios. Its approach is different, 
however, in that it is not necessarily concerned with 
one pair of periods but with an arbitrary number of 
complete periods. This approach is suitable for the 
estimation of an econometric model and can be applied 
on a current basis for price and quantity comparisons 
of different geographical units. 

In the second section the approach of the paper is 
demonstrated geometrically followed by an algebraic 
account in section three. The approach is linear in the 
sense that the price index is a linear function of the 
matrix of individual prices. Because a quadratic form 
relating to certain discrepancies is minimised the method 
of this paper is a “ best’ approach and closely related 
to principal component analysis. While sections two 
and three deal only with prices, the analogous situation 
for quantities is dealt with in section four. The two 
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VENEKAMP, P. E. (University of Amsterdam) 
Descriptive statistics—In Dutch 
Statist. Neerlandica (1961) 15, 171-174 


In this article the concepts of descriptive statistics and 
theoretical statistics are compared. It is concluded 
that descriptive statistics are to be considered an applied 
science, aimed at presentation and interpretation of 
figures. In this context the methods of theoretical 
statistics are also to be applied. 


(P. E. Venekamp) 
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situations are combined in sections five and six. The 
next section deals with a “ fitting-index ” which measures 
the statistical quality of the price and quantity vectors. 
The problems of partial index numbers and problems 
of aggregation are dealt with in section eight. 

The final section contributes some general comments 
by way of interpretation. For example, if the “ fitting- 
index ”’ is low then the second largest latent root of the 
price and quantity vectors is used: both being pairwise. 
If any particular year is an “‘ outlying” observation it 
is possible to apply a weighting procedure resembling 
Aitken’s method of generalised least squares [Proc. 
Roy. Soc. Edinb. (1934-35) 55, 42-48]. 


(W. R. Buckland) 
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